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FOREWORD 


This document is part of a three volume report prepared under NASA Ames 
Contract NAS2-10294^ Mathematical Modeling for Vertical Attitude Take-Off 
and Landing (VATOL) Simulation* Volume I: Model Description and implication 

provides background and details of a generic mathematical model for simulation 
of VATOL aircraft conc^ts* A six-degree-of-freedom off-line (non-piloted) 
digital simulation program incorporating this model was developed and applied 
to the Vought SF-121 VATOL concept. Volume I gives results of this application 
which included development and demonstration of a control system for terminal 
VATOL operations. Volume II: Model Equations and Base Aircraft Data gives 

all the model equations and SF-121 aircraft data in a simulation data package 
format. This volume facilitated the development of a piloted VATOL simu- 
lation at NASA Ames. Volume III: Users Manual for VATOL Simulation Program 

provides a description of the six-degree-of-freedom off-line digital simulation 
program, instructions for Its application, and examples of setup decks and out- 
put for several of the SF-121 application runs. 

Project Monitor for NASA Ames was Mr. Gary Hill. The Princxpal Investi- 
gator for the Vought Corporation was Robert L. Fortenbuagh. 
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1 . 0 INTRODUCTKM 


The objective of this document is to provide a description of and 
instructions for using VATLAS (Vertical Attitude Takeoff and lading Air- 
craft Simulation)* the digital simulation program for application to 
vertical attitude takeoff and landing (VATOL) aircraft developed for instal- 
lation on the NASA Ames CDC 7600 computer system* Details of the VATLAS 
mathematical models are provided in Volumes I and II of the contract 
documentation and will not be discussed here. The framework for VATLAS is 
supplied by Vought's OLSIM (Mf-Line Simu lation) routine depicted in figure 
1 - 1 . 


OLSIM is a multi-purpose digital computer program which is used for 
stability and control analysis and flight control system design analysis 
in addition to off-line flight simulation. It provides a flexible framework 
and standardized modules which facilitate the developement of off-line 
aircraft simulations. OLSIM runs under the control of VTOLTH, the main 
program* which calls the proper modules for executing user-specified options. 
These options include trim* stability derivative calci: .^tion, time history 
generation, and various input-output options. Standardized modules include 
TRIM which has a six DOF nonlinear trim capability, RUNCX which performs 4th 
order Runge-Kutta integration* and various table lookup and matrix computation 
routines. User supplied aircraft specific modules include FORCES for 
calculating forces and moments on the airplane* SYSEQS for control system 
equations and aircraft kinematics* and lERIVS which perturbs the aircraft 
about trim to generate stability derivatives. VATIAS is the OLSIM routine 
with the FORCES* SYSEQS* and DERIVS modules specialized to VATOL aircraft. 

The aerodynamic model programmed in VATLAS is a modified version of 
the aerodynamic model described in reference (a) . It calculates aerodynamic 
force and moment coefficients using DATCOM-type relations for all angles of 
attack and sideslip. To review the aerodynamic coefficients generated by 
this model, a special wind tunnel option is included in VATLAS. This option 
provides for the calculation of nondlmens iogal aerodynamic coefficients as 
functions of angle of attack and sideslip. 

This users manual continues in Section 2.0 with brief descriptions of 
VATLAS program modules. Section 3*0 provides a discussion of COMMON array 
segmenting which is the key to both understanding VATLAS input-output options 
and developing other aircraft simulations with the Vought OLSIM framework. 

The basic VATLAS program load and program options are detailed in Section 4.0. 
Examples of four commonly used VATLAS run setups are provided in Section 5.0. 
The examples cover a multi-trim run* a time history run, a wind tunnel run* 
and a plot run. i^pendix A contains a definition of COMMON parameters (by 
convention* the phase "COMMON parameters" encompasses program variables and 
constants) arranged in order by segment* while Appendix B contains listings 
of all VATLAS subroutines. Note that the program has been heavily commented 
to facilitate location of specific groups of equations or options. 

One final note of advice to the prospective user: Because of its 

flexibility* VATLAS is not a simple program to learn* but once mastered it 
provides the familiar user with a powerful analytical tool. Careful study 
of the examples in Section 5.0 should ease this learning process. 
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Figure 1-1. Framework of OLSIM Routine 
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2.0 VATLAS PROGRAM MODULES 

The functlone of the verloue VATLAS progrea aodulee are described in 
this section. VATLAS exists as an UPDATE* source file. This facilitates 
user Interaction to modify VATLAS modules (l.e., ’’decks*' In UPDATE terminol- 
ogy by the use of UPDATE directives. A summary of VATLAS decks, their type 
and purpose, and their relative frequency of modification, as well as reasons 
for modification, are given In Table 2-1. Note that many of the decks are 
designated as "never" modified for VATLAS application. Users not knovledge- 
cble In the application of UPDATE or equivalent utilities are referred to the 
UPDATE Reference Manual (Reference (b)). 

2.1 The Main Program 

VTOLIH is the main program module for VATLAS. It controls most of 
the input and output and all the calculations of the program based upca 
sequences of program options selected by the user from a group of twelve 
which are currently available. These options provide for loading or reading 
data in the COMMON array, specifying trim variables or variables for 
stability derivative calculations, directing a trim solution or a time history 
calculation, storing a list of variables for input or output, directing 
wind tunnel-type calculations with the aerodynamic model, and storing and 
printing summaries of a series of trim solutions. Specific loading instructions 
for each of these options and the Basic Data Load (BDL) are given in Section 
4.0. 


2.2 The Aircraft-Specific Program Modules 

The VATOL aircraft-specific program modules Include IKRIVS, ERROR, 
FORCES, SYSEQS, LIFDRG, WTRUN, and NAMES. Should the user wish to apply 
the Vought OLSIM structure to an aircraft other than a VATOL, he would have 
to develop, as a minimum, the DERIVS, ERROR, FORCES, SYSEQS, and NAMES 
modules. LIFDRG and WTRUN were specially developed for VATLAS. 

2.2.1 FORCES 


Subroutine FORCES calculates the forces and moments applied to a 
VATOL aircraft by aerodynamics, propulsion system, inlet ram, Coriolis effects, 
and reaction control system (RCS) . There is no reasv>n to modify existing 
FORCES coding for application to any current VATOL concept. Future modifi- 
cations might Include the addition of a capability to calculate aerodynamic 
forces and moments from tabularlzed data as an alternative to the current 
DATCOM-type calculations and the addition of propulsion-induced aerodynamics. 

The FORCES aerodynamics model calculates contrlliucions for an aircraft 
with three lifting surfaces - wing, horizontal stabilizing surface, and 
vertical stabilizing surface - and a fuselage. Each lifting surface can 
have trailing and leading edge flaps whose effects are modeled with changes 
in CjQ, C^, and Cq per unit flap deflection and/or control surface whose 


*UPDATE is the CDC utility which creates an easily referenceable source program 
and provides easily applied directives for modifying this program with a 
minimum of input. Equivalent utilities are available for other computer 
systems . 
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Table 2~1 (Completed) Summary of Decks in VATLAS UPDATE Source File 
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effects are modeled with a change in angle of attack of the lifting surface 
per unit control deflection. The FORCES propulsion system can have one or 
two jet engines. IDach engine has an inlet, thrust vectoring in tw directions, 
and an independent throttle control. The FORCES RCS can have up to ten in- 
dividually controlled jets located anywhere on the aircraft and oriented in 
any direction. Each jet can be specified as demand or continuous bleed. 

FORCES also mv>dels the RCS-propulsion system interactions (figure 2-1) ap- 
propriate to continuous and demand bleeding. The FORCES inlet ram model 
provides for changes in inlet ram magnitude, direction, and application point 
as functions of angles of attack and sideslip, airspeed, inlet mass flow 
rate and aircraft rotation rates. The FORCES Coriolis model provides for 
changes in forces and moments as functions of mass flow rate through the 
engine ducts and aircraft rotation rates 

2.2.2 SYSEQS 


Subroutine SYSEQS contains the equations for the aircraft kinematics, 
flight control system, actuation system and pseudo-pilot functions. SYSEQS 
will require modification more frequently than any other subroutine primarily 
because of the great variety of possible control systems and control system/ 
actuation system interfaces. The user should not shy away from exercising 
his 1' dividuality and modifying SYSEQS to his requirements; VATLAS was designed 
to be easily modified. 

The aircraft kinematics in SYSEQS use direction cosines to orient 
aircraft body axes to inertial space; this avoids the singularity at 6 * 90 
deg i.i the standard Euler rate equations. The Euler rate equations are 
available and integrated in SYSEQS for application as required but are not 
fed into the kinematics. The integration of the Euler roll and yaw rate equa- 
tions is suspended when 6 is in the range 90 + 0.00006 deg; this suspension 
introduces an error each time the appropriate range is encountered. 

The flight control system provided in SYSEQS Includes the generic roll, 
pitch, yaw and heave control systems depicted in figures 2-2 through 2-5. Feed- 
back possibilities for moment control Include stability and body axis angular 
rates, Euler angles and integrals of body axis angular rates; for heave control 
the onl;; feedback is inertial axis heave rate (z ). Inputs from cockpit control- 
lers can be shaped by combinations of gains and integrators to provide rate 
command, attitude command, or rate command-attitude hold type systems. Forward 
loop control can be proportional or proportional plus integral. The flight 
control system includes a switch programmed as a function of airspeed which 
can be used to change control laws or cockpit controller function as the air- 
craft enters the hover regime. System gains can be programmed as functions 
of airspeed or pitch angle. (Airspeed programming is indicated on the control system 
diagrams.) As shown in Figures 2-2 and 2-4, the control system switch reverses 
the conventional roles of pedals and lateral stick deflection for hover. In 
hover the pedals command body axis roll while the lateral stick commands body 
axis yaw. This role switchirg is required in a VATOL in which t.he cockpit is 
rotated 90 degrees to be parallel to the ground in lover. This is the first 
concept to be studied in the NASA-Ames piloted VAXOL simulations. 
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Fig ire 2-1. RCS-Propulsion System Interactions Model 
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Figure 2-5. Generic Heave Control System Model 



The flight control syst at/ actuation systa interface provided in STSEQS 
is depicted in Figure 2-6. The prograasied aircraft can be controlled via pitch 
and yaw thrust deflection^ thrust modulation^ eleirons which combine elevator 
and aileron functions^ rudder^ leading edge wing flaps» trailing edge horizontal 
stabilizer (canard) flaps^ and RCS. In addition to providing roll control as 
indicated in Figure the RCS will accept pitchy yaw, side force, and normal 

force coBunds. 

With the exception of thrust coosiands, all the actuation system comnands 
defined in Figure 2-6 are processed by the generic rate and position limited 
actuator model shown in Figure 2-7. This actuator 0 K>del is programed as 
subroutine PCACT. The thrust commands are processed by the thrust dynamics 
model depicted in Figure 2-8. 

Subroutine SYSEQS has two sections of equations. One section, entered 
by a call to SYSEQS^ is used only for trim and derivative calculations. The 
other section, entered by a call to DEQU, is used only during time history 
calculations and contains most of the differential equations in the simulation. 
As may be Inferred cany of the control system statics and control system/ 
actuation interface equations are common to both sections while aircraft 
kinematics and control system and actuation system dynamics equations appear 
only in the DEQU section. 

2.2.3 DERIVS 


Subroutine DERIVS generates aircraft body axis stability derivatives 
about trim. It is called automatically during the trim program option. 

Changes in body axis forces and moments (X, Y, Z, L, M, N) are calculated 
for unit changes in the following variables: body axis airspeed components 

(u^^., w ); body axis rotation rates (p, q, r); elevator; aileron; 

ru^er; wing leading edge flap; horizontal stabilizing surface (canard) 
trailing edge flap; RCS roll, pitch, yaw, normal force, and side force 
conmands; thrust of each engine; pitch thrust deflection; and yaw thrust 
deflection. Since these variables encompass a wide range of control concepts, 
subroutine DERIVS should not rec^ire modification during any application of 
VATLAS. Total derivatives and the contributions to these derivatives from the 
wing, horizontal and vertical stabilizers, thrust, inlet ram, and RCS are 
calculated and printed by DERIVS. 

2.2.4 ERROR 


Subroutine ERROR contains the fifteen equations which must be balanced 
to attain a six DOF trim solution. These equations include balances between 
applied forces and linear accelerations, applied moments and angular acceler- 
ations, Inertial and body axis translation accelerations and rates, and Euler 
angle rates and body axis rotation rates. ERROR also contains user-selected 
logic for tiisming with or without ambient winds and in turns. ERROR is called 
only during trinming. To generate the data required for the fifteen trim 
balances ERROR cycles through FORCES and SYSEQS. ERROR is a fairly general 
routine and should not have to be modified for any application of VATLAS. 
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Figure 2-6. Interface of Flight Control System and Actuation System Model 
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Figure 2-*7, Generic Actuator Model 
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Figure 2-8. Propulsion System Dynamics and Performance Model 






2.2.5 UFDRG 


Subroutine LIFDRG contains the generalized equations for calculating the 
lift* drag, moment, and center of pressure shift of the aero lifting surfaces. 
LIFDRG req* ^res the surface angle of attack, leading and trailing edge flap 
deflections, and a series of aero constants to perform these calculations. 

The module should not have to be modified for any application of VATIAS. 

2.2.6 uxm 


Subroutine VTRUN controls the application of the aerodynamic model 
equations to generate wind timnel-type data. Inputs for WTRUK specify the 
type of sweep to be run, angle of attack Cs(} or sideslip i initial 
and final angles in the sweep; and the Incremental angle between data points. 
At each data point, total and component Cdue to wing, fuselage, horizontal 
stabilizer, and vertical stabilizer) nondimensional stability axis aero 
coefficients are calculated and printed out. Initial values may be placed 
on any variables which appear in the aero model so that dynamic derivatives 
and aero characteristics as functions of ci and ^ can be generated. For 
example, Cmq as a function of o( might be generated as follows: Run 

sweeps with pitch rate (q) * 0 and 0.05 or 0.1 rad/ sec then calculate 




ft»o 
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,0.05 or cj 
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where . 

■ A pitch moment coefficient ato( with q 

^o ° 

= mean aerodynamic chord of the wing 


A 


V. = reference airspeed (for convenience V. = 100 ft/sec 
can be used) 


WTRUN should not have to be modified for any application of VATLAS unless 
the user wants other than «< or p sweeps. 

2.2.7 NAMES 


In UPDATE terminology, NAMES is a COMDECK (common deck). During a 
FORTRAN compile of the VATLAS UPDATE source file, NAMES is inserted in 
subroutines DERIVS, FORCES, SYSEQS, ERROR, and VTTRUN. NAMES is simply a 
series of COMMON statements which order VATLAS variables and constants in the 
blank COMMON reserved by VTOLTH., the main program. As indicated in Appendix 
E, NAMES is heavily commented to clearly identify COMMON segments. The 
Importance of COMMON segmenting to VATLAS applications will become evident 
in subsequent sections. Variables and constants listed in NAMES, which 
constitute over 95% of those used in VATLAS, are available for program 
input and output under the control of VTOLTH. 
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2.3 Standard Program Modules 

The standardized program nodules Include TRIM, RUNGE, BATAM2, SIMQ, 

ATMOS, PCACT, and the table lookup and interpolation package. These modules 
have withstood the severe testing environment of extensive, critical use 
at Vought and should not be blamed for VATLAS failures to produce results; 
the message Is "DON'T MODIFY 1H£ STANDARDIZED MODULES." Since several of the 
modules are available to the user In developing program modifications, the 
formal parameters in their calling statements are defined below. 

The VATLAS plotting package Is a standardized module provided the user's 
computing facility supports the CALCOHP plotting utilities. The NASA-Ames 
CDC 7600 computing facility does not support CALCOMP, thus modifications to 
the Vought -supplied plotting module are required. Since the Vought routines 
have simplified the Input data required to generate plots. It Is recommended 
chat Che Vought Input format be retained during these modifications. A 
description of this format Is given In Section 4.0 when VATLAS Program Options 
5 and 10 are discussed. An exaaiple of a plot deck and the resulting plot Is 
Included In Section 5.0. 


2.3.1 TRIM 

Subroutine TRIM controls the search for trim conditions. It has a six 
degree-of-freedom trim algorithm based on Newton's method of solving systems 
of non-linear equations. Instructions for trimming the aircraft in various 
standard trim conditions (e.g. Ig wings level, cunstcuit heading sideslip, 
coordinated turn, etc.) are provided in the discussion of Program Option 6 
In Section 4.0. TRIM Is not a user specified module; It Is called automati- 
cally by Program Option 8. 

2.3.2 RUNGE 


Subroutine RUNGE controls and applies a 4th order Runge-Kutta algorithm 
for integration of the VATLAS differential equations during time history 
runs. Entry point SETUP In RUNGE Is called once at the beginning of a time 
history to initialize all the integrators. RUNGE Is not a user specified 
module; It Is called automatically by Program Option 4. 

2.3.3 BATAN2 


BATAN2 Is a function subroutine which determines the inverse tangent 
of Y/X in the range -U < tan~^- CT/X) ^ Tf .It avoids the 

problems of Infinite or Indeterminate operands when X*0orX*Y“0. 
BATAN2 Is a user specified subroutine whose application Is demonstrated In 
the following exanqple: 

ANGLE » BATAN2 (Y,X) 

where ANGLE ■ tan ^(Y/X) in rad 
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2.3.4 SIHQ 

Subroutine SIHQ solves a system of linear equations using a Gauss 
elimination algorithm. It is not a user specified module; It Is called 
automatically during trim calculations. 

2.3.5 ATMOS 


Subroutine ATMOS provides air density, ambient temperature and pressure, 
temperature and pressure ratios, and speed of sound as functions of pressure 
altitude for standard and tropical day conditions. The atmosphere model Is 
the U.S. Standard Atmosphere. ATMOS Is user specified with the following 
statement: 

CALL ATMCS (ALT, KATMbS, RH6, TAMB, PRESS, TRATI(J, PRATI0, 

VSdUND) 

where 

ALT ■■ pressure altitude, ft 

KATMbS » 1 fur standard day, 2 for tropical day 

RH6 “ air density, slugs/ ft 

TAMB - ambient temperature, deg R 

2 

PRESS = ambient pressure, Ib/ft 

TRATIQI - temperature ratio 

PRATI0 pressure ratio 

VSI^UND a speed of sound, ft/sec 
ALT and KATMOls must be defined before ATMt^S Is called. 

2.3.6 PCACT 


Subroutine PCACT Implements the generic rate and position limited 
actuator depicted In figure 2-7. It Is a user specified module and Is 
called as follows: 

CALL PCACT (DIN, D(AJT, RLP, RLN, PLP, PLN, TAU, ITRIM, DfiiUTD) 


where 


DIN - actuator Input 
D0UT actuator output 
RLP * upper rate limit 
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RLN ■ lover rate limit 


PLP ■ upper position limit 

PIN ■ lower position limit 

TAU actuator time constant, sec 

IX^D > rate o£ change of actuator output 

ITRIM indicates whether program is trimming or calculating 
a rime history and is automatically controlled by ‘ ■ 

program. 

DIN, RLP, RLN, PLP, PLN, and TAU must be defined before PCACT is cal 

2.3.7 Table Lookup and Interpolation Package 

The table lookup and interpolation module contains seven subroutines; 
STLU, nnu, INTRP, 1NTRP2, INTRP3, INTRP4, SEARCH. Only STLU and DTLU are 
available for user specification; the other subroutines are called as 
required during applications of STLU or DTLU. STLU performs single table 
lookup and interpolation to determine y = f(x) while DTLU performs double 
table lookup and Interpolation to determine z = f(x,y). The call to STLU is 
as follows: 

CALL STLU (X, FX, L, XARG, U, KdDE) 

where 

X is the table of independent variable values. 

FX is the table of dependent variable values. 

L is a fixed point constant equal to the number of entries 

in the X and FX tables. L is the minimum dimension of the 
X and FX arrays. 

XARG is the current value of the Independent variable. 

U is the current value of the dependent variable. 

K(i^ral specifies the type of interpolation to be used; 2 for 
two point (linear) interpolation, 3 for three point 
(quadratic) Interpolation, 4 for four point (cubic) inter- 
polation. KQIde ® 2 for calls to STLU currently in 
VATLAS. 

The X array should be loaded in order from smallest co largest value. Further- 
more, to avoid indeterminate or infinite operands, two successive X array 
entries should not be equal. 


20 



The call to DTLU is as follows: 


CALL DTLU (X, Y. FXY, L, M, XARG, U, KX, KY) 

where 

X la the table of the Independent variable 1 values. 

Y is the table of the independent variable Z values 

FXY is the table of the dependent variable values 

Lisa fixed point constant equal to the number of entries 
in the X table. L is the minimum dimension of the X 
array. 

M is a fixed point constant equal to the number of entries 
in the Y table. M is the minimum dimension of the Y 
array* FXY must have a minimum of L times M entries. 

XARG is the current value of independent variable 1. 

YARG is the current value of independent variable 2. 

U is the current value of the dependent variable* 

KX and KY take on the same values as K0DE in STLU and 

describe the type of interpolation to be applied along 
the independent variable 1 and 2 axes» respectively* 

As for STLU 9 the X and Y arrays must be loaded in order from smallest to 
largest value. Similarly, two successive X or Y array entries should not 
be equal. The FXY table is loaded in the order: f(X. ,Y.) to f(X^,Y.), 
f(X^,Y^) to fCX^.Y^), etc to f(Xj^,Y^) to f(X^,Y^)* 
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3.0 


VATLAS CCmO» ARRAY 


VTOLTH reserves 12000 words of blank COMMON in a single array called 

R. The purpose of the NAMES module is to label and arrange the R array with 

variables and constants used in subroutines DERIVS, FORCES, STSEQS, ERROR, 
and WTRUN. These parameters are arranged in logical groupings called segments 
where, for example, one segment ccmtains all the inlet ram related variables 
while another contains all the control system variables. A reference to the 
NAMES portion of the subroutine FORCES listing in ^pendlx B Indicates that 
the VATLAS CCWMON currently has forty five segments %ffalch consume approx- 
imately 3600 words of the R array. 

When devel ping a new simulation within the OLSIM framework, the user 

is free to astablish as many segments as he desires so long as the number ''f 

COMMON parameters does not exceed t-.e 12000 allocated in the R array. With 
a few exceptions, C0MM(RI segments can be arranged and niud>ered at trill; the 
exceptions are as follows: 

1. Segment 1 always has the sane seven parameters which are in 

order: T, TO, DT, TMAX, DTPRNT, DTPI^, TI. These COMMON 

parameters are defined in Appendix A along with all other 
parameters of the VATLAS CCRMON. 

2. Segment 2 contains all the variables to be integrated. At 
present a maximum of seventy are allowed. VATLAS uses 
sixty nine of this allocation. 

3. Segment 3 contains the Integrals of the Segment 2 variables 
aligned in the same order. Thus there are an equal number of 
Segment 2 and 3 variables and the first variable in Segment 3 
is the integral of the first variable in Segment 2, etc. 

Segment 5 contains the ranges or tolerances to which each 
component of the error vector in subroutine ERROR must be driven 
to achieve trim. These are constants which must be specified 
in the VATLAS Basic Load. VATLAS allocates nineteen locations 
for trim tolerances even though only fifteen are required in 
the trim solution. 

For COMMON input and output as occurs in the Basic Load and many 
of ti>e Program Options, the C0^f^0N parameters are designated by segment 
number and location within that segment. For example, T is element 1 in seg- 
ment 1, TMAX is element 4 in segment 1. 
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4.0 VATLAS BASIC LOAD AND PROGRAM OPTIONS 


Thla Mctlon provldaa loading inatructiona and deacriptlona of the Baalc 
Load for VATLAS and each of Che prograa opClona. PacaMtera auch aa NTAP, NTAPD, 
NTAPW , etc. which refer to Input or output file nuaibera required by VATLAS appear 
many tlaea In thla report aection. The only file nunbera currently available 
for uae when quantifying theae paraaetera are 1, 2, 3, 4, and 8 through 17. 

File nuabera 5, 6. and 7 have preapeclfled purpoaea: file S la the INPUT file 

which contalna input data on carda, file 6 la the OUTPUT file which contalna 
data for the output printer, and file 7 la the PUNCH file which contalna data 
to be punched on carda. To aake other file nuabera available to VATLAS, the 
flrat atataaant In VTOLTH (Appendix B) la aodlfled to Include TAPE XX, where XX 
la the new file nuaber. In Ita liat of foraal paraaetera. 

4.1 The Basic Data Load 

The Basic Data Load (BDL) specifies how the R or COMMON array is to 
be segaanted, a list of desired output paraaeters for use with Options 4 
or 12, and umy of chore program constants designated as type N parameters 
in COMION (see Appendix A) . Data for Option 9 %rhlch specifies the variables 
to be used for stability derivative detecalnatlon could also be considered 
in the BDL since these variables are not normally changed in the course of a 
VATLAS application. The list of output parameters will vary according to the 
type of run; two separate lists have been developed for VATLAS applications, 
one for summarizing data from a series of trim cases (Example I in Section 
3.0) and one for specifying Che variables to be stored during a time history 
run (Example 2 in Section 5.0). Only non-zero type N COMMON parameters require 
specification since standard run procedure is to initialize conputer storage 
to zero. 

The first two cards in the BDL are as follows: 


Card 

Format 




1 

115 

1-5 

6-80 1 



NEQS 

Blank 

2 

215 

1-5 

0 

1 

11 - 80 



NTAP 

Intapw 

Blank 


where 


NEQS is the number (^70) of differential equations in Che problem 
- 69 for VATLAS 

NTAP la the number of the file vdiich contains the BDL. NTAP *• 3 
if the BDL is card input. If NTAP J* 5, the appropriate file 
must be ATTACHed in the ('onlrol Card stream before entering 
VATLAS. 
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NXAPW Is the nuaber of the file where the BDL Is to be copied. This 
allows the BDL to be stored oe e file and thus elieinete a 
significant nuaber of cards in the data deck. If MIAPW “ 0 
the BDL is not stored. If NTAPW ^ 0, file NTAPW must be 
AXTACHed before entering VATLAS, then CAXALOGed before ending 
the Control Card stream to insure that the data will be saved. 

These two cards will always be required even though the BDL is stored on 
a file. 

The next few cards in Che BDL specify the nuaber of segments in the R 


array and 
Card 

the nuaber of storage 
Format 

locations 

in each segment: 


3 

1615 

1-5 

6-10 

11 - 15 

« • • 

76 - 80 



NSEC 

LENGTH(l) 

LENGTH(2) 

« • * 

LENGTH(15) 

4 

1615 

1-5 

6-10 

11 - 15 

• • • 

76 - 80 

• 


LENGTH(16) 

• • « 

» » « 
• 

• • * 

LENGTH(31) 

2+NS 

1615 

1-5 

6-10 

11-15 

* « • 

76-80 



• • » 


! 

• • * 

LENGTH(NSEG) 



where 


NSEC is the nunbet (& 200) of segments in the R array. 

- 45 for VATLAS 

LENGTH(l) Is the number of storage locations In segment 1. 

rounded to next highest integer value. NS ■= 

- 2.875 “ 3 tor VATLAS 

The next subset of BDL cards specify the list of desired output 
parameters for Option 4 (Time history) or 12 (Trim data summary): 


Card Format 


NS+1 215.E10.3, 

1-5 

6-10 

11 - 20 

21 - 25 

26 - 30 

31 - 78 

2A5,12A4 

I 

J 

SCAL(l) 

LBL(l) 

LBL(2) 

IJdiEL 


NS 


NSEGt-1 

16 

45-t-l 

16 
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Card Format 

NS+NtfUT 2I5,E10.3, 

2A5,12A4 


NS4N0UT+1 


where 

N#UT is the number (^120) output parameters 

I is the segment number of the Kth output parameter 

J is the element number of the Kth output parameter 

SCAL(K) is the constant %ihich scales the Kth parameter before it is 
output. An example of the use of SCAL(K) follows: most 

angles in the computer are in radian measure; it is desirable 
trac these angles be output in degrees, therefore SCAL(K) 
for an angle parameter would be 57.296. In most cases, 

3CAL(K) will equal 1. since the computer variable has its 
. normal units. 

LBL(2*K“1) and LBL(2^K) contain the first five and second five 

alphanumerics, respectively, in the output parameter name. 

This name is not constrained to be the same as the coded 
name of the parameter given in Appendix A. 

LABEL is an optional field and is used to further identify an 

output parameter on the computer printout. LAP ’.L is not 
stored in the computer but is simply passed from input to 
output. 

The list of output parameters is ended by a blank card. Even if no output 
parameters are desired (i.e., N0UT =0), the blank card is required for the 
program to continue. 

The final and usually longest rubset of BDL cards directs the loading 
of the program constants. As noted above these are generally the type N 
parameters in COMMON. This subset and, equivalently, the BDL is ended by 
a blank card. Hie form for loading program constants is as follows: 
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Card 


Format 


NS4M(IUT-l-2 

4IS.12A4. 

1-5 

6-10 

11 - 15 


- 20 

21 - 68 

67 - 70 

71 - 80 

(An input/ 
output spec- 
ification 
card) 

2X.F10.6 

I 

J 

K 

L 

LABEL 

Blaak 

SCALE 

ms+n6ut+ 3 
(A para- 
nibrer value 
card) 

7E10.3 

1-10 

11 - 20 

« 

. . 

71-80 



Value of 
a»M)N 
parameter 
in Segment 
1 Element J 

Value of 

parameter 
in Segment 
I Element 
J+1 

* 

• 

• 



Blaak 


NS4M^UT+2 

7E10.3 

1-10 

11 - 20 

• > 

• • 

71 - 80 


-HJCARD 
(A para- 
TDeter value 
card) 


* « • 

• • • 

Value of CCXIMCXI 
parameter in 
Segment K 
Element L 

Blank 



NS+n6uT+3 

4I5,12A4 

1-5 

6 - 1C 

11 - 15 

16 - 20 

21 - 68 

69 - 70 

71 - 80 

-HICARD 

2X.F10.6 

I 

J 

K 


LABEL 

Blank 

SCALE 


Last Card 
of BDL 


where 


1-80 


Blank Card 


1, J, K, L specifies Chat the data contained on the next NCARD cards 

will fill all the locations in C0^t10N from Segment 1, Element J 
to Segment K, Element L. Thus (ILOC pieces if data must be 
contained on Che NCARD cards. 


NL0C is the number of COMMON locations from Segment I, Element J to 
Segment K, Element L. 

NCARD = rounded to next larger whole number, (i.e., to load 40 

(=NL0C) locations requires 6(=NCARD) cards. 
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LA3EL allows the user to define what data are being loaded. This 
trill appear on the co^iuter printout and will help identify 
what the data are. LABEL is sinply passed fron input to output 
and is lost tihen the next LABEL specification is loaded. 

SCALE specifies the constant which scales the data on the next NCARD 

cards before they are leaded into COtMON. For example, aircraft 
fuselage stations, buttlines, and waterlines are generally 
given in inches; the computer finds it more convenient to work 
in feet. Thus SCAIX would be 0.0833 and the card data tfould 
be in inches for aircraft locations. 

Note that if the user does not care to identify his data input other 
than it is CCBtlQN data, he can load all of COtMON with one specification 
card where, for VATLAS, I*J*1, K«45, L >170. This would be folloired 
by 513 . 3585 locations (NL^) in VATLAS COfCH _ ^ data cards 

and a final blank, card to end the BDL. 

A copy of BDL cards for the VATLAS application to the Vought SF-121 
airplane is given in Appendix C. The playback of this Information during a 
computer run is shown in the printouts for Examples 1, 2 and 3 contained in 
Appendices D, E, and F, respectively. 


4.2 Program Options 

VATLAS has 12 options for managing program execution, input, and output. 
These are li. ted and described briefly below; their use and loading are 
detailed in subsequent paragraphs: 

Option Function 

1 Create a file of input or output data specifications 

2 Read data into COMMON array 

3 Write data from COMMON array 

4 Ca..culate a time history 

5 Plot a time history 

6 Set up trim unknowns 

7 Load variable formats for trim sunnary or wind tunnel runs 

8 Calculate trim 

9 Specify variables for stability derivative calculations 

10 Plot a time history stored on a binary file during a previous run 

11 Calculate wind tunnel data 

12 Print trim data summary 
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The data cards for each option fora Independent input nodules. These 
modules are placed after the BDL in the order desired for program execution. 
Since some options depend on information supplied by other options this module 
ordering must adhere to certain guidelines. For example, options 4, 5, 8, 11 
and 12 cannot be used alone; their application is governed by the following 
rules ; 

1. Options 6 and 9 must be placed before Option 8 which means 

that the aircraft cannot be trlnmed until the trim and stability 
derivative variables have been specified. 

2. Option 8 oust be placed before Option 4; the aircraft must be 
trlnmed before a time history can be calculated. 

3. Option 4 must be placed before Option 5; the tine history 
must be available before it can be plotted. 

4. Options 7 and 8 must be placed before Option 12; the aircraft 
must be trlnmed and the formats for printing a data summar^ 
must be available before a trim summary can be printed. 

5. Option 7 must be placed before Option 11; the formats for printing 
wind tunnel data must be available before the data can be printed. 

Options 1, 2 and 3 can be used alone and can be placed almost anywhere in the 
sequence of options. The use of Option 10 requires only that a previously 
calculated and stored time history is available. It does not require a 
full BDL and frequently constitutes an entire computer run (i.e. a plot run). 

4.2.1 Option 1 - Store Lists of Input or Output Parameters 

Option 1 causes a list of input or output parameters to be stored on 
a local file. This Option is convenient when a series of trims or time 
histories are produced on a single computer run and the same set of parameters 
are to be input or output for each trim or time history. The format for 
Option 1 data cards except the first two and the last is that of an input/ 
output specification card which was described in the section: 


Card 

Format 



1 

115 

1-5 

6-80 



N^PT 

Blank 

2 

115 

1-5 

6-80 


1 

NTAP 

Blank 
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Card 


Fonnat 


All others 

415* 12A4, 

1-5 

6-10 

11- 15 

16 - 20 

21 - 68 

69 - 70 

71 - 80 

except 
last card 
(Input/ 

2X.F10.6 

I 

J 

K 

L 

LABEL 

Bl.-^nk 

SCAI£ 


output 

specifi- 

cation 

cards) 


Last card 


where 


1-80 


Blank. Card 


M0PT is the option number which is equal to 1 in this case. 

NTAF is the number of the local file where the input or output list 
is stored. 


Ihe parameters I, J, K, L LABEL* and SCALE were described in Section 
4.1. When Option 1 is applied* it must be placed before the trim (Option 8) 
or time history (Option 4) options. 

4.2.2 Option 2 - Read DAta Into COMMON A-ray 

Option 2 allows data to be read into the CCMdMON array under the direction 
of either a stored (by Option 1) list of input specifications or input 
specification cards in the input card stream. If a stored list is used then 
only the parameter value cards need to occur in the input data stream. If 
a stored list is not used then the parameter value cards are separated by 
input specification cards in exactly the same way that constants are loaded 
in the BDL. Option 2 is ordinarily used to initialize aircraft state and 
control variables; set up pseudo-pilot* control system, or aircraft con- 
figuration options; or change aircraft or control system constants before 
entering a trim, wind tunnel* or time history run. The sequence for Option 
2 cards is as follows: 


Card 

Format 



1 

115 

1-5 

6-80 



N0PT 

Blank 


2 

115 

1-5 

6-80 



NTAP 

Blank 
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Card 


Format 


3(NTAP-5) 

4I5,12A4, 

1-5 

6-10 

11 - 15 

16 - 20 

21 - 68 

69 - 70 

71 - 80 

(An input 
specifi- 

2X.F10.6 

I 

J 

K 

L 

LABEL 

Blank 

SCALE 


cation 

card) 


3(NTAPi‘5) 7E10.3 

1-10 

11 - 20 

• • • 

71 - 80 


Value of 

Value of 

• • • 

Blank 


first 
COMMON 
parameter 
specified 
by first 
input 
specifi- 
cation 
in file 
NTAP 

second . . . 




If NTAP > 5 , alternately load Input specification 
cards and parameter value cards as required until 
last card. 

If NTAP ^ 5, load parameter value cards until last 
card; insure that the parameter value for the first 
parameter In each Input specification on NTAP begins 
a new card. 


Last card 


1-80 


Blank, card 


where 


N0FT 2 indicating Option 2 

NTAP is the number of the local file where the input specifications are 
found. NTAP - 5 is the INPUT file. NTAP ^ 5 specifies the 
file where Option 1 stored the input specification list. The 
parameters I, J, K, L, LABEL and SCALE vere described in Section 
4.1. 


4.2.3 Option 3 - Write Data from COMMON Array 

Option 3 allows data to be written to OUTPUT from the COMMON array 
under the direction of either a stored (by Option 1) list of output 
specifications or output specification c. .ds in the input card stream. If 
a stored list is used then only the first two cards in the Option 3 load 
are required. If a stored list is not used, then the Option 3 load must 
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Include the appropriate output specification carda. Direct use of Option 3 
through the M<(PT route la rare in noraal applications of VATIAS; Instead 
automatic transfer to Option 3 capabilities Is provided In Options 4 and 8 
for printing out Initial conditions for a time history or trim data at the 
end of an individual trim case. Hie sequence for Option 3 cards Is as 
follows : 


Card 

Fonsat 



1 

115 

1-5 

6-80 



N0PT 

Blank 


2 

1L5 

1-5 

6-80 



NTAP 

Blank 


If NTAP ^5. the output specifications are found on the NTAP 
file and no more cards are required. 


3 4I5.12A4 

(NTAP-5) 2X,F10.6 
(An in- 
put/out- 
put spec- 
ification 
card) 


1-5 

6-10 

11 - 15 

16 - 20 

21 - 68 

69 - 70 

71 - 80 

I 

J 

K 

L 

LABEL 

Blank 

SCALE 


Last Card 
(NTAP-5) 


1-80 


Blank Card 


where 


N0PT - 3 indicating Option 3 

NTAP is the number of the local file where the output specifications 
are found. NTAP - 5 is the INPUT file. NTAP 0 5 specifies the 
file where Option 1 stored the output specification list. 

The parameters 1. J» K, L. LABEL, and SCALE were described In Section 

4.1 except in this instance SCALE is the constant which scales 

the COMMON array data specified by I, J, K, L before they are printed. 

4.2.4 Option 4 - Time History Calculation 

Option 4 controls the calculation of one time history. Multiple time 
histories during the same computer run require multiple sets of Option 4 
input cards. The sequence for one set of Option 4 input cards is as 
follows : 
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Card 

Format 



1 

115 

1-5 

6 - 80 



N0PT 

Blank 


1-5 

6-10 

11 - 15 

16 - 20 

21 - 25 

26 - 80 

HTAPE(l) 

MIAPE(2) 

NXAPE(3) 

NTAFE(4) 

MTAFE(5) 

Blank 


If 

3 

4 

5 

6 
7 

If 

3 


(NTAPE(4) ■ 5)t cards 3 to 7 are 



where 


N0PT 4 Indicating Option 4 

NXA.PE(1) specifies where Input specifications are found when Option 4 
automatically transfers to Option 2. 

NTAPE(2) specifies where output specifications are found when Option 4 
automatically transfers to Option 3. 

NTAPE(3) specifies the local file where time history calculations 
are to be stored for plotting 
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NTAPE(4) specifies vhere the variable formats (FKIl, FHT2, etc.) 
are found. If NTAPE(4) ■ S, the formats are found in the 
input card stream. If HTAPE(4) i S« the formats are found on 
local file NTAPE(4) . 

NTAPE(5) specifies where the variable formats are to be stored for 

future use during multiple time history runs. This parameter 
eliminates the requirements for cards 3 through 6 in. all but the 
first of several time histories. NTAPE(5) would be nonzero 
only for the first time history. For subsequent time histories, 
HIATC(S) would be zero and NTAPE(4) would be set to the nonzero 
NXAPE(S) value used in the first time history. 

FKIl is the format used to label the columns of the time history output. 
This format was the same for all VATLAS applications to date: 

C7X, 11C2A5)) 

FKT2 is the format of the data in the colums of the first page of the 
time history output. This format was also the same for all 
VATLAS applications to date: (5X, 11F10.3) 

FMT3 is the format of the data in the columns of the second and 
subsequent pages of the time history output. This format 
was also the sane for all VATLAS applications to date and 
is the same as FMT2: (5X, 11F10.3) 

FM14 ia not currently used. 

TITLE is the title which will appear at appropriate places in the 
computer printout. 

An automatic transfer to Option 2 occurs after cards 3 or 7 (depending 
on NTAPE(4)) of the Option 4 load. In this transfer NTAP in Option 2 Is 
set equal to NTAPE(l). This transfer provides for COMMON array changes 
before the Integrators are initialized. The program now expects to see 
Input data equivalent to that beginning at card 3 in the Option 2 load and 
continuing to the blank last card. If the user does not want to make R 
array changes he must load a blank card as card 4 or 8 (depending on NTAP£(4)) 
in the Option 4 card sequence. 

The transfer to Option 2 is followed by an automatic call to SETUP to 
initialize the integrators an' then an automatic transfer to Option 3 where 
NTAP in Option 3 is set equal to NTAPE(2) . This transfer provides for 
printing out any COMMON array parameters before the time history calculations 
are initiated. The program expects to see input data equivalent to that 
beginning at card 3 in the Option 3 load and continuing to the blank last card. 
These data will be found either in the Option 4 input card stream (NTAPE(2) 

5) or in a local file (NTAPE(2) 5) . If the user has not stored an appro- 
priate output specification list and does not want to printout any COMMON 
airay parameters he must load a blank card in the Option 4 card sequence. 
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Hote that the tloie history output variables which are priuted out and 
stored oa local file MTAPEC3) (i£ HTAFE(3) f 0) are those specified in the 
output parameter specification section of the BDL. 

4.2.5 Option 5 - Time History Plotting 

Option S controls the plotting of a single series of px./ts from a 
single time history run; it must be Invoked for each series of plots to be 
generated. The stibroutines which support Option S have standardised and 
simplified the Interface with CALCOlff routines. QfkTLAS will have to be 
modified if the user's computer facility does not support CALCOMP utilities.) 
Option 5 can be invoked in a computer run after a local tine history data 
file has been created by an application of Option 4. When the time history 
data file is stored on one computer run and then plotted during a subsequent 
computer run. the Option 5 support subroMtlnes are entered from Option 10. 

A sample data set and series of plots are provided in Example 4 in Section 5.0 
to illustrate the generalized Option 5 Inputs described beloc: 

Card Format 

1 115 


Z 8A10 


3 2F10.3 

4110 



3+1 


313, IX, 

2A10, 

2F10.2, 

4F5.1, 

215 
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Card Forma*: 


3+NVAR 313, IX, 
2A10, 
2F10.2, 
4F5.1, 
215 


where 


1-3 

4-6 

7-9 

10 

11 - 20 

21 - 30 

N0(NVAR) 

ISYM(NVAR) 

IPT(NVAR) 

Blank 

HEA0(1 ,NVAR) 

H£AD(2,NVAR) 

31 - 40 

41 - 50 

51 - 55 

56 - 60 

61 - 65 

66 - 70 


SV(NVAR) 

SDV(NVAR) 

b(NVAR) 

C(HVAR) 

YA0(NVAR) 

X2SHFT(NVAR) 



71 - 75 

76 - 80 

ITIME(NVAR) 

NTITLE(NVAR) 


N0FT •• 5 Indicating Option 5 


TlTLEl Is the first line of the plot title 


T1TLE2 .is the second line of the plot title 

XI directs plotter pen to move from Its reference to the left hand side 
of the first plot In the series. (A single "plot" consists of a 
suitably labeled and scaled ordinate axis and one or V o time 
history traces. The time (or abscissa) axis may or may not be 
present.) Positive XI (In Inches) moves pen to the right. 

YI directs plotte** pen to move from Its reference tc the top of the 

first plot In the series. Positive Yl (In Inches) movej pen down 
the page. The reference point for XI and YI (X'C ° Y.r ~ 0) Is 
generally at the upper left comer of the grid. 

NVAR is the number of variables (^0.. to be plotted. Time must be 

Included ^.s a variable In this count; for example a plot of elevator 
deflection versus time requires KVAR ■ 2. 

N0READ ^ 0 means that the next NVAR cards In the Option 5 card sequence 

should be read; l.e. read cards 3+1 to 3+NVAR. N0READ = 0 means 
that the current series of plots Is to be plotted according to the 
same Instructions as the previous series and cards 3+1 to 34NVAR 
are eliminated from the Option 5 card sequence. 

I0P “ 1 for first series of plots, ■ 0 for all other plot series. 

ICL * 1 for last series of plots, • 0 for all other plot series. 

N0(I) Is the number of the variable to be plotted as determined from Its 
position In the output specification list of the BDL. The first 
variable in the N0(I) array must be time. 

ISYM(l) Is the number of the CALCOMP symbol to be placed every IPT(I) 
data point along the time history trace. The fourteen symbols 
available and their Identification numbers are given In Table 4.1. 

If a symbol Is not desired, then ISYM(I) -1. 
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IPT(I) Is the spacing in terms of data points between the syebols on 
the tine history trace. If ISIMCI) “ -1* IPTCl) Is Ignored. 

HEAD(1,I) and HEAD(2,1) define the labeling for the plot axis. 

SV(I) is the value of the plotted variable at the lower end of its 
axis. 

SDV(I) is the scale factor of the plotted variable In units/ inch 

If SV(I) and SDV(I) are left blank, their values are detrriained 
by the CAL COMP automatic scaling option. 

B(l) is the length of the axis to be drawn for variable N^(l) in inches. 

C(I) is the distance in inches down the page bet%reen the lower left hand 
comer of the axis of the plot of variable H0(I) to the upper 
left hand comer of Che axis of the plot for variable N^Cl+l). 

C(I) is positive down the page. 

YA^(I) is the distance in inches across the page between XI and the 
ordinate axis of variable N0(I). YA0 Is positive to the 
right. 

X2SUFT(I) is the distance in inches across the page between the first 
point of variable M^(I) (i.e. the value of N0(I) at initial time) 
acd Che first point of variable H^(I-l-l) . X2SHFT(I) is positive to 
Che right. 

rriME(l) specifies whether Che time axis is to be plotted on Che plot 
for variable N0(I) . If ITIME(I) =* 1, the time axis is plotted, 

0, Che tine axis is not plotted. 

NTITLE^I) specifies whether TITLEl and T1TLE2 are to be printed at 
ti4e botccat of Che plot for variable N0(I) . 


Table 4-1 - Symbols Availablefor VATLAS Tine History Plots 


ISYM(I) 

SYMBOL 

1SYM<I) 

SYMBOL 

0 

□ 

7 

X 

1 

o 

8 

Z 

? 

A 

9 

Y 

3 

+ 

10 

X 

U 

X 

11 

)K 

5 

c> 

12 

X 

6 


13 

1 
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1^0 varlAhles can be plotted on the sane azla systen* This le eutonatlcally 
assmed by the ptograai vfaen the vertlcel scales Coc two successive vaclables 
are placed at cuLSCtly the sane place CCI) * CCl+1) ®nd YA^(I) “ 

YA0(I-»-l)). The second variable, the (1*1) th, will be plotted as a dashed line. 

4.2.6 Option 6 - Set gp Trln Ifalmowna 

Option 6 Is applied to specify the fifteen Independent variables which 
are to be perturbed to find a trla solution for the fifteen error equations 
In subroutine ERROR. Several standard sets of trln variables are specified 
by trla type In Table 4-2. These sets span a wide range of potential VATLAS 
applications. They do not begin, however, to exhaust all the variable sets 
which could be developed to trla the aircraft. Applications of two of the 
standard sets are deaonstrated In Exanples 1 (Transition Trim Variables) and 
2 (Ttla In a Wind Variables) In Section S.O. The Input card sequence for 
Option 6 la as follows: 


Card Fomat 


1 

115 

1-5 

6-80 



NtfPT 

Blank 


2 

115 

1-5 

6-80 



MEQU 

Blaak 


3 4I5,12A4 

1-5 

6-10 

11 - 15 

16 - 20 

21 - 68 

69-80 


I 

J 

K 

L 

LABEL 

Blank 


4 7E10.3 

1-10 

11 - 20 

« » • 

71 - 80 


DZ(1) 

DZ(2) 

• » % 

Blank 

4+«CARD 7E10.3 

1-10 

* • • 

0 
00 

1 


DZ(8) 

DZ(NL0C) 

Blank 


4-HlCARD 

415, 

1-5 

6-10 

11 - 15 

16 - 20 

21 - 68 

69 - 80 

+1 

12A4 

1 

J 

K 

L 

LABEL 

Blank 


4+NCARD 7B10.3 

1-10 

11 - 20 

• • • 

71 - 80 

+2 

DZ(NL0C4-1) 

DZ(NL0C+2) 


Blank 


Last Card 


1-80 


Blank Card 
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where 


N0PT " 6 indicating Option 6 

MEQU ■■ 15 - nunber of error equations in subroutine ERR0R 

1, J» fit L specifies that all the COIMOH variables froa Segaent 1 
Eleaent J to Segaent K. Eleaent L are to be tria variables. 

LABEL allows the user to identify these variables for printout. LABEL 
is not stored in the computer and is lost when the next LABEL 
specification is loaded. 

DZ(1) specifies the perturbation tc be applied to tria variable 1 

when tria aatrix sensitvities to trim variables are calculated. 
The tria variables are nuabered 1 to IS in the order they are 
loaded; for exaaple, tria variable 1 is the variable which is 
specified by card 3 to be Eleaent J in Segment I. 

NL0C is the nunber of COMIQN locations from Segaent I, Element J to 
Segment K, Element K. 

NCARD " — J — rounded to next larger whole number. 

4.2.7 Option 7 - Load Variable Formats 

The foromts for printing the «rind tunnel option (Option 11) outputs and 
the trim sunaary (Option 12) are loaded by Option 7. These sets of formats 
have been standardized for VATLAS applications and are given in Examples 1 
(trim summary formats) and 3 (wind tunnel option formats) in Section 5.0. 

The input card sequence for Option 7 is as follows: 


Card 

Format 




1 

115 

1-5 

6-80 



N0PT 

Blaak. 


2 8A10 

1-10 

11 - 20 

» • • 

71-80 

• 

FKTS(1,1) 

FMTS(1,2) 

* • » 

FMTS(1,8) 

13 8A10 

1-10 

11 - 20 

• • • 

71 - 80 


FKTS(12,1) 

FHTS(12,2) 


FMTS(12,8) 
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where 


N0PT “ 7 indicating Option 7 

FMTS(I,1)» FMIS(It2) .IMIS(I,8) specifies Ith variable fonat 

There mist be 13 cards in the Option 7 load. If less than t%ielve 
fo teats are required, load blank cards to main tain the total equal 
to 13. 

4.2.8 Option 8 - Trie Calculation 

Option 8 controls the calculation of tries. The sequence for Option 8 
input cards is as follom: 

Card Foreat 


1 


2 


3 


115 

815 


20A4 


1-5 

6-80 

N#PT 

Blank 


1-5 

6-10 

11 - 15 

16 - 20 

21 - 25 

26-30 

31 - 35 


HTAPE(l) 

NTAPEC2) 

NTAPE(3) 

M1APE(4) 

NCASES 

NTAPD 

IPUNCN 




36 - 40 

41 - 80 


IWRITE 

Blank 


0 

GO 

1 


TITLE 


where 


N0PT 8 indicating Option 8 

NTAPECl) specifies where input specifications are found when Option 8 
automatically transfers to Option 2. 

MTAPE(2) specifies where output specifications are found when Option 8 
automatically transfers to Option 3. 

NTAPE(3) specifies where the data for the trim summary option (Option 
12) are to be stored. 

NTAFE(4) is not used. 

NCASES is the number of trims to be calculated. 


<0 



MTAPD apaclflad «hare the verlebles for stability derivative calculations 
are stored by Option 9. 

IPUMCH ^ 0 causes stability derivatives snd initial conditions to be 
punched out for offline linear analyses. IPUNCH « 0 suppresses 
punched output. The list of data output by this directive are 
detailed In Table 4-3. This list uas developed specifically for 
the SP-121 airplane application. The liat is contained in DERIVS 
and can be easily altered by UPDATE directives. 

IWRITE ^ 0 causes Intemediate trie calculations to be printed out as 
an aid to debugging the program. IWRITE “ 0 suppresses the 
intemediate trie calculations. 

TITLE is the title which will appear at appropriate places in the 
computer printout. 

An automatic transfer to Option 2 occurs after the program has read 
card 3 of the Option 8 load. In this transfer NTAP in Option 2 is set equal 
to HTAPE(l) . This transfer provides for OCRMON array changes prior to 
triaming. The program now expects to see input data equivalent to that 
beginning at card 3 in the Option 2 load and continuing to the blank last 
card. If the user does not want to make COMMON array changes he must load 
a blank card as card 4 in the Option 8 card sequence. 

Option 8 next transfers control to the trim routines. When a trim has 
been found an automatic transfer to Option 3 occurs. In this transfer 
NTAP in Option 3 is set equal to NTAPE(2) . This transfer provides for 
printing out the trim results. The program expects to see input data 
equivalent to that be ginning at card 3 in the Option 3 load and continuing 
to the blank last card. These Option 3 data will be found either in the Option 
8 input card stream (NTAPE(2) « 5) or in a local file (MTAFE(2) # S). If 
the user has not stored an appropriate output specification list and does 
not want to print out any COMCRi array parameters he must load a blank card 
in the Option 8 card sequence. 

Also at the end of each trim, the program automatically loads NTAPE(3) 

(if it is ft 0) with the values of those parameters specified in the list of 
output parameters in the BDL. If NIAPE(3) » 0. this procedure is bypassed 
since the user does not intend to print a trim data summary later in the computer 
run using Option 12. 

The automatic transfer to Option 2 occurs before each trim case specified 
by the MCASES parameter. Similarly the transfer to Option 3 and loading 
NTAPE(3) occurs after each trim has been determined. 

4.2.9 Option 9 - Set Up Stability Derivative Variables 

Option 9 specifies the variables to be perturbed to generate stability 
derivatives. Since subroutine I^RIVS, which controls stability derivative 
calculations, is fairly general and is based on a specific set of variables. 
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Table 4-3 List of Tunebad Output Outa frai Optioo 8 (Tria Calcolatioii) 



1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

61-70 

1 

VAS 

HAS 

Q 

THETA 

THETD 

ALPH 

VA 

2 

IX 

n 

IZ 

IXZ 

VT 

JERG 

— 

3 

x/u 

. x/w 

X/Q 

X/DB 

X/DN 

X/DFC 

X/TBTl 

4 

X/1BT2 

x/n 

X/T2 

— 

— 

— 

— 

5 

Z/V 

z/w 

Z/Q 

Z/OB 

Z/DFW 

Z/DFC 

Z/THTl 

6 

Z/TBT2 

Z/Tl 

Z/T2 

— 

— 

— 

— 

7 

H/U 

M/W 

H/Q 

H/DB 

M/DFW 

H/OFC 

M/1BT1 

8 

H/THT2 

H/n 

M/T2 

— 

— 

— 

— 

9 

T/V 

T/P 

T/R 

T/OA 

T/nt 

T/psm 

T/PS1T2 

10 

T/Tl 

T/T2 

T/DRCSR 

— 

— 

— 

— 

U 

L/V 

L/r 

L/R 

1/04 

L/DR 

Wpsm 

1/PSIT2 

12 

I7I1 

L/T2 

L/DRCSR 

— 

— 

— 

— 

13 

■/V 

H/P 

H/R 

R/04 

H/DR 

R/psm 

B/P8IT2 

14 

H/Tl 

H/T2 

R/inCSR 

— 

— 

— 

— 


vhara: 

UASfWAS ■ trim airspeed components along X and Z body axes (ft/sec) 

Q^IHETD “ trim pitch rate (rsd/sec) 

THETA ■ ^rim pitch angle (rad) 

ALFB - trim pitch angle of attack (rad) 

VA - trim airspeed (ft/sec) 2 

IX, IT^1Z,IXZ • aircraft moments of Inertia (slug ft ) 

VT “ aircraft weight (lb) 2 

JENG * engine rotational moment of inertia (slug ft ) 

I/J ■ change in 1 acceleration per unit change In J variable where: 

1 - X,Y,Z,L,M, or N and J - U.W.Q,DE,DW,DFC,T!m,THT2.Tl,T2*V,P,R,DA, DR, 
?SIT1,PSIT2. or DRCSR 2 

X, Y,Z " translational accelerations along the X, Y, and Z body axes (|t/sec ) 

t,M,N "angular accelerations about the X»Y, and Z body axes (rad/sec ) 

U,V,W ■ airspeed components along X,Y, and Z body axes (ft/ sec) 

P,Q,R • angular rates about X,Y, and Z body axes (rad/sec) 

DE ■ aerodynamic pitch control deflection (rad) 

DFV ■ wing leading edge flap deflection (rad) 

DFC • canard or horizontal stabilizer trailing edge flap deflection (rad) 
THT1,THT2 • pitch thrust deflection of left and right engines (rad^ 

T1,T2 ■ thrust of left and right engines (lb) 

DA ■ aerodynamic roll control deflection (rad) 

DR “ aerodynamic yaw control deflection (rad) 

PS1T1,PS1T2 " yaw thrust deflection of left and right engines (rad) 

DRCSR • roll RCS normalized deflection 
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there should be no reason to ever change that set for any V ATLAS applications. 
Option 9 could be conaldered. therefore as a part of the BOL. The VATLAS 
Option 9 data cards are given In the examples of Section 5.0 and are formatted 
as described below: 


Card 

Format 




1 

115 

1-5 

6-80 



N0PT 

Blank 

2 

115 

1-5 

6-80 



NTAPD 

Blank 


3 

315 

1-5 

6 - 10 

11-15 

16-80 



I 

J 

NRD = 1 

Blank 


1+NVAR 315 

3-HlVAR — 

+1 

where 


1-5 

6-10 

11 - 15 

§ 

» 

vO 

I 

J 

NRD=NVAR 

Blank | 


1-80 


Blank. Card 


N0PT > 9 Indicating Option 9 

NTAPD is the local file where the variable list is to be stored 


I is the segment number of the NROth variable 
J is the element number of the NRDth variable 

NRD is Che variable nunber within the list of variables. Normally 
Che variables are loaded in order with NRD beginning at 1. 

NVAR is the number of variables to be perturbed to generate stability 
derivatives. NVAR *• 22 for VATLAS. 


4.2.10 Option 10 - Plot a Previously Stored Time History 

Option 10 controls the plotting of variables from a time history 
generated and stored during a previous computer run. Its input differs for 
chat of Option 5 only In Its first two cards: 
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Card 


Format 


1 

115 

1-5 

6-80 




N0PT 

Blank 


2 

315 

1-5 

6-10 

U - 15 

16-80 



NPTS 

N0UT 

NTAP 

Blank 


where 

N0PT - 10 ladlcatlag Option 10 

NPTS is the nuober of time samples in the time history ( ^ 302) 

N0UT is the number of output variables recorded at each time 
sample ( ib 120) 

NTAP is the local file containing the time history data. 

Starting with card 3, Option 10 input is Identical to the Option 5 
input which begins at card 2. In fact Option 10 transfers to Option 5 
as soon as the above two cards are read. 

Note that if the only purpose of a computer run is to plot a previously 
stored time history using Option 10, the output parameter list and the loading 
of program constants portions of the BDL can be deleted and replaced by two 
blank cards following that portion of the BDL which defines COMMON segment-^. 

4.2.11 Option 11 - Wind Tunnel Calculations 

Option 11 controls the application of the VATIAS aerodynamic model to 
calculate nondimenslonal coefficients in stability axes. Before Option 11 
can be invoked, the variable formats required to print the wind tunnel 
calculations must be loaded by Option 7. The standard formats adopted for 
VATLAS application are given in Example 3 of Section 5.0. The input card 
sequence for Option 11 is as follows: 


Card 

Format 



1 

115 

1-5 

6-80 



N0PT 

Blank 


2 

315 

1-5 

6-10 

11 - 15 

16 - 80 



IPITCH 

lYAW 

NTAP2 

Blank 


3 

7E10.3 

1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

61-70 

71-80 



A0A0 

DA0A 

A0AMAX 

PSI0 

DPSl 

PSIHAX 

WIND 

Blaak 
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trher« 


N0PT •• 11 indicating Option 11 

IPITCU ^ 0 indicates that an angle of attack sweep is to be simulated. 
IPITCH " 0 indicates angle of attack is to remain constant for 
tte run. 

lYAW 0 indicates that a yaw sweep is to be simulated. lYAW =* 0 

indicates yaw (sideslip) angle is to remain constant for the run. 

NTAP2 is the local file for storing the calculations for later printing. 

is the initial angle of attack in degrees for an angle of attack 
s%>eep run. 

DA0A is the angle of attack Increment in degrees for the angle of 
attack sweep. 

A0AMAX is the maximum angle of attack in degrees for an angle of attack 
sweep run. 

PSI0 is the initial yaw angle in degrees for a yaw sweep run. 

DPSI is the yaw angle increment in degrees for the yaw sweep. 

PSIMAX is the maximum yav angle in degrees for a yaw sweep run. 

VWIND is the wind speed in ft/sec. Since the aerodynamics are incom- 
pressible and the coefficients are nondlmenslonal VWIND can be 
set at any nonzero value. 

Note that on one application of Option 11, either a yaw sweep or an 
angle of attack sweeps not both, can be run. 

4.2.12 Option 12 - Trim Data Summary 

Option 12 controls the printing of appropriate data from a series of 
trim conditions. The data to be Included are designated by the list of 
output parameters in the BDL and are found on the local file specified as 
NTAPE(3) in the Option 8 input data. Hie formats for printing the data on 
each page (a maximum of 12 pages is allowed) of the data summary are loaded 
by Option 7. The format for the headings on each page is specified in the 
Option 12 input. The two cards required : to Initiate Option 12 are as 
follows: 


Card 

Format 



1 

115 

1-5 

6-80 



N0PT 

Blank 
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Card 


Format 


2 8A10 


where 


1-80 

FMTl 


NOPT ■ 12 indlcatiag Option 12 

FHfl Is the format for printing the page headings for the trim data 

suonary. For VATL&S applications to date DCIl has been set at 
(7X, 11(2A5)). 
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S.O EXAMPLES OF VATLAS APPLICATIONS 


This section presents four examples vhlch demonstrate the application 
of most of the VATLAS options. Discussions of the examoles emphasise the 
Input card deck arrangement and the format of the subsequent printed output 
and plots. Discussion of the example outputs as they relate to the SF-121 
baseline flight control system development is provided in Volume I. The 
four examples and the Options they apply are as follows: 

o Example 1 - Transition Trim Run (Basic Data Load plus Options 1, 2, 

6, 7, 8, 9, and 12) 

o Example 2 - Stationkeeping Time History Run (Basic Data Load plus 
Options 1» 2 , 6, 8, and 9) 

o Example 3 - Wind Tunnel Data Run (Basic Data Load plus Options 2, 

7, and 11) 

o Example 4 - Statioidceeplng Plot Run (0; cions 5 - called implicitly 
by Option 10 - and 10) 

An annotated listing of the Basic Data Load cards for the application of 
VATLAS to the SF-121 airplane is provided in Appendix C. Only the cards 
which follow Che Basic Data Load in the example input decks will be presented 
in the exanq>le discussions. 

5.1 Exanq>le 1 - Transition Trim Run 

The transition trim run demonstrates the application of VATLAS to 
calculate a series of trim points along a wings level, constant altitude 
transition of the SF-121 airplane. The trims are at four airspeeds - 200 kt. 
180 kt, 160 kt. and 140 kt - in that order. An annotated listing of the 
input deck is shown in FI re 5-1. Observe that Options 6 and 9 must precede 
Option 9 and that Options 7 and 8 must precede Option 12 in the input sequence. 
Also note the use of the SCALE parameter in the cards for Options 1 and 2 
to load or print the data in units different from those used in the program. 
Finally note the placement of Option 2 before Option 8 to establish starting 
values for the trim at 200 kt. If these values are too far from Che actual 
trim, a trim may not be obtained and the program will abort. In a multiple 
trim run the starting point for all cases after Che first is the lomiediately 
preceding trim point. Therefore large steps between two adjacent trim cases 
should be avoided. For example, Che program would likely abort if it is 
required to step from the 200 kt trim to a 60 kt trim. The user will have 
Co experiment to find starting values and trim increments which are appropriate 
for his application. 

A copy of the output of Example 1 is provided in Appendix D. Items of 
interest Include the playback of the input data for the Basic Data Load and 
each Option, the stability derivative listing by components, and the trim 
summary table. 
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Figure 5-1 • Listing of Example 1 Input Data Deck Without Basic Data Load (Sheet 2 of 5) 
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Figure 5-1. Listing of Example 1 Input Data Deck Without Basic Data Lo»J (Sheet 4 of 5) 




52 


Figure S-I. Listing of Exaiaple 1 Input Data Deck Without Basic Data Load (Sheet 5 of S) 




5.2 Exaople 2 • Stationkceplng Tlaa History Run 

Ths ststionkeeping tine history run dsnonstrstss ths sppllcstlon of 
VATLRS to cslculsts the tlas history of s pseudo— pilot— controlled turn over 
e spot while hovering In s 35 kt wind. An snnotsted listing of the Ivut 
deck is shown in Figure 5-2. The following Itsas should be observed: 

1. Options 6 snd 9 must precede Option 8 in the input data deck and all 
three Options must precede Option 4. 

2. The run uses pseudo-pilot-controlled stationkeeping logic thus the 
parameters AKU, AKV, AKH, AKXES, AKTES, and AKZBS (Segment 16, 
elements 7-12^ and ASKEEP (Segment 43, element 7) must be specified 
in the input deck before Option 4. Hote that several of these 
parameters were specified in the Example 1 input deck (Figure 5-*l); 
they were ignored, however, since the example did not Include a time 
history. 

3. An open- loop psuedo-pilot input to the pedal is used to rotate the 
airplane, thus the parameter ACMD(l) (Segment 43, element 1) and the 
entries in tables CMDT3 and TDfET (also in Segment 43) must also be 
specified before Option 4. Note that these parameters also had values 
but were Ignored in Example 1 (Figure 5-1) . 

4. Example 2 requires a trim in a wind, thus the appropriate set of trim 
variables (Table 4-2) are introduced by Option 6. Furthermore, the 
parameters AWIND, VWIND, and PSIWIND (Segment 41, elements 1-3) must 
be specified before or %d.th Option 8. 

5. Since Example 2 is a time history, appropriate values for the integra- 
tion control parameters in COMK)N Segment 1 must be specified before 
invoking Option 4. These parameters had values in Example 1 (Figure 
5-1) but were not required and thus ignored. 

A copy of output of Example 2 is provided in Appendix E. The infor- 
mation is similar to that of Example 1 with the notable addition of the 
time history printout. The time history data were stored during the 
Example 2 run on local file TAPE4. This file was CATALOG ed as a permanent 
file by the control cards for the run. The file was ATTACHed in Example 4 
to demonstrate the use of the VATLAS plotting capability. 


5.3 Example 3 - Wind Tunnel Data Run 

The wind tunnel data run demonstrates the application of VATLAS to 
calculate six degrees of freedom stability axis wind tunnel data for the 
SF-12J airplane. Example 3 directs the execution of four wind tunnel runs: 
two on each of two different airplane configurations. For each configuration 
an angle of attack sweep from 0 to 90 deg in 5 deg Increments at zero yaw 
angle and a yaw sweep from 0 to 90 deg in 5 dag incrancnts at zero 
angle of attack are made; these runs require separate calls to Option 11. 

The configuration was set up by application of Option 2: the first configura- 

tion allows for the determination of pitch control effectiveness at symmetrical 
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Figure 5-2. Listing of Example 2 Input Data Deck Without Basic Data Load (Sheet 2 of 5) 
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Figure 5-2. Listing of Example 2 Input Data Deck Without Basic Data Load (Sheet 4 of 5) 
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Figure 5-2. Listing of Example 2 Input Data Deck Without Basic Data Load (Sheet 5 of 5) 




10 deg trailing edge up (TEU) eleven define tlona with the canard trailing 
edge flap set at 25 deg and the wing leading edge flap at 30 deg. The 
second configuration aalntains the sane canard and wing flap settings but 
allows for the detemlnatlon of roll control effectiveness at differential 
eleven deflections of 10 deg TEU on the left wing and 10 deg trailing edge 
down (TED) on the right wing. 

The input data deck (minus the Basic Data Load) for Bxaoq>le 3 is listed 
with explanatory marginal notes on Figure 5-3. Note that Option 7 must 
precede Option 11 and that the Option 7 variable formats should be considered 
unchangeable for most wind tunned applications of VATLAS. The order of the 
six events directed by Example 3 is as follows: 

1. Set up pitch control effectiveness configuration (Configuration 1) 
with Option 2. 

2. Run angle of attack sweep on Configuration 1 with Option 11. 

3. Run yaw sweep on Configuration 1 with Option 11. 

4. Set tg> roll control effectiveness configuration (Configuration 2) 
with Option 2. 

5. Run angle of attack sweep on Configuration 2 with Option 11. 

6. Run yaw sweep on Configuration 2 with Option 11. 

The airspeed parameter on the third Option 11 card must be given a positive 
uon-zero value for the Option to work. Since the Option generates non- 
dimensional aerodynamic coefficients and because the VATLAS aerodynamic 
model is not a function of airspeed, the selected airspeed value is not 
critical; Example 3 used a value of 100 ft/ sec. 

The output generated by VATLAS for Example 3 is listed in Appendix F. 

Note that the time history output specifications were inadvertantly left in 
the Basic Data Load for Example 3. These constitute extraneous input cards 
and could have been replaced by one blank card (i.e.. a time history was 
not generated by Example 3). 

5.4 Example 4 - Stationkeeping Plot Run 

The plot run demonstrates the application of VATLAS Option 10 and. by 
implication. Option S to plot a time history which was created in a previous 
run (Example 2). The input data deck to produce the plot is listed with 
explanatory notes on Figure 5-4. Two features of this deck idiich warrant 
special emphasis are the following: 

1. Since Options 10 and 5 are the only VATLAS options used in Example 
4. the C(M10N variables are not addressed during the run. Thus 
only five cards are required for the Basic Data Load. VATLAS requires 
that at least one CCBMON segment be set up by the Basic Data Load 
thus the CONHON storage allocation card for Example 4 designates one 
segment with one element (any number of elements could have been 
designated) . 
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LISTING OF exAHPue 3 INFUT DATA OSCK WITHOUT BASIC DATA LOAD 
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Figure 5-3. Listing of Ezajaple 3 Input Data Deck Without Basic Data Load 
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Figure 5-4, Listing of Example 4 Input Data Dock 






2. Not all of the variables stored during the tine history run need be 
plotted. In Example 4, only 28 of the 82 available variables are 
plotted. The ordering in the variable list for Option 5 (10) is that 
the first variable is the independent variable and all remaining 
variables are plotted as functions of that variable. In Example 4, 
therefore, 27 varlebles are plotted versus time. 

Although not indicated in the input data deck, the dimensioning and 
spacing of plots for Example 4 were keyed to 11 x 17 in. blank CALCOMP 
paper. This choice of paper is indigenous to VATLAS. Should the user desire 
different paper he must indicate this to the computer operator by a PAUSE 
card in the control cards for the run. Also not indicated In the input data 
deck is the process of making the standard CALCOMP subroutines available for 
the VATLAS plot routines. This is contrwlled by the control cards for the 
run. The control cards used to produce Example 4 on the Vought CDC 6600 
Coiq>uter Facility are listed below: 

CONTROL CARD 


FORTEN, TlOa 

PROJECT, 1234xxxx9101, 253010, 
S1128*FORTENBAXH . 

ATTACH, TAPE4, STATIONKEEPING, 
ID-FORTENBGH, MR-1. 

ATTACH, ONE, BINARYVATLAS, ID-FORTENBGH, 
MR-1. 

LIBRARY (CALCOMP) 


LOAD (ONE) 


EXECUTE. 


7 

8 
9 

The user is cautioned against indiscriminate application of these control 
cards to attempt to reproduce the Example 4 results since control card formats 
vary slightly from facility to facility. 

The plotted results of Example 4 are shown in Figure 5-5. These arc 
photo-reduced copies of the original 11 x 17 output plots. To learn how to 
use the VATLAS plotting capability, the user is urged to study carefully the 
relations between plot spacing, dlmenblonlng and labelling on Figure 5-5; 
the Option 5 instructions in Section 4.2.5; and the Example 4 input data 


EXPLANATORY NOTES 
Job Card 
Project Card 

Attach previously run time history 
Attach binary version of VATLAS 


Make CALCOMP routines available 
to VATLAS 

Load the binary version of VATLAS 
preparatory to executing the 
program. 

Execute the prograir and create 
the plot. 

End of control card file. 
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Figure 5-5. Plot Generated by Example 4 of "Pceudo-Ptlot Flown Turn 
Over a Spot In a 35 kt Wind" 


(Sheet 1 of 5) 
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deck. For reference scaling in this comparison, the ordinate axes of the 
original plots are 2.0 In. long and the abscissa (tine) axes are 8.0 In. 
long. 
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APPENDIX A 

DEFIMlTKMi AND LOCATKNI OF COMHON PARAMETBBS 
FOR THE VATLAS PROGRAM 


A-1 



A • 1 IntroducCidii 


This Appendix provides a listing by coded naiae and defintion of all 
VAILAS CCMON paraoieters arranged in order by segment and element designation. 
Each parameter is indicated a type C or N where C means the parameter is 
calculated by the program and does not need to be quantified by program input 
and N means the parameter is not calculated and must be quantified by program 
input. Hany of the parameters in CC^OION are also used in Volume 11 of this 
report which describes the model equations and base aircraft data for VATOL 
simulation. The Volume II symbology for these parameters is Included 
in parentheses following the definitions. Note that every parameter in 
COMMON is a REAL program variable; therefore all non zero input data to COMMON 
must have a decimal point. 


A.2 Segment 1 - Time, Print and Plot Control Parameters 


Coded 

Segment Element Type Name 


Units 


Definition 


1 

1 


1 


1 

2 


3 


C T 

C TO 

N DT 


sec Time which is updated every DT 
seconds (t) 

sec Time controlled by RUNGE which 

takes on three intermediate 
values between T and T+DT. Since 
4th order Runge - Kutta integra- 
tion is used, the basic integra- 
tion Interval must be divided in 
four parts. 

sec Basic integration interval for 

RUNGE 


1 


4 


N IMAX sec Maximum value of I for a time 

history 


1 

1 


5 

6 


N DTPRNT sec Interval between times of printed 

output 

N DTPL^T sec Interval between times of plotted 

output 


1 7 N TI 


sec Initial time for time history; 
normally taken to be zero 


A-2 


A. 3 Segment 2 - Variables to be Integrated 


Coded 


Segment 

Element 

Type 

Name 

Italts 

Definition 

2 

1 

C 

UD^T 

ft/aec^ 

X body axis component of inertial 
acceleration (u) 

2 

2 

C 

VDI^T 

ft/sec^ 

Y body axis component of Inertial 
acceleration (v) 

2 

3 

C 

WDI^T 

ft/sec^ 

Z body axis component of inertial 
acceleration (v) 

2 

4 

C 

PDCn 

rad/Mc^ 

Body axis roll acceleration (p) 

2 

5 

C 

QDdi 

rad/ sec^ 

Body axis pitch acceleration (q) 

2 

6 

C 

BS0T 

2 

rad/ sec 

Body axis yaw acceleration (f) 

2 

7 

C 

Dili) 

l/aec 

Rate of change of 
cosine 

dj^ 2 ^ direction 

2 

8 

C 

D120 

1 /sec 

Rate of change of 
cosine (d^ 2 ^ 

dj ^2 direction 

2 

9 

C 

D13D 

l/sec 

Rate of change of 
cosine (^ 2 ^ 3 ) 

dj ^2 direction 

2 

10 

C 

D21D 

1 /sec 

Rate of change of 
cosine (“ 22 ^) 

d^^ direction 

2 

11 

C 

D22D 

l/sec 

Race of change of 
cosine (^ 22 ) 

d 22 direction 

2 

12 

C 

D23D 

l/sec 

Race of change of 
cosine (^ 23 ) 

d 23 direction 

2 

13 

C 

D31D 

l/sec 

Rate of change of 
cosine (^ 32 ^) 

d^j^ direction 

2 

14 

C 

D32D 

l/sec 

Rate of change of 
cosine (^ 32 ) 

d ^2 direction 

2 

15 

C 

D33D 

l/sec 

Rate of change of 
cosine (^ 33 ) 

d ^3 direction 

2 

16 

C 

xedc(t 

£t/sec 

X Inertial axis component of 
Inertial velocity (x^) 

2 

17 

C 

YEDJJT 

ft/sec 

Y Inertial axis component of 
Inertial velocity (y^) 


A-3 



Coded 


Segment 

Element 

Type 

Name 

Units 

Definition 

2 

18 

C 

ZED0T 

ft/sec 

Z inertial axis component of 
inertial velocity (i^) 

2 

m 

19 

C 

LELEVND rad/ sec 

Left elevon deflection rate 

2 

20 

C 

RELEVND rad/ sec 

Right elevon deflection rate 

2 

21 

C 

DRD 

rad/ sec 

• 

Rudder deflection rate 

2 

22-23 

C 

XABD(I) 

1/sec 

Normalized rate of change of. 
afterburner thrust; 1*1 for left 
engine, *2 for right engine (T^(l)) 

2 

24-25 

C 

PSITD(I) 

rad/ sec 

Yaw thrust deflection rate; 1*1 
f^r left engine, *2 for right engine 

2 

26-27 

C 

THTD(I) 

rad/ sec 

Pitch thrust deflection rate; 1*1 
f§r left engine, *2 for right engine 

(•tW) 

2 

28-37 

C 

DRCSDCI) 

l/sec 

Normalized area deflection rate of 
1th RCS jet; space for 10 RCS . 

Jets, only 2 used for VATLAS (i^Qg(l)) 

2 

38-47 

C 

FRCSAD(l) lb/ sec 

Rate of change of force of 1th RCS 
Jet; space for 10, RCS jets, only 2 
used for VATLAS (F^^g (1)) 

A 

2 

48 

C 

DFLAPHD 

rad/ sec 

Canard, trailing edge flap deflection 
rate ) 

n 

2 

49 

C 

DFLAPWD 

rad/ sec 

Wing leading edge flap deflection 

2 

50-51 

C 

imCD(l) 1/sec 

Normalized rate of change of non- 
afterbuming thrust; 1*1 for left 
engine, *2 for right engine (Tp(l)) 

2 

52 

C 

PCTBMD 

rad/ sec 

Rate of change of roll command 
path integrator output (figure 2.2) 

(Pp ) 

''TRM 

2 

53 

C 

PINTD 

rad/sec 

Body axis roll rate; also coded as 


P In program (p) 



Coded 


Segment 

Element 

_ Type. 

Name 

Units 

Definition 

2 

54 

C 

PETSMD 

rad/sec 

Rate of change of roll forward 
path integrator output (figure 2.2) 

(}, ) 

'tsm 

2 

• 

55 

C 

RCTSKD 

rad/ sec 

Rate of change of yaw command path 
integrator output (figure 2.4) 

> 

*^TRM 

2 

56 

C 

BIKTD 

rad/ sec 

Body axis yaw rate; also coded as 
R in program (r) 

2 

57 

C 

RETBMD 

rad/ sec 

Rate of change of yaw forward path 
integrator output (figure 2.4) 

^^e > 

®TRM 

2 

58 

C 

ytrmd 

rad/ sec 

Rate of change of yaw trim Input 
(figure 2.4) (r^j^) 

2 

59 

,C 

QCTSMD 

rad/ sec 

Rate of change of pitch comsand path 
intagrator output (figure 2.3) 

(9c ) 

2 

60 

C 

QINTD 

rad/ sec 

Body axis pitch rate; also coded 
as Q in program (q) 

2 

61 

C 

QETBMD 

rad/ sec 

Rate of change of pitch forward 
path Integrator output (figure 2.3) 

2 

62 

C 

PTBMD 

rad/ sec 

Race of change of pitch trim input 
(figure 2.3) 

2 

63 

C 

ZEDT8MD 

Ib/sec 

Race of change of heave forward 
path Integrator output (figure 2.5) 
(T- ) 

TRM 

2 

64 

C 

THETD 

rad /sec 

• 

Rate of change of pitch angle (0^ 

2 

65 

C 

PHID 

rad/ sec 

• 

Race of change of roll angle (0) 

2 

66 

C 

PSID 

rad/ sec 

s 

Rate of change of yaw angle (90 

2 

67 

C 

LNGSTKD 

1/sec 

Rate of change of pitch input shaping 
filter output (q^) 

2 

68 

C 

ROLLCD 

1/sec 

Rate of change of roll input shaping 
filter output (p^) 

2 

69 

C 

YAWCD 

. 1/sec 

Rate of change of yaw input shaping 
filter output (r/) 
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A. 4 Sftgneat 3 - Int^rated Varlutdea 

Coded 


Segment 

Element 

■ST9I 

Name 

Units 

Definition 

3 

1 

c 

u 

ft/sec 

X body axis component of Inertial 
velocity (u) 

3 

2 

c 

V 

ft/aec 

Y body axis coBq>onent of Inertial 

• 





/eloclty (y) 

3 

3 

c 

w 

ft/sec 

Z body axis component of Inertial 
velocity (w) 

3 

4 

c 

p 

rad/ sec 

Body axis roll rate (p) 

3 

5 

c 

Q 

rad/ sec 

Body axis pitch rate (q) 

3 

6 

c 

R 

rad/ sec 

Body' axis yaw rate (c) 

3 

7 

c 

Dll 

— 

dj^j^ direction cosine (dj^j^) 

3 

8 

c 

D12 

— 

direction cosine 

3 

9 

c 

D13 

— 

d^^ direction cosine (d^^^) 

3 

10 

c 

D21 

— 

d 2 j^ direction cosine (^£ 1 ^ 

3 

11 

c 

D22 

— 

d 22 direction cosine (<* 22 ^ 

3 

12 

c 

D23 

— 

direction cosine (^ 22 ) 

3 

13 

c 

D31 

— 

d^^ direction cosine (d^^) 

5 

14 

c 

D32 

— 

dj 2 direction cosine (dj 2 ) 

3 

15 

c 

D33 

— 

djj direction cosine (d^^) 

3 

16 

c 

XE 

ft 

X position of eg In inertial 

coordinates Cx ) 

' e 

3 

17 

c 

YE 

ft 

Y position of eg In inertial 
coordinates (y^) 

3 

18 

c 

ZE 

ft 

Z position of eg in Inertial 
coordinates (z^) 

3 

19 

c 

LELEVN 

rad 

aft eleven deflection (1)) 

TEFy 

3 

20 

c 

RELEVM 

rad 

Right eleven deflection (f-pp (2)) 

^ W 

3 

21 

c 

DR 

rad 

Rudder deflection ( 


A-6 



Se^menc 

Element 

Type 

Coded 

Name 

Units 

Definition 

3 

22-23 

C 

TAB(I) 

— 

Normalized afterburner thrust; 1 • 
1 for left engine, *2 for right 
engine (T^(I)) 

• 3 

• 

24-25 

C 

PSITCI) 

rad 

Yaw thrust deflection; I*1 for 
left engine, •■2 for right engine 
(•Kjd)) 

3 

26-27 

C 

THT(I) 

red 

Pitch thrust deflection; I«1 for 
left engine, -2 for right engine 
(9^(1) ) 

3 

28-37 

C 

DRCS(I) 


Normalized area of Ith RCS Jet; 
space for 10 RCS jets, only 2 
used for VATLAS 

3 

38-47 

C 

FRCSA(I) 

lb 

Force of 1th RCS jet; space for 
10 RCS jets, only 2 used for 
VAIUS (P^ ) 

A 

3 

48 

C 

DFLAPU 

rad 

Canard trailing edge flap deflection 

(Itef^^ 

3 

49 

C 

DFLAPU 

rad 

Wing leading edge flap deflection 

3 

50-51 

C 

TFRAC(I) 

•• 

Normalized non-afterbumlng thrust; 
1-1 for left engine, *2 for right 
engine (Tp) 

3 

52 

C 

PCTRM 

rad 

Roll command path Integrator output 
(figure 2.2) (p ) 

TRM 

3 

53 

C 

PINT 

rad 

Integral of body axis roll rate 
(figure 2.2) (Pjj,^) 

3 

54 

C 

PETRM 

rad 

Roll forward path’ Integrator 
output (figure 2.2) (p ) 

®TRM 

3 

55 

C 

RCTBM 

rad 

Yaw command path Integrator output 
(figure 2.4) (r ) 

'^TRM 

3 

56 

C 

RINT 

rad 

Integral of body axis yaw rate 
(figure 2.4) 

3 

57 

C 

RETRM 

rad 

Yaw forward path integrator output 
(figure 2.4) (r ) 

®TRM 

3 

58 

C 

YTRM 

rad 

Yaw trim input (figure 2.4) (r.ppjj) 

3 

59 

C 

QCTRM 

rad 

Pitch command path Integrator output 
(figure 2.3) (qcxRM^ 
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Coded 

Segaent Element Type Name 


Ohlts 


Definition 


3 

61 

C 

QETRM 

rad Pitch forward path integrator out- 
put (figure 2.3) (q ) 

®TEM 

3 

62 

C 

FIRM 

rad Pitch trim Input (figure 2.3) (<l^;|y|) 

• 

3 

63 

C 

ZEDTBM 

lb Heave forward pach Integrator output 

(figure 2.5) (T ) 

T8M 

3 

64 

C 

THETA 

rad Pitch Angle (6) 

3 

65 

C 

PHI 

rad Roll Angle (0) 

3 

66 

• C 

PSl 

rad Yaw angle (^ ) 

3 

67 

C 

LN6STK. 

- Output of pitch ''nput shaping filter 

3 

68 

C 

ROLLC 

— Output of roll input shaping filter 

(P4> 

3 

69 

C 

TAHC 

Output of yaw Input shaping filter 

(r') 

A. 5 Segment 4 - 

Aircraft 

Hass, 

Geometry and Inertia Constants 

Segment 

Element 

Type 

Coded 

Name 

Units Definition 

.4 

1 

N 

WT 

lb Aircraft weight (W) 

4 

2 

N 

IX 

2 

slug ft Aircraft roll body axis Inertia (I^) 

4 

3 

M 

lY 

2 

slug ft Aircraft pitch body axlr Inevc’.a (ly) 

4 

4 

N 

IZ 

2 

slug ft Aircraft yaw body axis inertia (1^) 

4 

5 

N 

IXZ 

2 

slug ft Aircraft cross product of inertia (l^^^ 

4 

6 

. N 

FSCG 

ft Fuselage station of reference eg 

4 

7 

N 

BLC6 

ft. Buttllne of reference eg (BL^q) 

4 

8 

N 

■ULCG 

ft. Waterline of reference eg 

4 

9 

N 

GC 

2 

ft/sec Acceleration due to gravity (g) 

4 

10 

C 

EAC 

— Inert! ratio - (IZ - lY) » X 

4 

11 

C 

FAC 

Inertia ratio - (lY - IX)/’^Z 

4 

12 

c 

GAC 

— Inertia - ••io - (IZ - IX) /lY 

4 

13 

c 

JX 

— Cross Product of inertia ratio ■ 

- IXZ/ IX 
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Coded 


Segment 

Element 

Type 

Name 

QalCs 

Definition 

♦ 

14 

C 

JY 


Cross Product of inertia ratio * 
- IXZ/IT 

4 

15 

C 

JZ 

— 

Croaa Product of inertia ratio • 
- IZZ/IZ 

4^ 

Id 

c 

HASS 

slug 

Aircraft mass (m) 

4 

17 

c 

^Hjxjz 

— 

1 - JX*JZ 

A. 6 Segment 5- Tolerances for Errors In Trim Solutions 

Segment 

Eleim' t 

type 

Coded 

Name 

Units 

Definition 

5 

1-19 

N 

1^(1) 

Varies 

Specifies the acceptable toletance 
on the trln solution of the 1th 
trim equation. The trim equations 
are given in Subroutine ERBOR. 

A-7 Segment 6 - 

Aero Li^ 

t and Drag Constants 

Segment 

Element 

Type 

'led 

Marne 

Units 

Definition 

6 

1 

N 

AAEB0 



Specifies whether aero data are to 


be calculated <AAER^ = 0) or are 
stored in table format (AAE^ / 0). 


Note: For all sebscrlpted constants in Segment 6; I » 1 correspo' 

to wing data, I « 2 t*) horizontal stabilizer (includes fotw^, 
canard-andnft-conventional-locatlons) , I * 3 to vertical 
stab 11 izer. 


6 

2,25,48 

N 

AER^(1,I) 


Geometric or reference aspect ratio 
of Ith lifting surface (AR) 

6 

3,26,49 

N 

AER^(2,1) 

— 

Aspect ratio of exposed area of Ith 
lifting surface (ARe) 

6 

4,27,50 

N 

aee^O.i) 

ft^ 

Geotjetric or reference area of 1th 
lifting surface (S) 

6 

5,28,51 

N 

AE8J|(4,I) 

ft^ 

Exposed area of Ich lifting surface 
(Se) 

6 

6,29,52 

N 

A»!R#(5,I) 

rad 

Sueep of geometric or reference 
quarter chord of Ith lifting surface 

6 

7,30,53 


AER0(6,I) 

rad 

Sweep of leading edge of 1th lifting 
surface (Aj^) 

6 

8,31,54 

N 


— 

Reference tapei ratfo of 1th lifting 
surface (X) 
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Coded 


Seoeent 

Eleaent 

_Tj®« 

Naae 

Udlts 

Definition 

6 

9.3Z.55 

N 

AEI0(8,I) 

— 

Mazlnua lift coefficient of 1th 
lifting surface (Ci^y) 

6 

• 

10.33,56 

N 

AE80(9,I) 

rad 

Angle of attack for nxlmin lift 
coefficient of Ith lifting surface 

<V 

6 

11,34,57 

N 

AB80(1O,I) 

l/rad 

Sectional lift curve slope of Ith 
lift Surface (• 

o 

6 

12,35. -8 

M 

AEB^(ll.I) 

— 

Zero lift drag coefficient of Ith 
lifting surface (Cd^) 

6 

13,36,59 

H 

AEB^C12.I) 


Ratio of fuselage width at lifting 
surface intersection to lifting 
surface span for Ith lifting surface 
Cd/b) 

6 

14,37,60 

N 

AE8^(13,I) 

l/rad 

Change in lift coefficient at zero 
angle of attack per unit trailing edge 
flap deflection for Ith lifting surface 

6 

15,38,61 

N 

AER<^(14.I) 

l/rad 

Change in lift coefflclrat at angle 
of attack for naxtaua lift per unit 
trailing edge flap deflection for 
Ith lifting surface. ^ ^ ^ lEF^ 

6 

16,39,62 

N 

AER^(15,I) 

l/rad 

Change in drag coefficient per unit 
trailing edge flap deflection for 
Ith lifting surface (4k C^/ 

6 

17.40,63 

H 

AESil(16,l) 

l/rad 

Change in f itching nonent coefficient 
per unit trailing edge flap deflection 
for Ith lifting surface (ftC /^^c) 

6 

18.41,64 

a 

AER#(17,I) 


Change in angle of attack for maxlonim 
lift per unit leading edge flap 
deflection for Ith lifting surface 

MAX^^LEF^ 

6 

19,42,65 

H 

AE8^(18,I) 

l/rad 

Change in lift coefficient at zero 
angle of attack per imlt leading 
edge flap deflection for ith lifting 

surface /^LEp) 

o 

6 

20,43,66 

M 

AE8^(19,I) 

l/rad 

Change in lift coefficient at angle 


of attack for maxlaua lift per unit 
leading edge flap deflection for Ith 
lifting surface 
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Coded 

Segment Element Type Naae 


Units 


Definition 


6 21,44.67 M AEB^(20,1) 1/rad 

t 22,45.68 N AeM( 21,I) 1/rad 

6 23,46,69 1, AB8^(22,1> rad 

6 24,47,70 N AEB^(23,I) rad 

6 71,74,77 N SURFCPCl.I) ~ 

6 72,75,78 N SURFCPC2,I) -- 

6 73,76,79 N SURFCP(3,I) — 

6 80-82 C CSAAO(l) 1/rad 


Change in drag coefficient per unit 
leading edg'* flap deflection for 
Ith lifting surface ( 

Change In pitching acment coefficient 
per unit leading edge flap deflection 
for Ith lifting surface 

Difference between angle of attack 
for zero leading edge flap effectiveness 
and angle of attack for maxlimm lift, 
l.e. when (ALFA.(X. (AER0(22,I) 
AEB0(9,I)) the leading edge flap 
effectiveness equals zero (gp< 

Difference betwe^ angle of attack 
fo..- zero trailing edge flap effective- 
ness and angle of attack for naxiaun 
lift, l.e. when (ALFA.GE. (A£RO(23,l) 

+ AER0(9,I)) the trailing edge flap 
effectiveness equals zero (gdl^p) 

Noraallzed (to Ith lifting surface 
aean aerodynamic chord) const int for 
calculating center of pressure shift 
of Ith lifting surface as a function 
of angle of attack. Appears as 
Bultlplier of cosine terms in the 
equations (CP^7 

MoraaliTed (to Ith lifting surface 
aean aerodynamic chord) constant for 
calculating center of pressure shift 
of Ith lifting surface as a function 
of angle of attack. Appears as 
oultiplier of sine terms in the 
equations (CP^) 

Normalized (to Ith lifting surface 
mean aerodynamic chord) constant 
for calculating center of pressure 
shift of Ith lifting surface as a 
function of angle of attack. Appears 
only in the high angle of attack 
equation for center of pressure shift 
(CP3) 

Constant calculated by Subroutine 
FORCES for use by Subroutine LIFDRG 
to generate nonlinear Ith lifting 
surface lift coeff icient(C^, ) 
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Coded 

Segaent Element Type Name Ualts Definition 


6 

83-85 

C 


1/rad 

Constant calculated by Subroutine 
FORCES for use by Subroutine LIFDRG 
to generate nonlinear Itb lifting 
surface lift coefficient (K) . 

6 

86-88 

C 

E#S(I) 


Oswald efficiency factor for 1th 
lifting surface calculated in 
Subroutine F(ttCBS (e) 

6 

89-91 

C 

CLA(l) 

1/rad 

Linear lift curve slope of 1th 
lifting surface calculated in Sub- 
routine FORCES (C^ ) 

e 

A. 8 Segaent 7 - 

Wing Anro Constants 



Segaeiit 

Eleaent 

Type 

Coded 

Haw 

Units 

Definition 

7 

1 

■N 

FSW 

ft 

Fuselage station of wing reft^rence 
center of pressure (?S^) 

7 

2 

N 

Bun. 

ft 

Buttline of reference center of 
pressure of left vring (Bly(l)) 

7 

3 

N 

B1W2 

ft 

Buttline of reference center of 
pressure of right wing (BI^(2)) 

7 

4 

N 

WLW 

ft 

Waterline of wing reference center 
of pressure 

7 

5 

N 

CW 

ft 

Mean aerodynamic chord of wing (c^) 

7 

6 

N 

BW 

ft 

Reference wing span (b^) 

7 

7 

N 

EKBAl 


Change in left %rtng anple of attack 
per unit deflection oi left wing 
control surface (K^ (l)i 

7 

8 

M 

EK0A2 


Change in right wing angle of attack 
per unit deflection of right wing 
control surface (K r (2)) 

* c 

7 

9 

N 

EIW 

rad 

Wing incidence angle 

7 

10 

N 

CLBOW 

1/rad 

Basic C| of wing (C| ) 

r P 0 

7 

11 

N 

CLBaW 

1/rad 

Change in wing Cj|a wit >f 

wing lift coefficient ) 
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Coded 


Seijnent 

Elenent 

Type 

Marne 

Units _ 

Definition 

7 

12 

M 

CLRCLH 

1/rad 

Change in wing C|^ per unit of 
wing lift coefficient (C^ r^ ) 

7 

13 

N 

CLPW 

1/rad 

Wing Cj (C<_) 
P 

7 

14 

N 

CMBOW 

1/rad 

Basic C of vlng (CL ) 

7 

15 

N 

CMBCL2tf 

1/rad 

Change in wing C^^ per unit of 
irlng lift coeffic&nt squared. 

7 

16 

N 

CNRCL2W 

1/rad 

Change in wing Cq^ per unit of wing 
lift coefflclenc squared (<^ 2^ 

7 

17 

N 

CNRCDtf 

1/rad 

Change in wing per unit of wing 

drag coefficient ) 

7 

18 

N 

CNPCLW 

1/rad 

Oiange in wing Cq- per unit of wing 
lift coefficient (Cji_ ) 

PCl 

7 

19 

N 

CNPCDAW 

— 

Not used 

7 

20 

M 

(MOW 

— 

Zero lift pitching ooment coefficient 
of wing 

7 

21 

N 

CMDIW 

— 

Net used 

7 

22 

N 

AEKEL 

— 

Nuni>er of entries in elevon 
effectiveness table (^7) 

7 

23 

N 

CYBCL2W 

1/rad 

Change in wing Cy. Per unit of wing 
lift coefficirat squared .,) 

7 

24 

C 

XtfINGl 

ft 

X body axis distance between aircraft 
reference eg and left wing reference 
center of pressure 

7 

25 

c 

mMG2 

ft 

X body axis distance between aircraft 
reference eg and right wing refe'ence 
center of pressure 

7 

26 

c 

YWINGl 

ft 

Y body axis distance between aircraft 
reference eg and left wing reference 
center of pressure C'tvflNG^^^^ 
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Coded 

Segment Element Type Naie 


Units 


Definition 


7 

27 

C 

TUING2 

ft 

T body axis distance between aircraft 
reference eg and right wing reference 
center of pressure 

7 ; 

28 

C 

ZHIMGl 

ft 

Z body axis distance between aircraft 
reference eg and left wing reference 
center of pressure 

7 

29 

C 

ZWIMG2 

ft 

Z body axis distance between aircraft 
reference eg and right wing reference 
center of pressure (ZyjjjQ(2)) 

A. 9 Segment 8 - 

Horizontal Stabilizer Aero 

Constants 

Segment 

Elesmnt 

T?pe _ 

Coded 

Name 

Units 

Definition 

8 

1 

N 

FSK 

ft 

Fuselage station of horizontal 
stabilizer reference center of 
pressure (FS{j) 

8 

2 

N 

Bia 

ft 

Buttline of horizontal stabilizer 
reference center of pressure (Bl^) 

8 

3 

K 

von 

ft 

Waterline of horizontal stabilizer 
reference center of pressure (WI^) 

8 

4 

N 

CH 

ft 

Mean aerod)mainic chord of horizontal 
stabilizer (Cjj) 

8 

5 

N 

BH 

ft 

Reference span of horizontal 
stabilizer (by) 

8 

6 

N 

EKASe 


Change in horizontal stabilizer angle 
of attack per unit deflection of 
horizontal stabilizer control 
surface (1^^ ) 

8 

7 

C or N 

EIH 

rad 

Horizontal stabilizer Incidence angle 

(i«) 

8 

8 

M 

E|^CLW 

rad 

Change in downwash an,.le per unit 
of wing lift coefficient (6/Cj^y) 

8 

9 

N 

CMDFH 

— 

Not used 

8 

10 

N 

am 

— 

Zero lift pitching moment coefficient 
of horizontal stabilizer (Cm ) 

8 

11 

N 

DXHTO 

— 

Not used 
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Coded 


Segment 

Element 


Name 

Units 

Definition 

8 

12 

c 

XH 

ft 

X body axis distance between aircraft 
eg and horizontal stabilizer ref- 
erence center of pressure (X^) 

8 . 

13 

c 

YH 

ft 

Y body axis dlstauce between aircraft 
eg and horizontal stabilizer 
reference center of pressure (Y^) 

8 

14 

c 

ZH 

ft 

Z body axis distance between aircraft 
eg and horizontal stabilizer reference 
center of pressure (21^) 

8 15 
8 16 
8 17 
8 18 

A. 10 Segment 9 - 
Segmen' Element 

c 

c 

c 

c 

Vertical 

Type 

XMEL ft 

ZWEL ft 

AUEL rad 

DUEL — 

Stabilizer Aero 
Coded 

Name Units 

Constants related to relative wlng- 
horlzontal stabilizer location and 
which are used to calculate dynamic 
pressure at the horizontal stabilizer 

^*WEL’ ^L* ^L’ 

Constants 

Definition 

9 

1 

N 

FSV 

ft 

Fuselage station of vertical stabilizer 
reference center of pressure (FS^) 

9 

2 

N 

BLV 

ft 

Buttline of vertical stabilizer 
reference center of pressure (BL^) 

9 

3 

N 

WLV 

ft 

Katzrline of vertical stabilizer 
reference center of pressure (Wly) 

9 

4 

H 

CV 

ft 

Hear aerodynamic chord of vertical 
stabilizer (cy) 

9 

5 

N 

BV 

ft 

Reference span of vertical stabilizer 
(by.) 

9 

6 

N 

EKADR 


Change in vertical stabilizer angle 
of attack per unit deflection of 
vertical stabilizer control surface 

9 

7 

C or N 

EIV 

rad 

Vertical stabilizer incidence angle (i^) 

9 

8 

M 

EKV 

— 

Sldewash coefficient (K ) 

V 

9 

9 

N 

CNOV 

— 

Zero alee force yawing moment co- 
efficient of vertical stabilizer (C ) 

n 


o 

V 
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Coded 


Seenent 

Elenent 


Name 

Units 

Definition 

9 

10 

N 

AEERDD 

— 

Nuober of entries In rudder effective- 
ness table (£7) 

9 

• 

11-20 

N 

TALFA(l) 

rad 

Angles of attack for entires in 
vertical stabilizer effectiveness 
table (I-l to AALFA) 

9 

21 

N 

AALFA 

— 

Number of entries In vertical 
stabilizer effectiveness table ( £ 10) 

9 

22-31 

N 

TVTEFFCt) 


Vertical stabilizer effectiveness 
table entries (I ■= 1 to AALFA) 

^^3aER0 

9 

32 

C 

XV 

ft 

X body axis distance between aircraft 
eg and vertical stabilizer reference 
center of pressure (Xy) 

9 

33 

C 

YV 

ft 

Y body axis distance between aircraft 
eg and vertical stabilizer reference 
center of pressure (^y) 

9 

34 

c 

ZV 

ft 

Z body axis distance between aircraft 
eg and vertical stabilizer reference 
center of pressure (Zy) 

A. 11 Segmeat 10 - 

I^selage Aero Constants 


Segaent 

Elenent 

lype 

Coded 

Name 

Units 

Definition 

10 

1 

N 

FSF 

ft 

Fuselage station of fuselage 
reference center of pressure (FS^) 

10 

2 

N 

BLF 

ft 

Buttllne of fuselage reference 
center of pressure (BLp) 

10 

3 

N 

WLF 

ft 

Waterline of fuselage reference 
center of pressure (WLp) 

10 

4 

N 


ft^ 

Maxlnpua fuselage axial force normal- 
ized to dynamic pressure((A/q)„. „) 

riAA 

10 

5 

N 


ft 

Maximum fuselage side force normalized 
to dynamic pressure((S/q)j^) 

10 

6 

N 

NF|fQM 

ft^ 

Maximum fuselage normal force normalized 
to dynamic pressure ((**f/9)j^) 

10 

7 

N 

tffQMXL 

ft^ 

Maximum fuselage pitch moment In 
region ( If /1 8 + AOF) i ot ^ AIF - 
normalized to dynamic pressure 




A-ie 



Coded 

Segment Element Type Name Units Definition 


10 

8 

N 

M0QMX2 

ft^ 

Naximum fuselage pitch moment In 
region - AIF ^ - normalized 

to dynamic pressure 

10 . 

9 

N 

H0QHU 

ft^ 

Maximum fuselage yaw moment In region - 
1T/18 i *^Fy * ~ tiojmalized to 

dynamic pressure where *^Fy ** 
fuselage sideslip angle ((«/q)nAXi^ 

10 

10 

N 

H0QHX2 

ft^ 

Maximum fuselage yaw moment In region - 
A2F'< - normalized to 

dynamic pressure ((N/q)„.^ ) 

10 

11 

N 

AOF 

rad 

Angle of attack where fuselage normal 
force equals zero ( Op) 

10 

12 

M 

AIF 

rad 

Angle of attack where fuselage pitch 
moment eqxials zero (^ip) 

10 

13 

N 

A2F 

rad 

Fuselage sideslip angle where fuselage 
yaw moment equals zero 

10 

14 

N 

ENl 

— 

Exponent of cosine shaping function 
for fuselage axial force (n^) 

10 

15 

H 

EN2 

— 

Exponent of sine shaping function for 
fuselage side force (n2) 

10 

16 

N 

EH3 

— 

Exponent of sine shaping function for 
fuselage normal force (n^) 

10 

17 

N 

ENA 


Exponent of sine shaping function for 
fuselage pitch moment in the region 
( It'/lS + AOF) 6 bl. • AIF (n^) 

10 

18 

N 

EN5 



Exponent of sine shaping function for 
fuselage pitch moment in the region 
(n^) 

10 

19 

N 

EN6 


Exponent of sine shaping function for 
fuselage yaw moment in the region 
ir/18 £ «Py ^ A2F (n^) 

10 

20 

H 

EN7 


Exponent of sine shaping function for 
fuselage yaw moment in the region 
A2F £ ^ Fy ^^7^ 

10 

21 

M 

DFSFO 

ft 

Maximum aft shift of fuselage center of 
pressure with angle of attack ( FSp^) 
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Segment 

Element 

Tyjje 

Coded 

Name 

Units 

Definition 

10 

22 

N 

ADFSFl 

rad 

Angle of attack where fuselage center 
of pressure begins to shift aft 
(Afsj^) 

10 

23 

N 

ADFSF2 

rad 

Angle of attack tdiere fuselage center 
of pressure shift equals DFSFO 

(AfS2^ 

10 

24 

N 

XKXF 

" 

Constant to control fuselage side 
force contribution to fuselage yaw 
moment (Kyj.) 

10 

25 

C 

XF 

ft 

X body axis distance between aircraft 
eg and fuselage reference center of 
pressure (Xp) 

10 

26 

C 

YF 

ft 

Y body axis distance between aircraft 
eg and fuselage reference center of 
pressure (Yp) 

10 

27 

c 

ZF 

ft 

Z body axis distance between aircraft 
eg and fuselage reference center of 
pressure (Z^) 


A. 12 Segment 11 - Direct Thrust Force and Moment Constants 


Coded 

Segment Element Type Same Units Definition 


11 

1-2 

N 

FSSW(l) 

ft 

Fuselage station of engine nozzle 
swivel point; 1*1 for left engine, 
“ 2 for right engine (FSsy(I)) 

11 

3-4 

N 

BLSW(l) 

ft 

Buttline of engine nozzle swivel 
point; 1*1 for left engine, ■ 2 
for right engine (BLgyf(I)) 

11 

5-6 

N 

WLSW(I) 

ft 

Waterline of engine nozzle swivel 
point; 1*1 for left engine, * 2 
for right engine (WLgwt^)) 

11 

7 

N 

SIGY 

rad 

Inclination of engine centerline to 
fuselage refererce line ( 

11 

b 

N 

ELN 

ft 

Length of engine nozzle (1 {jq2) 

11 

9 

N 

KAFTl 

1/rad 

Change of engine nozzle thrust 
efficiency per unit thrust turning 


angle 
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Coded 


Segment 

Element 

Type 

Name 

Units 

Definition 

11 

10 

N 

KAFT2 

— 

Engine nozzle thrust efficiency 
at zero thrust turning angle 

11 

11 

N 

AENG 

— 

Number of engines ( Sfe 2) (n ) 

11 

12 

C 

CSIGY 

— 

cos ( y) 

11 

13 

C 

SSIGY 

— 

sin ( <T y) 

11 

14-15 

c 

XT (I) 

ft 

X body axis distance between 


aircraft eg and nozzle swivel point; 
1*1 for left engine, » 2 for right 
engine (Xi]^(I)) . 


11 16-17 C YT(I) ft Y body axis distance between aircraft 

eg and nozzle swivel point; 1*1 for 
left engine, * 2 for right engine 
(Yj(D) 

11 18-19 C ZT(I) ft Z body axis distance between aircraft 

eg and nozzle swivel point; 1*1 for 
left engine, ■ 2 for right engine 
(Z^(D) 

A. 13 Segioent 12 - Inlet Ram Force and Moment Constants 

Coded 

Segment Element Type Mane Units Definition 


12 

1-2 

N 

FSIN(I) 

ft 

Fuselage station of engine inlet force; 
1*1 for left engine, * 2 for right 
engine (FSjjj(I)) 

12 

3-4 

N 

BLIN(I) 

ft 

Buttline of engine centerline at 
inlet face; 1*1 for left eBgine, 
- 2 for right engine (BLj^(I)) 

12 

5-6 

N 

WLIN(I) 

ft 

Waterline of engine centerline at 
Inlet face; 1*1 for left engine 
■ 2 for right engire (WL^jj(I)) 

12 

7 

N 

RTNLET 

ft 

Equivalent inlet radius (r^j) 

12 

8 

N 

A\.'V1 

— 

Hunger of entries in V/Vj table (6 5) 

12 

9 

N 

AATURM 

— 

Number of entries in ram moment arm, 
ram efficiency and incremer “^al turning 


angle tables ( £ 8) 
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Coded 

Segment Element Type Mane 


Units 


Definition 


12 

10-11 

C 

xmcD 

ft 

X body axis distance between aircraft 
eg and engine inlet; 1*1 for left 
engine* ■ 2 for right engine (Xjjj(I)) 

U 

12-13 

C 

YIN (I) 

ft 

Y body axis distance between aircraft 
eg and engine inlet; 1*1 for left 
engine* * 2 for right engine (Yj|)(I)) 

12 

14-15 

C 

ZIN(I) 

ft 

Z body axis distance between aircraft 
eg and engine Inlet; 1*1 for left 
engine* ■ 2 for right engine (Zuj(I)) 

A. 14 Segment 13 - 

Coriolis Force 

and Moment Constants 

Segment 

Element 

Type 

Coded 

Name 

Units 

Definition 

13 

1 

N 

ELDUCT 

ft 

Length of engine duct (Ipuc-j) 

A. IS Segment 14 - 

Propulsion System Constants 

Segment 

Element 

Type 

Coded 

Name 

Units 

Definition 

14 

1 

N 

AMACU 

— 

Number of entries in Mach number 
table ( A 4) 

14 

2 

N 

KBTl 

Ill 

Fractional change in thrust per 
percent bleed 

14 

3 

N 

KBT2 

— 

Fractional thrust at zero bleed (Rgj^^ 

14 

4 

M 

KMl 


Fractional change in inlet mass flow 
rate per fractional change in thrust 

14 

5 

N 

KM2 

— 

Fractional inlet mass flow race at 
zero fractioru'i thrust (K„ ) 

14 

6 

H 

ilMEGMX 

rad /sec 

Maximum engine rotational speed in 
rad/sec 

14 

7 

N 

RPMMX 

rpm 

)iaximum engine rotational speed in 
revolutions per minute 

14 

8 

H 

RCS 


Indicates whether reaction control 


system is active; RCS i 0, RCS is 
active; -• 0, RCS not active 
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Segment 

Elamnt 

Type 

Coded 

Nom 

Oalts 

Definition 

14 

9 

N 

AFRAT 

— 

Nunf>er of entries In fractional thrust 
table for determining engine rotational 
speed ( ^ 6) 

14 

10 

N 

JENG 

slug £t^ 

Engine rotational inertia per engine 
^■^ENG^ 

14 

11 

N 

AFRAC 


Number of entries in fractional thrust 
table for engine dynamics parameters 
( 6 10) 

14 

12 

N 

XAUAB 

sec 

Time constant for afterburner thrust 
dynamics ( 

14 

13 

N 

ABACCEL 

1/sec 

Acceleration limit of afterburner 
thrust 

14 

14 

N 

ABDECEL 

1/sec 

Deceleration limit of afterburner 
thrust (TAijjy,) 

14 

15 

N 

TAB^ 


level of TABIN when afterburnei: la 
lit'(T ) 

^ON 

14 

16 

N 

TABOFF 


Level of TAB (I) when afterburner is 
turned off ,T^ ) 

“off 

14 

17 

C 

KBT 

— 

Fractional thrust level at reference 
bleed (K^^) 


A. 16 Segment 15 - RCS Force and Moment Constants 


Coded 

Segment Element Type Name Units Definition 


15 

1 

N 

ARCS JET 

— 

Uund>er of RCS Jets ( 10) 

15 

2 

N 

ABLDR 

— 

Number of entries in bleed tables (^6) 

i5 

3 

N 

BLDBEF 

% 

Reference bleed (B^^p) 

15 

4-13 

N 

FRCSMX(I) 


Maximum force of Ith Jet (^rcSmx^^^^ 

15 

14-23 

V 

IMSMAMD(l) 


Indicates whether Jet I is demand 
bleed; DEMAND(I) ■ 1, jet I is demand; 

■■ 0, Jet I is continuous bleed tDMD(I)) 

15 

24-33 

N 

BLDM^R(l) 


Indicates whether Jet I requests 
additional bleed; BLDM0R(1)-1, Jet 1 
requests additional bleed; 0, Jet 

1 must operate with available bleed 
(BLDM't)) 
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Coded 

Seitae; . Elecent Type Name 


Unite 


Definition 


15 

34-83 

N 

SCSL(1,J) 


Definee whether Jth control drives 
Ith jet; RCSL(I,J) ^ 0, Jth control 
drives Ith jet; ■ 0, Jth control does 
not drive Ith jet. J ■ 1 roll 
control; J ■ 2, pitch control; J > 3 
yav control; J 4, side force control 
J ■ 5» normal force control 

15 

84-93 

N 

psfcscD 

ft 

Fuselage station of Ith jet (FSjg^(I)) 

15 

94-103 

N 

BLRCS(I) 

ft 

Buttline of Ith jet (BLjg^(I)) 

15 

104-113 

N 

WLRCSCI) 

ft 

Vaterline of Ith jet (WLr:.^(I)) 

15 

114-123 

N 

THTRCS(I) 

rad 

Pitch orientation of jet thrust line 
(ejET(i)) 

15 

124-133 

N 

PSIB.CSCI) 

rad 

Yaw orientation of jet ihrus" line 

15 

134 

C 

FRCSMXR 

— 

Mot used 

15 

135 

C 

FRCSMXO 

lb 

Total RCS force produced at reference 


conditions (reference bleed and 
zero control inputs) ; If all jets are 
demand bleed, then FRCSMXO > 0 


(FrcSmxo> 


15 

136 

C 

FRCSMXl 

— 

Mot usee 


15 

137-146 

C 

CTHTRCSCD 

— 



15 

147-156 

C 

STHTRCS(I) 

— 

•in (9 jei(D) 


15 

157-166 

C 

CPSIRCS(I) 

— 



15 

167-176 

C 

SPSIRCS(I) 

— 

.In 


15 

177-186 

C 

XJxI(l) 

ft 

X body axis distance between 
eg and Ith let (Xj^^'I)) 

aircraft 

15 

187-196 

C 

YJET(I) 

ft 

Y body axle distance between 
eg and Ith jet (Xjg^(I)) 

aircraft 

15 

197-206 

C 

Z JET CD 

ft 

Z body ' Is distance between 
eg and i.ch jet (Zj^^(I)) 

aircraft 
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A. 17 Segment l6 - Control Systoi Constants 


Coded 


.egment 

Element 


Name 

Uhits 

Definition 

16 

1 

N 

DUM 

— 

Eot used 

16 

2 

1! 

AKZED 

— 

Rot used in final version of program. 

io 

3 

N 

AKZEEl 

— 

Rot used in final version of program. 

16 

U 

N 

AKZEE2 

— 

Rot used in final vctrsion of program. 

16 

5 

N 

AKYED 

— 

Hot used in final version of program. 

1 ; 

# 

H 

AKYEE 

— 

Rot iised in final version of program. 

26 

7 

N 

AKU 

sec 

Pilot gain on ZEDOT feedback for 
stationkeeping option 

Id 

8 

N 

AKV 

sec 

Pilot gain on YHXJT and XEDOT feedback 
in lateral posicion loop for station- 
keeping option 

16 

9 

N 

AKW 

sec 

Pilot gain on XEDOT and YEDOT feed- 
back in longitudinal position loop 
for stationkeeping option. 

16 

10 

N 

AKXES 

l/sec 

Pilot heave rate controller gain for 
stationkeeping option 

16 

11 

N 

AKYES 

rad/ft 

Pilot lateral stick gain for station- 
keeping option 

16 

12 

N 

AKZES 

rad/ft 

Pilot longitudinal stick gain for 
stationkeeping option 

16 

13 

N 

CSV 

— 

Hot used in final version of program. 

l6 

lU 

N 

VASW 

ft/sec 

Control system function switch speed 

<’'a > 

SB 

16 

15 

If 

CSl 

— 

Rot used in final version of program. 

16 

16 

N 

THETSW 

— 

Rot used in final version of program. 

l6 

17-2U 

N 

TAKPC(I) 

rad 

Roll control system input gain table 
if the gain is a function of pitch emgle 

16 

25-32 

N 

TTHETA(I) 

rad 

Pitch cuigle table for determining 
control systm gains. 

16 

33-hO 

N 

TAKPE(I) 

1/rad 

Roll control sj'^stem error gain table if 
the gain is a function of pitch emgle. 

l6 

ia-U8 

N 

TKPCIl(l) 

rad/sec Roll control system input integrator 


gain table if the gain is a function 
of pitch angle. 

6 - 1 


Coded 


SeicMnt 

Elencat 


Ham 

Ubtta 

Definition 

16 

49 

N 

A3HE1A 

— 

Musbar of entries in control systaa 
gains pitch angle table ( ^ 8) 

16 

50-57 

H 

TKPBUCD 

rad/cAc/ 

rad 

Soil control syatus error integrator 
gain table* 

16 

50-65 

H 

TKP(l) 

rad/rad/ 

sac 

Boll control syaten roll rate 
feedback gain table* 

16 

66 

H 

PKLIM 

— 

Output linit on roll control 
systen error integrator (P«j^jn) 

16 

67 

8 

AKOToa 

rad 

Boll control to yaw control inter- 
connect gain (R^ j ) 

16 

68 

M 

CIHTR 

rad 

Pitch thrust deflection to roll 
control gain g ) 

16 

69 

M 

KOA 

rad 

Differential wing trailing edge 
flap (elevon) deflection to roll 
control gain (R# ) 

® A 

16 

70-77 

H 

XAKRC(l) 

rad 

Yaw control systen input gain cable* 

16 

78-85 

M 

TAKBB(l) 

1/rad 

Yaw control systen error gain table* 

16 

86-93 

H 

IKECUd) 

rad/ sac 

Yaw control systen input integrator 
gain table* 

16 

94-101 

N 

TKBEll(l) 

▼ad/ sac/ 
rad 

Yaw control systen error integrator 
gain table* 

16 

102-109 

N 

TKR(l) 

rad/ rad/ 
sac 

Yaw control syaten yaw race feedback 
gain table* 

16 

110 

M 

8ELIN 

— 

Output Unit on yaw control system 
error integrator (*^enn) 

16 

lU 

N 

AiaCTUf 

rad/ sec 

Yaw control systen trimer gain 

16 

112 

H 

KPSIT 

rad 

Yaw thrust deflection to yaw control 
gain (R^.j ^ y) 

16 

113 

H 

KDR 

tad 

Rudder deflection to yaw control 
gain 

16 

114-121 

H 

TAKAYCD 

rad/ ft/ 
sec2 

Lateral acceleration feedback gain 
table; loaded as a function of 
TQBAR(l) 


•Provided the pnin is a ttinot. it-n 
nf pitch angle. 

A-24 


Halts 


Dsfinltlon 


Coded 

Segeeet Eleeent Type Met 


16 

m 

H 

KIT 

lb 

Thrust connsid to yaw control gain 0^^) 

16 

U3-130 

M 

TAKQCCD 

rad 

Pitch control system input gain table* 

16 

Ul-138 

M 

xyava) 

rad/eec 

Pitch control system error gain, 
table* 

16 

139-1*»6 

H 

TKiCnd) rad/sec 

Pitch control syston ininit integration 
gain table* 

16 

147-154 

M 

TKOBUCD rad/sec/ 
rai 

Pitch control ays ten error integrator 
gain table* 

16 

155-162 

H 

TKKD 

rad/rad/ 

sec 

Pitch control system pitch race feed- 
back gain table* 

16 

163 

H 

QKLIM 

— 

Output limit on pitch control system 
error Integrator 

16 

164 

N 

AKPTKM 

rad/ sec 

PiC^ control sys.tam trimser gain 

16 

165 

N 

DSTK 

— 

When the absolute value of longitud- 


Inal stick deflection is greater 
than DSTK, autoaatic pitch tr inning 
is activated 


16 

166 

M 

KDE 

rad 

Symmetrical wing trailing edge flap 
Celevon) deflection to pitch control 
gain (K^^) 

16 

167 

N 

KTHTP 

rad 

Pitch thrust deflection to pitch 
control gain (Kg^4 p) 

16 

168-175 

H 

XARZE(l) 

Ib/ft/sec 

Heave control system error gain table* 

16 

176-183 

N 

IKZEIl(l) 

Ib/sec/ 

ft/sec 

Heave control system error integrator 
gain table* 

16 

184 

N 

AKIHB^T 

lb 

Thrust command to manual throttle 
input sun O^ggoj) 

16 

185 

N 

ZELIN 

lb 

Output limit on heave control system 
error integrator 

16 

186 

N 

ELEVNKF 

rad 

Reference elevon deflection (^xgpvj^) 

16 

187-188 

N 

PSITRF(I) 

rad 

Reference yaw thrust deflection; 
1*1 for left engine, "2 for right 
engine 

16 

184 

M 

KHFLAPl 

•— 

Change in horisontal stabilizer 


trailing edge flap deflection per unit 
change in angle of attack 


•Provided the gain is a function of 
pitch angle. 


Coded 

Segwttt Eleeeat Type Ma— Cdlf 


16 

190 

M 

XHPUP2 

rad 

16 

191 

H 

DiVLAPl 

— 

16 

192 

N 

IHFIAra 

rad' 

16 

193 

H 

AQEAR 

— 

16 

19‘»-201 

M 

TORARCt) 

Ib/ft^ 

l6 

202-209 

R 

TTPrrcB(i) 

sec 

16 

210-217 

R 

TTROLL(I) 

sec 

l6 

218-225 

R 

TTTAW(I) 

sec 


A. 18 Segpent 17 - Actuation Systee Constants 

Coded 


Segment 

E}.ement 


Maas 

Units 

17 

1 

H 

LBOIATL 

rad/ sec 

17 

2 

H 

LSLRATL 

rad/ sec 

17 

3 

H 

leuNsl 

rad 

17 

4 

H 

LELP^SL 

rad 

17 

5 

M 

TLEL 

sec 

17 

6 

N 

REURAXL 

rad/sec 

17 

7 

N 

RELBATL 

rad/ sec 

17 

8 

M 

REUP0SL 

cad 

17 

9 

H 

RE1JP9SL 

rad 

17 

10 

N 

TREL 

sec 

17 

11 

H 

DRURATL 

rad/ sec 

17 

12 

H 

ORLRATL 

rad/ sec 

17 

13 

H 

DRUPEL 

rad 

17 

14 

H 

DBLP0SL 

rad 

17 

15 

H 

TDR 

sec 


Definition 

Horlsontal stabiliser trailing edge 
flap deflection at leco angle of 
attack 

Change In wing leading edge flap 
deflection pec unit change In 
angle of attack 

King leading edge flap deflection 
at seco angle of attack 

Hunber of entries in dynanic pressure 
table for lateral acceleration 
feedback gain ( dk 8) 

Dynanic pressure t'^ble for lateral 
acceleration feedback gain 

Pitch input sbaping filter tine 
constant table if the tine constant 
is a function of pitch angle 

Roll input shaping filter time constant 
table if the time constant is a function 
of pitch angle 

Yav input shaping filter time constant 
table if the time constant is a function 
of pitch angle 


Definition 


Left wing trailing edge flap (elawon) 
actuator constants In order from 
element 1 to element 5: upper rate 

limit, lower cate limit, 'pper 
position limit, lower position limit, 
time constant 

Right wing trailing edge flap (eleven) 
actuator constants In order from 
element 6 to element 10: upper rate 

limit, lower rate limit, upper position 
limit, lower position limit, time 
constant 

Vertical stabiliser control surface 
(Rudder) actuator constants in order 
from element 11 to eleamnt 15: upper 

rate limit, lower rate lidit, upper 
position limit, lower position limit, 
time constant 
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Coded 

Segffont Element Type Na»e 


Units 


Defioltlop 


17 

16 

H 

HFURm 

rad/ sec 

17 

17 

N 

HFLRATL 

rad/ sec 

17 

18 

H 

HFLAPIC' 

rad 

17 

19 

H 

HPUPMr 

rad 

17 

20 

N 

THP 

sec 


17 

21 

K 

WFUSATL 

rad/sec 

17 

22 

H 

WFLBATL 

rad/ sec 

17 

23 

N 

tfFLAPHX 

rad 

17 

24 

M 

mAPMI 

rad 

17 

25 

H 

TWP 

sec 

17 

26-27 

M 

THTURLCI) 

rad/ sec 

17 

28-29 

N 

IHTLBL(I) 

rad/sec 

17 

30-31 

H 

THRUPL(I) 

rad 

17 

32-33 

N 

THTLPL(I) 

rad 

17 

34-35 

N 

TTHT(I) 

sec 


17 

36-37 

N 

PSIUFL(I) 

rad/ sec 

17 

38-39 

N 

PSILRL(I) 

rad/ sec 

17 

40-41 

N 

PSIUPL(I) 

rad 

17 

42-43 

M 

PSTT.PL(I) 

rad 

17 

44-45 

N 

TPSIT(I) 

sec 


17 

46-55 

N 

RCSDRL(I) ~ 

17 

56-65 

N 

RCSLRL(I) — 

17 

66-75 

N 

RCSUPL(I) — 

17 

76-85 

N 

RCSLPL(I) — 

17 

36-95 

N 

XAURCS(l) sec 


Borlsontal stabiliser trailing edge 
fli^ actuator constants In order froit 
eleaent 16 to elesant 20: upper rate 

lialt» lover rate llnlt» upper 
position lielt* lower position llodt, 
tlae constant 

Wing leading edge flap actuator 
constants in order from eleient 21 
to eleaent 25: upner rate Unit, 

lover rate limit, upper position 
Halt, lower position limit, time 
constant 

Pitch thrust deflection actuator 
constants in order (two elements' per 
constant where I " 1 is left engine 
and 1 « 2 is right engine) from 
element 26 to element 35: tq>per 

rate limit, lover rate limit, upper 
position limit, lower position liaiit, 
tiae constant 

Yaw thrust deflection actuator 
constants in order (two elements per 
constant where 1 » 1 is left engine 
and I ~ 2 is right engine) from 
eleaent 36 to mlement 45: upper 

rate limit, lover rate limit, upper 
position limit, lover position limit, 
time constant 

RCS area actuator constants in order 
by groups of ten (up to ten RCS jets 
can be specified) from elemeDt 46 to 
element 95: upper rate limi*: (elements 
46 to 55) , lower rate limit (elements 
56 to 65), upper position limit (elements 
66 to 75) , lower position llotlt 
(elements 76 to 85) , time constants 
(elements 86 to 95) 
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Coded 


Secment 

EleacnC 

Ki'JJi 

Naae 

Uaita 

Definition 

17 

96-lOS 

H 

VSCSl]BL(l) Ib/aec 

KCS" force dynaalcs conetanta In 

17 

106-U5 

H 

PBCSLILCI) Ib/aec 

order by groups of ten (up to ten 

17 

U6-125 

H 

FRCSDPLCI) lb/ 

RCS Jeta can be specified) froa 

17 

126-135 

H 

ISCSLPLCI) u> 

eleaent 96 to eleaant 145: i^per rate 

17 

136-145 

H 

t40FBCS(I) 

sec 

Halt (eleaents 96 to 105) » lower 
rate Halt (eleaents 106 to 115), 
upper position Halt (eleaents 116 
to 125), lower position Halt (eleaents 
]26 -to 135) , tlae constant (eleaeuta 
136 to 145) 

17 

146 

N 

ULIM 

rad 

MaxlauB coaaanded syaactrlcal wing 
trailing edge flap (eleven) 
deflection ( 

17 

147 

H 

DALIM 

rad 

Maxlaua cMaanded differential 
wing trailing edge flap (eleven) 






deflection ( 

17 

148 

N 

DBLIM 

rad 

Maxlaua coaaanded vertical stabilizer 
control surface (rudder) deflection 







A. 19 Segaeat 18 - 

Lift-Drag Variables 





Coded 



Segaent 

Bleaent 


hoae 

Halts 

Definition 

18 

1 

c 

LAMBl 



Mlscel\-'neou3 intemiedlate parameters 

18 

2 

c 

CLAE 

1/rad 

associated with calculaf s of 

18 

3 

c 

EKHB 

— 

CNAAO(l), AEROK(l), E0S(i>. and 

18 

4 

c 

a 

— 

CLA(l) . In Volume 11 symbology. 

18 

5 

c 

C2 

— 

these parameters are In order: 

18 

6 

c 

EJ 


taaA]^^2> » *W(B)» ^2 

and J ^ 

A.20 Segment 19 - 

tflng Aero Variables 





Coded 



Sesaent 

Elatatait 


Naas 

ftilts 

Definition 

19 

1 

c 

DXtf 

lb 

X body axis wing aerodynamic force 

(kV 

19 

2 

c 

DW 

lb 

Y body axis wing aerodynaalc force 
(hV 


19 


3 


DZH 


lb 


Z body axis win ' aerodyaaolc force 


0(*d«4 

Sttttt 11 — w t Ptftnitloo 


19 

4 

G 

SW 

ft lb 

Wag aaro47BaaAc bod^ axia toll 
-D-mt (Al^ 

19 

S 

C 

aMT 

ft lb 

Wag aarodfaaale bodj axia pltcb 
MMt (AV 

19 

€ 

C 

UM 

ft lb 

Wag aarodvBaalc bodr axla jm 
mamat (Ai^) 

19 

7 

C 

cur 

— 

fatal wlag lift coafflclaat (Cy^^) 

19 

8 

C 

Aunr 

rad 

Aforaga local aagla of attack of 
Wag ( aliij) 

19 

9 

C 

cnr 

— 

fatal vlag drag coafflclaat (Cpg^) 

19 

10 

C 

w 

ft/aac 

X body axla alr^aad at wlag caatar 
of v<-'^aaata (oy) 

19 

11 

C 

HH 

ft/aac 

Z body axla alrapaad at wlag coiter 
of praaaara (a^ 

19 

12 

C 

V2Z1 

ft/aac 

Total alrapaad aaaaad by left Wag 
la aa X-Z body axla plane at the 
left Wag canter of preaaure 

(▼xz(i)) 

19 

13 

C 

VXZ2 

ft/aac 

Total alrapeed aenaed by right Wng 
la aa X-Z body axia plane at the 
right wing center of preaaure 
(Vn(2)) 

19 

14 

C 

VXZ 

ft/aac 

Average total alrapaad aenaed by 
left and rl|At Wag halvea (Vj2» ~ 
0.5 (Vjb(1) + Vjj,(2))) 

19 

15 

C 

nATUI 

rad 

Change la aagle of attack of left 
wing due to left wing' control surface 
deflection ( 

19 

16 

C 

nfPA7 

rad 

Change in angle of attack of right 
wing dua to right wing control 
outface deflection (Am A 

19 

17 

C 

4LIH1 

rad * 

locAl angle of attack of left wing 

(mJ(D) 

19 

18 

C 

Aina 

rad 

Local angle of attack of right wing 

19 

19 

C 

CLHl 

— 

Lift coefficient of left wing * 
(CLy(l)) 

19 

20 

C 

CUI2 

— 

Lift coefficient of right wing 
<CLm(2)) 
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Codtd 


S«cMnt 

BlMMt 



Ukiit* 

Dafinltlon 

19 

a 

C 

cncL 

— 

Drag eoafflclaat of laft wing 
(C^,(l)) 

19 

22 

C 

CW2 


Dcag eoafflclaat of right %ilng 

(Cdj,(2)) 

19 

23 

c 

CD4 

— 

Hot uaad 

19 

24 

c 

BBW 

rad 

Hlng aldaallp angle (|^) 

19 

25 

c 

QMS 

Ib/ft^ 

Dynamic preaaure baaed on VXZ(^) 

19 

24 

c 

QHSB 

Ib/ft 

qHS*Bfr 

19 

27 

c 

BOI? 

MC 

0.5 BU/VZZ(b/V) 

19 

29 

c 

IMP 

ft lb 

Wing stability axla roll aonent 

19 

29 

c 

NNP 

ft lb 

Wing stability axis yaw noaeat (N^ ) 

19 

30 

c 

YHP 

lb 

Wing aide force 

19 

31 

c 

ALFl 

rad 

Left wing angle of attack (Vy(l)) 

19 

32 

c 

SALFl 

— 

sin (ALFl) 

19 

33 

c 

CALFl 

— 

cos (ALFl) 

19 

34 

c 

ALF2 

rad 

Right %rlng angle of attack («(^C2)) 

19 

35 

c 

SALF2 

«ii* ^ 

sin (ALF2) 

19 

36 

c 

CALF2 

— 

cos (ALF2) 

19 

37 

c 

jori 

lb 

X body axis left wing aerodynaalc 
force (Iw(l)) 

19 

38 

c 

XM2 

lb 

X body right wing aerodynamic 

force (X^(2)) 

19 

39 

c 

ZHl 

lb 

Z body axis left wing aerodynaalc 
force (Z^(l)) 

19 

40 

c 

ZH2 

lb 

Z body axis right wing aerodynaalc 
force (Z|f(2)) 

19 

41 

c 

ALP 

rad 

Wing angle of attack (o(f ) 

19 

42 

c 

SALF 

— 

sin (ALF) 

19 

43 

c 

CALF 

— 

cos (ALF) 
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Segment 

Element 

TyM 

Coded. 

Name 

Units 

Definition 

19 

44 

C 

MI 

ft Ih 

Hlng aerodynamic pitch moment due 
to flaps and zero lift pitch 

19 

45 

C 

Ut 

ft lb 

Vtag body axis roll moment without 
aaynmatrlc lift effects (ly^) 

19 

46 

C 

NH 

ft lb 

Vlng body axis yaw moment without 
aayHMtric drag effects 

19 

47 

C 

DAI 

rad 

Left wing trailing edge flap 
deflection (elevon) corrected for 
large deflection effects (4^py(l)) 

19 

48 

C 

DA2 

rad 

Right wing trailing edge flap 
deflection (elevon) corrected for 
large deflection effects 

19 

49 

C 

EKBLEL 

— 

Large deflection correction factor 
for left wing elevon (K^pp^Cl)) 

19 

50 

C 

EKELER 

— 

Large deflection correction factor 
for right wing elevon 

19 

51 

C 

CMIl 


Left wing pitch moment coefficient 
due to flaps and zero lift pitch 

(c^(D) 

19 

52 

C 

CMH2 


Right wing pitch moment coefficient 
due to flaps and zero lift pitch 
(C^(2)) 

19 

53 

C 

DmNGl 

" 

Normalized (to CH)dhl£t of left wing 
center of pressure as a function of 
angle of attack ( 4 X^jjjgd)) 

19 

54 

C 

DXHIN62 


Normalized (to CW )shlf t of right 
wing center of pressure as a function 
of angle of attack (& 


A. 21 Segnent 20 - Horizontal Stabilizer Aero Variables 

Coded 

Segment Element Type Name Unlta Definition 


20 

1 

C 

DXE 

lb 

X body axis horizontal stabilizer 
aerodynamic force (AX^) 

20 

2 

C 

DYE 

lb 

Y body axis horizontal stabilizer 
aerodynamic force (^^Yh' 
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Coded 

Segiaent Element Type Name 


Unite 


Definition 


20 

3 

C 

OZH 

lb 

Z body axis horiasontal stabilizer 
aerodynasLc force (bZ^) 

20 

4 

C 

SLH 

ft lb 

Horizontal stabilizer body axis 
roll noaent ( 

20 

5 

C 

Di& 

ft lb 

Horizontal stabilizer body axis 
pitch noaent ( 

20 

6 

C 

dhh: 

ft lb 

Horizontal stabilizer body axis 
yaw moaent ( iiN^) 

20 

7 

C 

UH 

ft/sec 

X body axis airspeed at horizontal 
stabilizer center of pressure (u^) 

20 

8 

C 

MB. 

ft/sec 

Z body axis airspeed at horizontal 
stabilizer center of pressure (w^) 

20 

9 

C 

ALFG 

rad 

Aircraft angle of attack (e( ) « also 
coded as ALPH In segment 29 

20 

10 

C 

EPS 

rad 

Oownwash sncle ( ^ ) 

20 

11 

C 

OALFtf 

rad 

Incremented angle of attack for 
hacking angle of attack rate effects 
( or 2^ 

20 

12 

C 

CLl 


Wing lift coefficient at DALFW where 
a ^pMMktlve 0( increment Is used 
(Cli) 

20 

13 

C 

CL2 


Wing lift coefficient at DALFW where 
a negative o( increment Is used 

20 

14 

C 

CDl 


Wing drag coefficient at DALFW where 
a positive increment is used 

(Cdi) 

20 

15 

C 

CC2 


wing drag coefficient at DALFW where 
a negative o( increment is used 

(CP2) 

20 

16 

c 

OAK 

rad 

Change in horizontal stabilizer angle 
of attack due to horizontal stabilizer 
control surface deflection ) 

20 

17 

c 

ALFR 

rad 

Horizontal stabilizer local angle 
of attack (*<^) 

20 

18 

c 

CIH 

— 

Horizontal stabilizer lift coefficient 
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Coded 


Segment 

Element 

KT9I 

Name 

Units 

Definition 

20 

19 

c 

COB 


Horizontal stabilizer drag co- 






efficient (Cq^) 

20 

20 

c 

AN6 

rad 

Variables related to relative vlng- 

20 

21 

c 

DIS 


horizontal stabilizer location and 

20 

22 

c 

XI 

— 

which are used to calculate 

20 

23 

c 

ZU 


pressure at the stabilizer CA j. i>a ^ 

20 

24 

c 

DQoqo 

— 

Normalized dynamic pressure loss at 
horizontal stabilizer center of 






pressure ( h 9/5^) 

20 

25 

c 

ALFHl 

rad 

Horizontal stabilizer angle of attack 






C- tan”l (MH/uh)) 

20 

26 

c 

SALFHl 

— 

sin (ALm) 

20 

27 

c 

GALFHl 

— 

cos (ALIHl) 

20 

28 

c 

DAE 

rad 

Change in horizontal stabilizer 
angle of attack due to angle of 
attack rate 

20 

29 

c 


rad 

Horizontal stabilizer control surface 






(i.e. elevator) deflection 

20 

30 

c 

QEFF 

Ib/ft^ 

Effective dynamic pressure at 
horizontal stabilizer (q^pp) 

20 

31 

c 

DAV6 

rad 

Average deflection of left ^Uld right 






wing trailing edge flaps ( ^xef^) 

20 

32 

c 

Qfi 

— 

Horizontal stabilizer pitch moment 
coefficient due to flaps (Cm^) 

20 

33 

c 

DXHT 


Normalized (to CE)shlft of horizontal 
stabilizer center of pressure as a 
function of angle of attack (hX^^) 

A. 22 Segment 21 

- Vertical Stabilizer Aero 

Variables 




Coded 



Segment 

Element 


Marne 

Units 

Definition 

21 

1 

c 

DXV 

lb 

X body axis vertical stabilizer 
aerodynamic force (AXy) 

21 

2 

c 

DYV 

lb 

Y body axis vertical stabilizer 
aerodynaiulc force ( t^Y^) 
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Coded 

Segment Element Type Name 


Units 


Definition 


21 

3 

C 

DZV 

lb 

Z body axla vertical stabilizer 
aerodynamic force ( LZy) 

21 

4 

C 

DLV 

ft lb 

Vertical stabilizer body axis roll 
moment (&Ly) 

21 

5 

C 

im 

ft lb 

Vertical stabilizer body axis pltc'ti 
moment (A My) 

21 

6 

C 

BNV 

ft lb 

Vertical stabilizer body axis yaw 
moment (^Ny) 

21 

7 

C 

UV 

ft/sec 

X body axis airspeed at vertical 
stabilizer center of pressure (uy) 

21 

8 

C 

W 

ft/aec 

Y body axis airspeed at vertical 
stabilizer center of pressure (vy) 

21 

9 

C 

DADR 

rad 

Change in vertical stabilizer angle 
of attack due to vertical stabilizer 
control surface deflection 

21 

10 

C 

ALFV 

rad 

Vertical stabilizer local angle of 
attack ( y) 

21 

11 

C 

CLV 

— 

Vertical stabilizer lift coefficient 
CCLy) 

21 

12 

C 

CDV 

*— 

Vertical stabilizer drag coefficient. 
(Coy) 

21 

13 

C 

ALFVl 

rad 

Vertical stabilizer angle of attack 
(- tan-^ (W/UV))( 

21 

14 

C 

CALFVl 

— 

cos (ALFVl) 

21 

15 

C 

SALFVl 

— 

sin (ALFVl) 

21 

16 

C 

QOV 

Ib/ft^ 

Effective dynamic pressure at vertical 
stabilizer (qy) 

21 

17 

C 

DRP 

rad 

Rudder deflection corrected for 
large deflection effects ( j ^) 

21 

18 

C 

EKB1<0 


Large deflection correction factor 
for rudder (i^py ) 

21 

19 

C 

CMV 

— 

Vertical stabilizer pitch moment 
coefficient due to flaps (C_ ) 
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Coded . 


Segment 

Element 

Type 

Name 

Units 

Definition 

21 

20 

C 

DXVT 

— * 

Normsllced ^o CV)shift of vertical 
stabiliser center of pressure as a 
function of angle of attack ( 

21 

21 

C 

VTKFP 

— 

Vertical tall effectiveness factor 
^'^T£pfJ 

A.23 Segment 22 - 

Fuselage Aero Variables 


Segment 

Element 

Type 

Coded 

Name 

Ublta 

Definition 

22 

1 

C 

SEF 

lb 

X body axis fuselage aerodynamic 
force (^Xp) 

22 

2 

C 

DYP 

lb 

Y body axis fuselage aerodynamic 
force (bYp) 

22 

3 

C 

DZF 

lb 

Z body axis fuselage aerodynamic 
force itkZf) 

22 

4 

C 

DLF 

ft lb 

Fuselage body axis roll moment (AL^,) 

22 

5 

C 

DMF 

ft lb 

Fuselage body axis pitch moment ANp) 

22 

6 

C 

ONF 

ft lb 

Fuselage body axis yaw moment (i^Np) 

22 

7 

C 

ALFF 

rad 

Aircraft angle of attack («(), also 
coded as -ALPH In segment 29 and ALFG 
in segment 20 

22 

8 

C 

ALFYF 

rnd 

Fuselage yaw angle of attack 

22 

9 

C 

ALFTF 

rad 

Fuselase total anele of attack 
(- tan“l ( JvAs2 + WAS-^ VUAS))(«f^) 

22 

10 

C 

QXT 

Ib/ft^ 

Dynamic pressure based on UAS and 
VAS (- 1/2 ^(UAS2 + VAs2))(TJj^) 

22 

11 

C 

QO 

Ib/ft^ 

Total dynamic pressure (q) , also 


coded as QBAR in segment 29 
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Coded 


Seeaenc 

Element 

Type 

Hane 

Units 

Definition 

22 

12 

C 

QXZ 

lb/£t^ 

Dynamic pressure based oa UAS 

and H4S C* 1/2 f (UAS* + HAs2))Ci^) 

22 

13 

C 

A#Q 

ft^ 

Fuselage axial force coefficient CA/q) 

22 

14 

C 

AS 

lb 

Ftoelage axial force (Ap) 

22 

15 

C 

SAP 

— 

sin CAP) - AF is element 17 in 
Segm«it 22 

22 

16 

C 

S#Q 

ft^ 

Aiselage side force coefficient (S/q) 

22 

17 

C 

AF 

rad 

Angle redefined man^ times Hhile 
calculating fuselage aero forces 
and moments 

22 

18 

C 

MI]^ 

ft^ 

F^elage normal force coefficient 

22 

19 

C 


ft 

Fiiselage pitch moment coefficient 
(M/q) 

22 

20 

C 

*#ll 

ft^ 

Fuselage pitch moment coefficient 
wnen AF -If2/(18(A1F - AOF)) 

22 

21 

C 

s 

lb 

Fhselage side force (Sp) 

22 

22 

C 

NF 

lb 

Fuselage normal force (Mp) 

22 

23 

C 

PITCBF 

ft lb 

Fuselage pitch moment about reference 
eg (H) 

22 

24 

C 

ALFTFD 

— 

Not used 

22 

2S 

C 


ft^ 

Fuselage yaw moment coefficient (N/q) 

22 

26 

C 

H0Q1 

ft^ 

Fuselage yaw moment when 
AF - 

22 

27 

C 

TAHF 

ft lb 

Fuselage yaw wowent about reference 
eg (N) 

22 

28 

C 

DFSF 

ft 

Shift of fuselage center of pressure 
as a function of angle of attack 


C AFSp) 
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A.24 S«tMBt 23 - Direct Thxuet Force and KoMst Variables 

MOTE: For all sobscrlpted variables In segpent 23, 1*1 Indicates left engine 

and I«2 indicates right engine 

Coded 


Seeaenc 

Elenent 


Mane 

Units 

Definition 

23 

1 

c 

DXT 

lb 

X body axis direct thrust 

‘W 

23 

2 

c 

DTT 

lb 

T body axis direct thrust 
tore. 

23 • 

3 

c 

DZT 

lb 

X body axis direct thrust 
*»«• <*nsT^ 

23 

4 

c 

DLT 

ft/lb 

Direct thrust body axis roll 

23 

5 

c 

DMT 

ft/lb 

Direct thrust body axis pitch 

. 23 

6 

c 

DMT 

ft/lb 

Direct thrust body axis yaw 

‘"nsi’ 

23 

7-8 

c 

AFT(I) 

rad 

Flow turning angle of exhaust 

(Apj(D) 

23 

9-10 

c 

TCfril) 

lb 

Net thrust applied to aircraft 
after correction for flow 
turning losses (Tqqp(I)) 

23 

11-12 

c 

DX(I) 

lb 

X body axis force applied by Ith 
engine (AX(I)) 

23 

13-14 

c 

DT(I) 

lb 

T body axis force applied by Ith 
engine ( A T(D) 

23 

15-16 

c 

DZ(I) 

lb 

Z body axis force applied by Ith 
engine ( A Z(I)) 

23 

17-18 

c 

DFSSW(I) 

ft 

X body axis distance between 
nosxle swivel point and exit 
plane ( AFS^jCD) 

23 

19-20 

c 

DBLSW(I) 

ft 

T body axis distance between 
nosxle swivel point and exit 
plane ( A BLgvjd)) 

23 

21-22 

c 

PHLSV(I) 

ft 

Z body axis distance between 
nosxle swivel point and exit 
plane ( A WLsyd)) 

23 

23-24 

c 

SPSIT(I) 

- 

sin (PSIT(D)-PSIT(I) located 
In segment 3, elements 24 and 25 

23 

25-26 

c 

CPSIT(I) 

— 

cos (PSIT(D)-PSIT(I) located 
In segment 3, elements 24 and 25 

23 

27-28 

c 

srar(i) 


sin (THT(D) THT(I) located in 
segment 3, elements 26 and 27 

23 

29-30 

c 

CTHT(I) 

— 

cos (THT(I))-THT(I) located In 
segment 3, elements 26 and 27 

23 

31 

c 

KAPT 


Flow turning correction factor 

(Kaft^ 


A-37 


A.2S taipMt 24 - Inlet Baa Force and Monent Variables 

MOTE: For all aabacrlptad varlablao la aagpant 24, 1*1 Indlcacaa left aaslaa, 

1*2 ladlcacaa right angina. 


Coded 

Saint Elanant Mhlta Definition 


24 

1 

C 

DZl 

lb 

X body axis inlet ram force 

24 

2 

C 

DTE 

lb 

T body axis inlet ram force 

^\am^ 

24 • 

3 

C 

DZR 

lb 

2 body axis inlet ram force 

24 

4 

C 

DLR 

ft lb 

Inlet ram body axis roll OKmient 

24 

S 

C 

DMR 

ft lb 

Inlet ram body axis pitch mcmient 

24 

6 

C 

mat 

ft lb 

Inlet ram body axis yaw moment 

24 

7 

C 

VO 

ft/sac 

True airspeed (V^) , also coded 
as VA in segment 29 

24 

8 

C 

ATURNO 

rad 

Ceonatric inlet flow turning 
angle in X-Z inlet plana (A^ygu ) 

24 

9 

C 

BTDSM 

rad 

4J 

Geometric inlet flow turning angle 
in X-T inlet plane (b^ugn) 

24 

10 

C 

CBTURH 


cos (BTUBN) 

24 

11 

C 

SBTDRN 

- 

sin (BTURN) 

24 

12 

C 

VTN 

ft/s€C 

Inlet flow velocity (Vjj|) 

24 

13-14 

C 

V^VI(I) 

- 

Ratio of true air^eed to inlet 
flow velocity (VOV^d)) 

24 

15-16 

C 

ATUBM(I) 

rad 

Actual inlet flow turning angle 
in X-Z inlet plane (A.j^im(I)) 

24 

17-18 

C 

ELBAM(I) 


Inlet ram monent am normalized 

to equivalent inlet diameter 

24 

19-20 

C 

DATDRN(I) 

rad 

Difference between actucl and 
geometric Inlet flow turning 
angles ( ^ Ayur{((1)) 

24 

21-22 

C 

BTAR(I) 

— 

Ratio of actual to theoretical inlet 
ram force magnitude (9|g(I)) 

24 

23-24 

C 

FRAM(I) 

lb 

Actual inlet ram force magnitude {F_ (I)) 

RAM 

24 

25 

C 

ELR 

ft 

Inlet ram force moment aim (distance along 

engine centerline from inlet face 

to inlet ram force application point) 

24 

26 

C 

VLRAM 

ft 

Waterline of inlet ram force application 

point 
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Patf 


Peftnitlon 


Scgaent Elcweat 


type 


Coded 


24 

27 

C 

FSK4M 

ft 

Fdaelage etetion of inlet ram force 
a^llcetlan point (FSg^^ 

24 

28 

C 

SATDKH 

- 

Ueed to r^resent sin (ATOKHCI)) 
for ell velnee of I 

24 . 

29 

c 

CATD8H 

— 

Ueed to repreaent cos (ATDBll(I)) 
for all values of I 

24 

30-31 

c 

nno(i) 

lb 

X body axis inlet ram force due 

to freestreae air flow ( A 

24 

32-33 

c 

TSD(I) 

lb 

T body axis inlet ram force due 

to freestream air flow ( ^T|^(I)) 

24 

34-35 

c 

DZSO(I) 

lb 

Z body axis Inlet ram force due 

to freestream air flow ( ii Z|^(I)) 

24 

36-37 

c 

Dm(i) 

lb 

X body axis inlet ram force due 

inlet flow velocity Induced by 
aircraft rotation rates (Al^(I)) 

24 

38-39 

c 

DTRI(I) 

lb 

T body axis inlet ram force due 

inlet flow velocity Induced by 






aircraft rotation rates ( (I)) 

24 

40-41 

c 

DZRI(I) 

lb 

. body axis inlet ram force due 

inlet flow velocity induced by 
aircraft rotation rates ( ^ Z^(D) 

24 

42 

c 

UI 

ft/sec 

X body axis Inlet flo*' velocity 
Induced by aircraft rotation rates 

(Ui) 

24 

43 

c 

VI 

ft/sec 

T body axis l 'it flow velocity 
induced by aircraft rotation rates 

<»i) 

24 

44 

c 

VI 

ft/eec 

Z body axis inlet flow velocity 
induced by aircraft rotation rates 

A.26 Segeent 25 - 

Corlolla Force and Hoaent Verlebles 




Coded 



Segnent 

Eleaeat 

Type 

Naae 

Onlts 

Definition 

25 

1 

C 

DXC0k 

lb 

X body axis Coriolis force (Yqq|{^) 

25 

2 

C 

DTCAt 

lb 

T body axis Coriolis force (Yqq|^) 

25 

3 

C 

DZC0R 

lb 

Z body axis Coriolis force (Zqq^^) 

25 

4 

C 

DLC08 

ft/lb 

Coriolis body axis roll moment 

25 

5 

C 

1MC0R 

ft/lb 

Coriolis body axis pitch moment 

25 

6 

C 

lHfC#R 

ft/lb 

Coriolis body axis yaw moment (^qq^) 
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Uttlts 


Dtfinltlon 


Co4ed 

Seitot El—ent Type Ma»e 


25 

7 

C 

KLQ^ 

ft 

X distance fron aircraft eg to 
angina lalat face naasurad along 
angina cantarllne. 

25 

8 

• 

C 


ft 

Z dlstanca fron aircraft eg to 
angina eanterline naasured 
parallel to engine Inlet face. 

15 

9 

C 

DLp^tD 

ft/lb 

Coriolis roll nonent about engine 
centerline (1 qo|^) 

25 

10 

C 

OMC^ 

ft/lb 

Coriolis body axis pitch nonent; 
sane as DHC(Hl 

25 

11 

C 

DRq^ 

ft/lb 

Coriolis yaw nonent perpendicular 


to plane defined by engine 
centerline and T b^y axis 


A. 27 i>igBcnt 26 - Propulalon Syatcn Variables 

I'OTB: For all anbacrlpted variables in this aegnent, !•! indicates left engine, 

I*2 Indicates right engine. 

Coded 


S^«it 

BleMnt 

Type 

Naae 

Units 

Definition 

re 

1-2 

C 

MDOT(l) 

Ibn/sec 

Inlet nass flow rate (*(!)) 

26 

3-4 

C 

RPM(I) 

rev/aln 

Engine rpn 

26 

5-6 

C 

TRST(I) 

lb 

Engine thrust corrected for RCS 
bleed (T^pl(D) 

2‘ 

7 

C 

BLDAVL 

Z 

Bleed available 

26 

8 

C 

FGMAX 

lb 

Maxlnua afterburning thrust 

26 

9 

C 

FOIIN 

lb 

Maxlnua nun-afterbumlng thrust 

(Fn ). Also ninlnua afterburning tnrust 
HIN 

26 

10 

C 

FGIDL 

lb 

Idle thrust (Fp ) 
IDL 

26 

11 

C 

MDTHX 

Ibn/sec 

Maxlnua Inlet nass flow rate 

26 

12 

C 

TRCSREF 

— 

Not used 

26 

13 

C 

TREQAV 

lb 

Average uncorrected (for RCS and 
flow turning effects) thrust (T^a ) 

^AV 

26 

’4-15 

C 

TREQ(l) 

lb 

Uncorrected (for RCS and flow turning 
effects) thrust of Ith engine (T^(l)) 

26 

16 

C 

FRATIO 

- 

Ratio of uncorrected tlirust to 
aaxlnntn non-af ter burning thrust 

26 

17 

C 

FRPM 

- 

Fractional rpn (%p(|) 
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ttllf 


Patinltloii 


Coded 

Segamt El— tat Type Mat 


26 

26 

18 

19 

C 

C 

XMDT 

KBTP 

- 

Fractional Inlet mass flow rate, 
llitttst correction factor for RCS 
coapling (Kg^ 

26 

26 

• 

20 

21-22 

C 

C 

TRCS 

TC(I) 

lb 

lb 

Tbrust loan due to RCS coupling 
Coananded tbruat corrected for 
RCS reference bleed (T^CD/K^,) 

26 

23-24 

. C 

HE(I) 

lb ft sec 

Engine angular noaentua (R^(I)) 

26 

25-26 

C 

TFBACC(I) 


Ratio of TC(I) to maximum non- 
afterburning tbruat (Tp ) 

*nl 

26 

27 

C 

TF8AC1N 

• 

Uaed to repreaent TFRACC(I) for 
all valnea of I (Tp ) 

^IN 

26 

28 

C 

FRACL 


Ratio of waxlmum non-afterbuming 
tbruat to Idle tbruat (F« /F^ ) 

'*IDL 

26 

29 

C 

TEHG 

sec 

Engine time conatant (T£hq(I)) 

26 

30 

C 

TFSACB 

• 

Non-afterbuming thmst dynamics 

loop error (Tp ) 

E 

26 

31 

C 

TABIN 

- 

Input to afterburning thrust 
dynamics loop (T^n ) 

IN 

26 

26 

32 

33 

C 

C 

BL0T3 

TABE 

- 

Not used. 

Afterburning thrust dynamics loop 
error (T^ ) 

A 

26 

34-35 

C 

150(1) 

lb 

Sate as TREQ(I) (T (1)) 


A. 28 Segnent 27 - RCS Force and Moaent Variables 


Coded 

Segment Element Type Name Units Definition 


27 

1 

C 

DXRCS 

lb 

X body axis RCS force (Xrcs> 

27 

2 

C 

DTRCS 

lb 

T body axis RCS force (Yrcs) 

27 

3 

C 

DZRCS 

lb 

Z body sxls RCS force (Zrqs) 

27 

4 

C 

DLRCS 

ft/lb 

RCS body axis roll moment 

27 

5 

C 

DHRCS 

ft/lb 

RCS body axis pitch moment (M ) 

RCS 

27 

6 

C 

DNRCS 

ft/lb 

RCS body axis yaw moment (Np^g) 

27 

7 

C 

OFRCSDB 

lb 

RCS force produced by those demand 
jets which require additional 
bleed ( 4 F - 1 

27 

8 

C 

DFRCSD 

- 

Not used. 

27 

9-18 

C 

FRCS(I) 

lb 

Uncorrected RCS force produced 
at Ith jet (Fj^jjgd)) 
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Coded 


Sesiwat 

Rleaeat 

Type 

Marne 

Units 

Definition 

27 

19 

C 

BLDRBQ 

X 

Bleed retired by RCS (B^^q) 

27 

20 

C 

FRCSAL 

lb 

Total RCS force available within 






bleed Halts ) 

All 

27 . 

21 

C 

FRCSLD 

lb 

Total RCS force desired ) 

27 

22-31 

C 

TRCSCd) 

lb 

Coeasnded RCS force at Ith jet 
corrected for bleed available 
and MSS flow rate 

27 

32-41 

C 

DXRC(I) 

lb 

X body axis force produced by 
1th Jet ( ^ Xj^cs^D) 

27 

42-Si. 

C 

DTRC(I) 

lb 

Y body axis force produced by 
1th Jet ( A Trcs(I>' 

27 

52-61 

C 

DZRC(I) 

lb 

Z body axis force produced by 
1th Jet ( A ZrcsU)) 

27 

62-71 

C 

DRCSC(I) 

— 

CoBaanded notMlised area of 






1th RCS Jet ( i Rcs^(I)> 

27 

72 

C 

BID 

Z 

Bleed actually used by RCS (BIA) 

27 

73 

C 

FRCSDB 

• 

Not used. 

A. 29 Segaeat 28 - 

Control System Vsrlsbles 





Coded 



Senent 

Bleaent 

Type 

N«ae 

Units 

Definition 

28 

1 

C 

ZEE 

_ 

Not used . 

28 

2 

C 

ZEIC 

ft 

Initial value of ZE 

28 

3 

C 

TEE 

ft 

Error in pseudo-pilot closed sero 
lateral drift loop 

28 

4 

C 

YEIC 

ft 

Initial value of YE 

28 

5-8 

C 

DXE.DTE, 

DZE.XES 


Not used 

28 

9 

C 

YES 

ft 

Y body axis position error in 
pseudo-pilot stationkeeping option 

28 

10 

C 

ZES 

ft 

Z body axis position error in 
pseudo-pilot stationkeeping option 

28 

11 

C 

CSSU 

- 

Control system function switch 






gain (eSg^) 

28 

12 

C 


rsd or 

Cockpit input to roll control 





red /sec 

systea; coablnea lateral stick 
and pedal Inputs in accordance 
with control systea function 
switching 

28 

13 

C 

DLATSTK 

rsd or 
rsd/sec 

Lateral stick deflection ( ijj^xSTK,' 

28 

14 

C 

AKPC 

rsd 

Roll control systea Input gain(K 


A-A2 



Codttd 

S«gi«tt Kli<nt Tnf H— 


Vwdtm 


D«flnlticm 


28 

15 

C 

AKPa 


28 

16 

C 

PCMD 

rad 

28 

17 

C 

AXPR 

1/rad 

28 

18 

C 

AXPEI 

— 

28 

19 

C 

PE 

rad 

28 

20 

C 

DROLL 

- 

28 

21 

C 

TAHCMD 

rad or 


T«d/sec 


28 

22 

C 

DPED 

tad or 
rad/sec 

28 

23 

C 

AKRC 

rad 

28 

24 

C 

AKRGI 

— 

28 

25 

C 

TTRMBDT 

- 

28 

26 

C 

RE 

rad 

28 

27 

C 

AKRE 

1/rad 

28 

28 

C 

AKREI 


28 

29 

C 

DTAW 

* 

28 

30 

C 

AKAT 

rad/ft/ 

sec‘ 

28 

31 

C 

DLNGSTK 

rad or 
rad/sec 

28 

32 

C 

AKQC 

rad 

28 

33 

C 

AKQCI 


28 

34 

C 

PTRMBUT 

- 

28 

35 

C 

QE 

rad 

28 

36 

C 

AKQB 

1/rad 


Roll input integrator selector 
gain (Kj ) 

Shaped in^t to roll control 
ayateai (p^) 

Boll control systea error gain 

(KPe> 

Roll error integrator selector gain 
for roll control system (K^ ) 

®I 

Roll control systaL error (p^ 
Normalised roll control (igoLL^ 

Cockpit input to yaw control system; 
combines lateral stick and pedal 
inputs in accordance with control 
syaten function switching (Rqq^ 

Pedal deflection ( S 

Taw control system Input gain (Kr^) 

Taw input integrator selector 

gain (K^ ) 

Taw control system trim input (1 ) 

'but 

Taw control system error (re) 

Taw control system error pain (K^) 

Yaw error ‘•egrator selector gain 
for yaw cc system (K ) 

Normalised yfM control 

Lateral acceleration feedback gain 

(K. i 

Longitudinal stick deflection 

Pitch control system input pain (Kq^) 

Pitch input integrator selector 

gain (K ) 

Pitch control system trim input 

Pitch control system error (<^) 

Pitch control system error gain 

<*q.> 


A-43 



Coded 


SetWAt 

SlraMt 

.-SzesL 

Haw 

Ualts 

Definition 

as 

If 

c 

MUffL 


Mtch error istogretor eeleetor 
gain for pitch control synten (K 

% 

as 

3S 

c 

onxGB 

— 

HOnelised pitch control ((pxxcH^ 

as 

39 

c 

ZBDTC 

ft/MC 

Beevn control eyetra input 

as 

40 

c 

ZBOTB 

ft/MC 

Heeoe control eyston error (Ip.) 

a 

as 

41 

c 

4KZB 

Ib/ft/ 

•ec 

Hoeve control eyeten error gein 

as 

4a 

c 

4KZBI 

Ib/Mc/ 

ft/sttc 

Beene control eyeten error integntor 

gain 

as 

43 

c 


— 

■nual throttle input 

as 

44 

c 

DSUBF 

— 

Momeltaed aide force control 

as 

'45 

c 



Homallaed nomel force control 

as 

46-47 

c 

PSIIC(l) 

rad 

Gonnanded yaw thruet deflection for 
Ith engine 1*1 for left engine* 

1 " a for right engine (I)) 

as 

48-49 

c 

1HTC(I) 

rad 

Gonnanded pitch thruet deflection for 
1th engine; I 1 fpr left engine* 

I >■ a for right engine (1)) 

as 

50 

c 

DSC 

rad 

Coianded ayneetrical elenon 
deflection ( 

as 

51 

c 

DBC 

rad 

Gonnanded rudder deflection ( 

as 

5a 

c 

04C 

rad 

Conannded differential elewon 
deflactlon 

as 

53 

c 

OKCSR 

— 

Homaliced coanand for RCS roll control 

as 

54 

c 

DBCSP 

— 

Homallaed coneend for RCS pitch control 

as 

55 

c 

D8CSY 

— 

Homallaed connend for RCS yaw control 

as 

56 

c 

Dicssr 

— 

Homallaed connand for RCS side force 
control 

as 

57 

c 

OSCSMF 

~ 

Homallaed coanand for RCS noraal 
force control 

as 

58 

c 

TQfl) 

lb 

Gonnanded thruac (iQfp) 

as 

59 

c 

BILAPC 

rad 

Coananded canard trailing edge flap 
deflection ( ^ 


ArAA 



Coded 

Segment Element Type Name 


Unite 


Definition 


28 

60 

C 

UILAPC 

red 

OoMMaded wing leading edge flap 
deflection 

28 

61 

c 

LELEVHC 

red 

Gonnanded left elevon deflection 

28 

62 

c 

BBLBVHC 

red 

Cn—nded right elevon deflection 
( ^ (2)) 

28 

63 

c 

AKPCIl 

i;e4/sec 

Roll control systan input integrator 
gain (Kp^ ) 

28 

64 

c 

AKPBIl 

1/eec 

liiltlpUee^ll control system 
error gain to form error integrator 

8®^ CKpeu^ 

28 

65 

c 

AKP 

rad/ rad/ 
sec 

Roll control system roll rate feed- 
back gain (Kpp) 

28 

66 

c 

AKRCIl 

rad/sec 

Yaw control system input integrator 

Hultlj^lies yaw control system error 
gain to form error integrator gain 

(Kr^ ) 

®I1 

28 

67 

c 

AKREll 

1/sec 

28 

68 

c 

AKR 

rad/rad/ 

sec 

Yaw control system yaw rate feedback 
gain (Kr^) 

28 

69 

w 

AKQCIl 

rad /sec 

Pitch control syst«n input integrator 

28 

70 

c 

AKQEIl 

1/sec 

Multiplies pitch control system error 
gain to form error Integrator 
gain (Kq ) 

^®I1 

28 

71 

c 

AKQ 

rad/ rad/ 
sec 

Pitch control system pitch rate 
feedback gain (Kq^) 

28 

72 

c 

AKZEIl 

1/sec 

Multiplies heave control system 
error gain to form error Integrator 

28 

73 

N 

csvsw 

— 

Not used in final version of program. 

28 

7U 

c 

TPITCH 

sec 

Pitch input shaping filter time constant 

28 

75 

c 

TROLL 

sec 

Roll input shaping filter time constant 

28 

76 

c 

TYAW 

sec 

Yaw input shaping filter time constant 
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A. 30 Segaant 29 - Klsc«llaneous Air Data* Forca and Ifamt and Klnenatic Variables 


Coded 


CO 

& 

1 

ft 

Elcaent 

Type 

Nane 

Units 

Oeflnltlen 

29 

1 

C 

UAS 

ft /sec 

X body axis coaponent of airspeed 
(“AS> 

29 

2 

C 

VAS 

ft/aec 

Y body axis component of airspeed (v^) 

29 

3 

C 

HAS 

£t/sec 

Z body axis component of airspeed (w^) 

29 

4 

C 


slug/ft^ 

Air density (^) 

29 

5 

C 

IAMB 

deg R 

Aitblent air temperature (Available 
but not used in program) 

29 

6 

C 

VS^D 

£t/sec 

Speed of sound (Vgg) 

29 

7 

C 

TRAIi;0 

— 

Temperature ratio (Available but 
not used In program 

29 

8 

C 

PRATI|I( 

— 

Pressure ratio (Available but not 
used in program) 

29 

9 

C 

VA 

ft/ sec 

Total airspeed (V^) 

29 

10 

C 

VK 

kt 

Total airspeed (V^) 

29 

11 

C 

Ml 


Mach number ()^) 

29 

12 

C 

QBAR 

lb/£t^ 

Dynamic pressure (q) 

29 

13 

c 

ALFE 

rad 

Angle of attack («< ) 

29 

14 

c 

BETA 

rad 

Sideslip angle (^ ) 

29 

15 

C' 

ALFAD 

lad/sec 

Rate of change of angle of attack 
(«) 

29 

16 

c 

SALPH 

— 

sin (oO 

29 

17 

c 

CALPH 

— 

cos (et) 

29 

18 

c 

SBETA 

— 

sin (^ ) 

29 

19 

c 

CBETA 

— 

cos (^ ) 

29 

20 

c 

XAEBU 

lb 

Total X body axis aerodynamic force 
^AERO^ 

29 

21 

c 

YAEM 

lb 

Total Y body axis aerodynamic force 
^AERO^ 

29 

22 

c 

ZAER0 

lb 

Total Z body axis aerodynamic force 
^^AERO^ 

29 

23 

c 

LAER0 

ft lb 

Total aerodynamic body axis roll 
moment (L^^q) 
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Coded 

Element Type Name 


Units 


Definition 


29 

24 

C 

MAEB# 

ft lb 

Total aaiodynaaic body axis pltcti 

moaent (M._.0 
' AERO' 

29 

25 

C 

MAES0 

ft lb 

Total aarodynaalc body axis yav 
Boasnt 

29 

26 

c 

SUKFX 

lb 

Total X body axis applied force 

29 

27 

c 

SUMFf 

lb 

Total Y body axis applied force 

29 

28 

c 

SUMFZ 

lb 

Total Z body axis applied force 

29 

29 

c 

SUML 

ft lb 

Total body axis applied roll 
oonent 

29 

30 

c 

SUMM 

ft lb 

Total body axis applied pitch 
aoment 

29 

31 

c 

SUIM 

ft lb 

Total body axis applied yaw moment 

^tot^ 

29 

32 

c 

DllP 

— 

Corrected d^^ direction cosine (dj^^) 

29 

33 

c 

D12P 

— 

Corrected d^^ direction cosine (djj^) 

29 

34 

c 

D13P 

— 

Corrected d^^ direction cosine (dj^^) 

29 

35 

c 

D21P 

— 

rorrected d^j^ direction cosine (^ 21 ^ 

29 

36 

c 

D22P 

— 

Corrected ^^2 direction cosine i.^ 2 ' 2 } 

29 

37 

c 

D23P 

— 

Corrected d^^ direction cosine 

29 

38 

c 

D31P 

— 

CorT*ected d^^ direction cosine (d^^^) 

29 

39 

c 

D32P 

— 

Corrected d ^2 direction cosine (d^ 2 ) 

29 

40 

c 

D33P 

— 

Corrected d^^ direction cosine (d^^) 

29 

41 

c 

DCHAG 

— 

Determinant of direction coslue 
matrix 

29 

42 

c 

UAIR 

ft/sec 

X body axis component of alrmass 
velocity 

29 

43 

c 

VAIR 

ft/ sec 

Y bod'/ axis component of alrmass 
velocity (v^j^) 

29 

44 

c 

WAIR 

ft/sec 

Z body axis component of alrmass 
velocity (w. j-J 
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Coded 


Seement 

Element 

Type 

Hame 

Units 

Definition 

29 

45 

C 

ALT 

ft 

Altitude ( ■ - *^) 

29 

46 

C 

RS 

rad/ sec 

Stability axis yaw rate (r ) 

29 

47 

C 

PS 

rad/ sec 

Stability aid.8 roll rate (p_) 

29 

48*50 

C 

BL^1C3) 

— 

Not used 

29 

51 

C 

GTI^L 

ft lb 

Gyroscopic body axis roll noment 

(^gyro^ 

29 

52 

C 

GYR0N 

ft lb 

Gyroscopic body pitch laoment 

^**GYR0^ 

29 

53 

C 

GYB^ 

ft lb 

Gyroscopic body yaw moment 

®GYR0^ 

29 

54 

C 

KXCG 

g 

X body axis load factor (<^xcG^ 

29 

55 

C 

. NYCG 

g 

Y body axis load factor (uyqq) 

29 

56 

C 

NZCG 

g 

Z body axis load factor 

29 

57 

C 

XDH 

ft/ sec 

X inertial axis alrmass velocity 
component 

29 

58 

C 

/TDH 

ft/sec 

Y Inertial axis airmails velocity 
component 

29 

59 

C 

ZOU 

ft/sec 

Z inertial axis alrmass velocity 
component 

4 ^ 

60-62 

C 

BI^T2(3) 

— 

Not used 

29 

63 

C 

LFRIME 

rad/ sec^ 

Body axis roll moment adjusted to 
Include cross product of Inertia 
coupling of yaw moment. LPRIME 
equals PDOT at time zero with zero 
initial conditions on angular rates 

29 

64 

C 

NPRIME 

2 

rad/sec 

Body axis yaw moment adjusted to 
Include cross product of Inertia 
coupling of roll moment. NPRIME 
equals ROOT at time zero with zero 
initial conditions on angular rates 
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Coded 


Segnent 

Element 


Name 

Units 

Definition 

29 

65 

C 

VK 

ft/sec 

Magnitude of inertial speed 

( - 4uz + v^ + 

29 

66 

c 

VBDI 

2 

ft/sec 

Rate of change of magnitude of 
Inertial speed. 

29 

67 

c 

VhEST 

ft/ sec 

Predi^cted magnitude of inertial 
speed (VEEST • VEDT*DT + VE) 

A. 31 SegMsnt 30 - 

Jalet 

Ram Forces 

and Moments 

Tables 

Segment 

Element 

Typ* 

Coded 

Name 

Units 

Definition 

30 

1-5 

N 

V/vtT(I) 

— 

Table entries for independent 
variable VOVI used to detetnmlne 
ram drag parameters (1 Is defined 
by AVOVI and must be ^ 5) 

30 

e-13 

M 

. ATUFNT(J) 

rad 

Table entries for Independent 
A 9r.j.able ATUBN used to determine 

inlet ram parameters (J is defined 

by AATURN and must be it 8/ 

30 

14-53 

N 

ELRAKr(I,J)— 

Table entries for normalized (to 
equivalent inlet diameter) inlet 
ram uoment arm. Table Is addressed 
by defining values for VOVI and 
ATURN. ELRAMT(I,J) Is the value 
of ELRAH when VOVI - VOVIT(I) and 
ATURN ATURNT(J) 

30 

54-93 

N 

DATURNT(l.J) rad 

Table entries for difference between 
actual and geometric Inlet flow 
turning angles. Table Is addressed 
by defining values for VOVI and 
ATURN. DAXURNT(I.J) is t)-. value 
of DATURN when VOVI “ VOVIT(I) and 
ATURN - ATURNT(J) 

30 

94-133 

N 

ETARI(I.J) — 

Table entries for ratio of actual 


to theoretical IrJet ram force »nap:ni- 
tude. Table is addressed by defining 

values for VOVI and AIUKN. ETART 
(I,J) Is the value of ETAR when 
VOVI - VOVIT(I) and ATURN - 
ATURNT(J) 
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A.32 8( 


iC 31 « 1C8 Fotcm and HoMats TablM 


Sccaanc Elaaeat 


Coded 

Type Sam'S 


Jfelta 


Dafinltlon 


31 


l-« 


FIC8T(I) lb btrlaa for KCS forca tabla. 

Trbla la uaad vltb BIMT to 
dataralaa BCS forca aa a 
fuaetiott of blaad or vlea varaa. 
I la dafinad by ABIiMl aad auat 
ba « 6. 


31 


7-12 


II BlDBTd) 


Eatrlaa for blaad tabla. Tabla 
la uaad with FBCST to dataralaa 
blaad aa a function of BCS forca 
or vice aarsa. I la dafinad by 
ABLOR and mat ba ^ 6. BLDRT 
(I) aqiuala blaad whan ICS forca 
a<inala TBCST(I). 


A. 33 Sagnant 32 - Thraot Ralatad Tables 

Coded 


SaanAt 

tl«BM&t 


Bane 

Units 

Dafinition 

32 

1-4 

— 1 — 

N 

MHT(I) 


bitrlea In Kacb unaber table. 
Tabla used with TPCattX, TPGHni, 
TF6IDL and 1MDT to obtain 
■artafcw and nlninua afterburning 
thrust, idle thrust and naxlana 
Inlet asss flow rate functions 

of Mach neaber. I Is defined by 
AMdCE and mst be ^ 4. 

32 

5-8 

H 

TFGHAX(I) 

lb 

Entries in aaainnn afterburning 
thrust table. I amst be £ 4. 
TFGMAX(I) equals FGMAX at MK« 
MHT(I) . 

32 

9^12 

N 

TFQflN(I) 

lb 

Entries In mlniaua afterburning 
thrust Cable. 1 nust be It 4, 
TFGMIN(I) equals FOfIN at MM? 
MHT(I). 

32 

13-16 

N 

TPGH>L(I) 

lb 

Entries in Idle thrust table. I 
oust be £ 4. TFGIDL(I) equals 
FGIDL at MH-MNT(I) . 

32 

17-20 

N 

TMDT(I) 

Ibm/sac 

Entries in naxlnutn inlet nass 
flow rate table. I oust be £ 4. 
TMDT(I) equals MDTHX at MN-MNT 
(I). 



Ufttf 


DeCinltiofi 


Scewent Element 


Coded 

Im Sss-. 


32 21-26 

32 27-32 

32 33r42 

32 43-52 

32 53-62 

32 63-72 


■ TFIdld) 


II nnf(i) 

II TFIACTd) 


latrlM In fractinonl nen-nfter- 
baralat thruat table. Table 
vaed with HPN to obtain fractional 
eatina rpn. I ia dafiaad by AFIAT 
and anat be # 6. 

Bntriaa In fractional rpn table. 

I anat be ^ 6. nPM(I) aquala 
FIPM at mni^TFUTd). 

btriea in fractional non-after- 
bomlng tbxuat table. Table 
need with TAOC&L, TMCEL, and 
TTERG to obtain engiaa accalaration 
and dacaleratioa llnita and tine 
conatant. I ia defined by AFBAC 
and auat be < 10. 


II T40CEL(I) l/aoc Entriaa in acceleration Unit 

table. I auat be 4 10. TACCEL 
(1) aqoala acceleration linit at 
TmOTFKACTd) . 

M TSBCEL(t) 1/aoc Bntriaa in dacaleratloc. linit 

table. I auat be ^ 10. TDEOn. 
(I) eqinala deceleration linit at 
TFRAOTFRACT(I) . 


N TTBNU(I) see Bntriea in engine tine constant 

table. I auat be 4 10. TTENG 
(I) equals TENG at TFRAOTPRACT(I) . 


A. 34 Segnents 33 tfarfiugh 39 - Stability Derivative Arrays and Control Parameters 

X. T and Z force and roll, pitch and yaw aonent stability derivatives 
are stored in that order in segnents 33 through 38. l.e. . X force derivatives 
are contained In segaent 33 .Y force derivatives In segment 34. etc. to yaw 
aonent derivatives in segaent 38. Stored within each segaent are the total 
derivatives along with the contributions to those derivatives froa each source 
of applied forces and aonants. 

Derivative contributions froa wing, horisontal stabiliser, vertical 
stabiliser, fuselage, total aero, direct thrust, ran drag. Coriolis, and RCS- 
are included. Derivatives are calculated for perturbations In 22 (space for 
25 la available) state and/or control variables. These variables designeted 
by I In the segaent listing below are as follows: 
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1 

Perturbed Varlabxe 

Seaeent 

Kloient 

Units 

1 

VAS 

29 

1 

ft/sec 

2 

WAS 

29 

3 

ft/aec 

3 

Q 

3 

5 

rad/sec 

4 

OBC 

28 

50 

rad 

5 

DFLAPW 

3 

49 

rad 

6 

DFLAFR 

3 

48 

rad 

7 

DRCSP 

28 

54 

- 

8 

DRCSNF 

28 

57 

- 

9 

1)0(1) 

26 

34 

lb 

10 

m2) 

26 

35 

lb 

11 

TST(l) 

3 

26 

rad 

12 

THT(2) 

3 

27 

rad 

13 

VAS 

29 

2 

ft/aec 

14 

P 

3 

4 

rad /sec 

15 

R 

3 

6 

rad/aec 

16 

Ml 

3 

21 

rad 

17 

DAC 

28 

52 

rad 

18 

MICSR 

28 

53 

* 

19 

DRCST 

28 

55 

- 

20 

DRCSSP 

28 

56 


21 

psrr(i) 

3 

24 

rad 

22 

PS1T(2) 

3 

25 

rad 


Note In the segment listing that PDPIIV ties adopted as a shorthand notation for 
"per unit perturbation In Ith variable" 

Coded 


SeRuent 

Eleeent 


Naee 

Units 


Definition 

33 

1-25 

C 

DXMG(I) 

ft/sec^/ 
units of 

I 

Change in wing X body axis force 
PUPIIV (per unit perturbation in 
Ith variable) . 

33 

26-50 

C 

DXHS(I) 

ft/sec^/ 
units of 

I 

Change in horizontal stabilizer 
X body axis force PUPIIV. 

33 

51-75 

C 

DXVS(I) 

ft/sec^/ 
units of 

I 

Change in vertical stabiliser 
X body axis force PITPIIV. 

33 

76-100 

C 

DXFS(I) 

ft/sec^/ 
uults of 

I 

Change in fuselage X body axis 
force PUPIIV. 

33 

101-125 

C 

DXAER(I) 

ft/sec^/ 
units of 

I 

Change in aerodynamic X body 
axis force PUPIIV. 

33 

126-150 

C 

DXTS(I) 

ft/sec^/ 
units of 

I 

Change in direct thrust X body 
axis for^i PUPIIV. 

33 

151-175 

C 

DRRAHd) 

ft/sec^/ 
units of 

I 

Change in ijilet rrua X body axis 

force PUPIIV. 
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Coded 


SccBksnt 

EleEicnt 

Type 

Maine 

Units 

Def inition 

33 

176-200 

c 

nzcLs(i) 

ft/s«c^/ 
units of I 

Change Coriolis X body 

axis force PUPIIV. 

33 

• 

201-225 

c 

UKCSFd) ft/MC^/ 

imlts of I 

Change in RCS X body axis 
force PDPIIV. 

33 

226-250 

c 

X(I) 

ft/s«c^/ 
units of I 

Change in total X body axis 
force PUPIIV. 

34 

1-25 

c 

DTUGd) 

ft/sec^/ 
units of I 

Change In vlng Y body axis 
force per PUPIIV. 

34 

26-50 

c 

DTHSd) 

ft/sec^/ 
units of I 

Change in horisontsl stsbili».ur 
T body axis force PUPIIV. 

34 

51-75 

c 

DTVS(I) 

ft/sec^/ 
units of I 

Change in vertical stabiliser 
T body axis force PUPIIV. 

34 

76-100 

c 

DTFS(I) 

ft/sec^/ 
units of I 

Change in fuselage T body axis 
force PUPIIV. 

34 

101-125 

c 

DTAER(I) 

ft/sec^/ 
units of I 

Change in serodyniAic T body 
axis force PUPIIV. 

34 

126-150 

c 

DTTSd) 

ft/sec^/ 
units of I 

Change in direct thrust T body 
axis force PUPIIV. 

34 

151-175 

c 

DTSAMd) 

ft/sec^/ 
units of I 

Change in inlet ram Y body axis 

force PUPIIV. 

34 

176-200 

c 

DTCLSd) 

ft/sec^/ 
units of I 

Change in Coriolis Y body axis 
force PUPIIV. 

34 

201-225 

c 

.DTRCSF(I) 

ft/sec^/ 
units of I 

Change in RCS Y body axis force 
PUPIIV. 

34 

226-250 

c 

Td) 

ft/sec^/ 
urits of I 

Change in total Y body axis 
force PUPIIV. 

35 

1-25 

c 

DZWG(I) 

ft/sec*"/ 
units of I 

Change I’t wing Z body axis 
force PUPIIV. 

35 

26-50 

c 

D2HS(I) 

ft/sec^/ 
units cf 1 

Change in horisontal stabilizer 
Z body axis force PUPIIV. 

35 

51-75 

c 

DZYSd) 

ft/sec^/ 
units of I 

Change in vertical stabilizer Z 
body axis force PUPIIV. 

35 

76-100 

c 

DZFS(I) 

ft/sec^/ 
units of I 

Change in fuselage Z body axis 
force PUPIIV. 
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Coded 


Se spent 

Elcnenc 

Type 

Noae 

Uhits 

Definition 

33 

101-125 

C 

DZAER(I) 

ft/eec^/ 
units of I 

Change in ssrodynaaic Z body 
axis fores PUPIIV. 

35 

• 

125-150 

C 

DZTS(I) 

ft/sec^/ 
units of I 

Change in direct thrust Z 
body axis force PUPIIT. 

35 

151-175 

C 

DZXAM(I) 

ft/ssc^/ 
units of I 

Change in inlet ram Z body 

axis force FUPIIV. 

35 

176-200 

c 

DZCLS(I) 

ft/ssc^/ 
caits of I 

Change in Coriolis Z body 
axis force PUPIIV. 

35 

201-225 

c 

DZRCSF(I) ft/MC^/ 

unite of I 

Change In RCS Z body axis 
force ^tTIIV. 

35 

226-250 

c 

Z(I) 

ft/sec^/ 
units of I 

Change in total Z body axis 
force PUPIIV. 

36 

1-25 

c 

DU)G(I) 

rsu/sec^/ 
units of I 

Change in wing body axis 
roll nonent PUPIIV. 

36 

26-50 

c 

DLHS(I) 

rsd/sec^/ 
units of I 

Change in horlsontal stabiliser 
body axis roll nonent PUPIIV. 

36 

51-75 

c 

DLVS(I) 

rsd/sec^/ 
units of I 

Change in vertical stabiliser 
body axis roll nonent PUPIIV. 

36 

76-100 

c 

DLFS(I) 

rsd/sec^/ 
units of I 

Change in fuselage body axis 
roll nonent PUPIIV. 

36 

101-125 

c 

DLAER(I) 

rsd/sec^/ 
units of I 

Change in aarodynanic body 
axis roll nonent PUPIIV. 

36 

126-150 

c 

.DLTS(I) 

rad/sec^/ 
units of I 

Change In direct thrust body 
axis roll nonent PUPIIV. 

36 

151-175 

c 

DLRAM(I) 

rad/sec^/ 
units of I 

Change In inlet ram body axis 

roll nonent PUPIIV. 

36 

176-200 

c 

DLCLS(I) 

rad/sac^/ 
«’nits of I 

Change in Coriolis body axis 
roll nonent PUPIIV. 

36 

201-225 

c 

DLRCSF(I) 

rad/sec^/ 
units of I 

Change in RCS body axis roll 
nonent PUPIIV. 

36 

226-250 

c 

AL(I) 

rad/sec^/ 
units of I 

Change in total body axis roll 
nonent PUPIIV. 


A-54 


Coded 



Elcaent 

Type 

Nome 

Oalcs 

Definition 

37 

1-25 

C 

DMIG(I) 

rad/sec^/ 
unite of I 

Change in ving body axis pitch 
Mcnt PUPIIV. 

37 

• 

26-50 

C 

DI«8(I) 

red/eec^/ 
units of I 

Change in horisontal stabiliser 
body axis pitch nonent POFIIV. 

37 

51-75 

c 

SMVS(I) 

rad/sec^/ 
units of I 

Change in vertical stabiliser 
body axis pitch nonent PUPIIV. 

37 

76-100 

c 

OMFS(I) 

rad/sec^/ 
units of I 

Change in fuselage body axis 
pitch nonent PDPICV. 

37 

101-125 

c 

OHAERd) 

rad/sac^/ 
units of I 

Change in aerodynanic body axis 
pitch nonent PUPIIV. 

37 

126-150 

c 

OHTS(I) 

rad/sec^/ 
units of I 

Change in direct thrust body 
axis pitch nonent PUPIIV. 

37 

151-175 

c 

OHRAM(I) 

rad/sec^/ 
u"tts of I 

Change in inlet ram body axis 

pitch nonent PUPIIV. 

37 

176-200 

c 

DHCLS(I) 

rad/sec^/ 
units of I 

Change in Coriolis body axis 
pitch nonent PUPIIV. 

37 

201-225 

c 

IMRCSF(I) 

rad/sec^/ 
units of I 

Change in RCS body axis pitch 
nonent PUPIIV. 

37 

226-250 

c 

AM(I) 

rad/sec^/ 
units of I 

Change In total body axis pitch 
nonent PUPIIV. 

38 

1-25 

c 

DNWG(I) 

rad/ sec^/ 
units of I 

Change in wing body axis yaw 
nonent PUPIIV. 

38 

26-50 

c 

mmsd) 

rad/sec^/ 
units of I 

Change in horisontal stabiliser 
body axis yaw nonent PUPIIV. 

38 

51-75 

c 

DNVSd) 

rad/sec^/ 
units of I 

Change in vertical stabilizer 
body axis yaw moment PUPIIV. 

38 

76-100 

c 

DNFSd) 

rad/sec^/ 
units of I 

Change in fuselage body axis 
yaw moment PUPIIV. 

38 

101-125 

c 

DNAER(I) 

rad/sec^/ 
units of I 

Change in aerodynamic body axis 
yaw moment PUPIIV. 

38 

126-150 

c 

DNTSd) 

rad/sec^/ 
units of I 

Change in direct thrust body axis 
yaw moment PUPIIV. 

38 

151-175 

c 

DNIlAMd) 

2 

rad/sec / 
units of I 

Change in inlet ram body axis yaw 

moment PUPIIV. 

38 

176-200 

c 

DNCLS(I) 

rad/sec^/ 
units of I 

Change in Coriolis body axis yaw 
moment PUPIIV. 
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Coded 


Sceaent 

Element 

Type 

Nome 

Units 

Definition 

38 

201-225 

C 

DNRCSFd) 

rad/sec^/ 
units of 1 

Change in ECS body axis yaw 
moment PUPIIV. 

38 

• 

226-250 

C 

AH(I) 

rad/see^/ 
units of I 

Change in total body axis yaw 
moment PUPIIV. 

SogMnt 39 ttoros infoiaatlon rogardlng the nuaber of verlablee to be 
perturbed to fore atebUlty derivatives and the else of these perturbations. 

Seseent 

Element 

Tyte 

Ceded 

Name 

Units 

Definition 

39 

1 

M 

ADER 


Number of variables to be 
perturbed to form stability 
derivatives ( < 25) . In 
current program set-up, ADER* 
22. 

39 

2-26 

H 

DELV(I) 

Varies 

Slse of perturbation in Ith 
variable to be used in com- 
puting stability derivatives. 

A. 35 Segment 40 - 

Scratch COMMON 



Fifty locations labeled WORK have been set aside as segment 40 for 
use In temporary program changes or for further program development. 

The WORK array la not used in the current program set-up. 

A. 36 Segment 41 - 

Miscellaneous Constants 


Seament 

Element 

IZEl. 

Coded 

Name 

Units 

Definition 


41 

1 

N 

AWIND 


AWIND not equal to zero indi- 
cates a moving alrmass and 
that the trim in a wind option 
is to be invoked. 

41 

2 

N 

VWIND 

tt/sec 

Magnitude of alrmass velocity 
(VWIND) 

41 

3 

N 

PSIWIND 

rad 

Direction of alrmass velocity. 
For headwind, PSIWIND-180 deg* 
3.14159 rad. «VlND) 

41 

4 

N 

KE2D 

ft/aec^ 

Initial value of aircraft 
inertial acceleration along 
X inertial axis 
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Coded 


Secaent 

Clectcnt 

Type 

ilome 

Uhlts 

Definition 

41 

5 

M 

TE2D 

ft/sac^ 

Initial value of aircraft 
inertial acceleration along T 
inertial axle 

41 

• 

6 

H 

EE2D 

ft/sec* 

Initial value of aircraft 
Inertial acceleration along Z 
Inertial axle 

1*1 

7 

N 

TRMTURN 


TRMTURK not equal to zero indi- 
cates that aix plane is to be 
trimmed in a steady turn. GAMMA 
VADOT, and TURNRAD must be 
specified if TRMTURN f^O. 

l»l 

8 

N 

GAMMA 

rad 

Airplane flight path angle 

1*1 

9 

N 

VADOT 

ft/sec^ 

Airplar.e deceleration along 
flight path 

1*1 

10 

N 

TURNRAD 

ft 

Radius of turn 


A. 37 Segment 42 - Aero Surface Effectiveness Tables 


Coded 

Sewteot Blenent Type Nsae Units 

42 1-7 N ELEVT(I) rad 


42 8-14 N EKELET(I) 


42 15-21 N RUDT(I) red 


42 22-28 N EKRUDT(I) - 


Definition 


Entries In eleven table. Table 
used %rith EKELET to obtain eleven 
effectiveness as a function of elevoi 
deflection. I is definr^d by AEKEL 
snd oust be ^7. 

Entries in eleven effectiveness 
table. I niust be ^ 7 . EKELET 
(I) equals elevon effectiveness 
at ABS (LELEVN) or ABS (RELEVN)- 
(ELEVT(I) . 

Entries in rudder table. Table 
used with EKRUDT to obtain rudder 
effectiveness as a function of 
rudder deflection. I is defined 
by AERRUD and vust be & 7. 

Entries in rudder effectiveness 
table. I Bust be ^ 7. EKRUDT 
(I) equals rudder effectiveness 
at ABS(DR)«RUDT(I) . 


originai. pace b 

OF POOR 
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A. 38 Segment 43 - Peeudo-Pilot Input Control Data 


CcKled 


Segment 

Element 


Name 

Units 

Definition 

43 

1-5 

N 

ACMD(I) 


Each nonzero ACMD(l) specifies a 
cockpit control input time history 
vhich is to be imposed on the air- 
craft. ACMD(I) “ 1. indicates 
DLNGSTK input, « 2. indicates 
DLATSTiC input, ® 3. indicates 
DPED input, =4. indicates ZEDTC 
input, =5. indicates THR0T input. 
For example, if ACMD(l) = 3. and 
ACMD(2) « 1., simultaneous DPED 
and DLNGSTK inputs are imixased on 
the aircraft. 

43 

6 

N 

ATIME 

—— 

Number of entries in the cockpit 
control time history tables and 
must be ^ 20 

43 

7 

N 

ASKEEP 

— 

ASKEEP not equal to zero activates 
pseudo-pilot stationkeeping loops. 

43 

8 

N 

AYEHLD 

— 

Not used in final versi n of program* 

43 

9-28 

N 

CMDTI(I) 

varies 

Entries in DLNGSTK input time 
history table. CMDTI(I) equals 
DLNGSTK when T(time) = TIMET(l). 

43 

29-48 

N 

CMDT2(I) 

varies 

Entries in DLATSTK input time 
history table. CMDT2(I) equals 
DLATSTK when T(time) =TIMET(I). 

43 

49-68 

N 

CMDT3(I) 

varies 

Entries in DPED input time history 
table. CMDT3(l) equals DPED when 
T(time) =■ TIMET(I) 

43 

69-88 

N 

CMDT4(I) 

ft/sec 

Entries in ZEDTC input time history 
table. CMDT4(I) equals ZEDTC when 
T(time) =TIMET(T) 

43 

89-108 

N 

CMDT5(I) 

fractional 

throttle 

deflection 

Entries in THROT input time history 
table. CMDT5(I) equals THROT 
when T(time) STIMET(I). 

43 

109-128 

N 

TIMET(I) 

sec 

Entries in time table. Table used 
with CMDTl, CMDT2, CMDT3, CMDT4, and 
CMDT 5 . I is defined by ATIME and 
must be £ 20 . 

43 

129 

N 

ATRANS 

— 

ATRANS not equal to zero activates 
a pseudo-pilot flown transition. 

43 

130 

N 

AVAT 


Number of entries in transition 
tables (THROTT, VAT, and THETCT). 
Must be ^10. 

43 

131 

N 

TTRANS 

sec 

Time at which transition is initiated 
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Coded 

Segment Element Type Name- Units Definition 


43 132-141 

43 142-151 

43 152-161 


N THROTT(I) fractional Entries in transition throttle table. 

throttle THROfT(l) equals manual throttle 
deflection setting when VE or VEEST »VAT( I) 


N VAT(I) ft/sec Entries in transition speed table. 

Table used with THROTT and THETCT 
to obtain throttle setting auid pitch 
angle as a function of Inertial speed 
during transition. 


N THEPCT(I) rad Entries in transition pitch angle 

table. THETCT(I) equals desired 
pitch angle when VE or VEEST = 
VAT(I) 


A. 39 Segment 44 - Tables, Constants, and Variables for Varying Feedback Control Laws 

Coded 


Segment 

Element 

Type 

Name 

Units 

Definition 

44 

1 

N 

ACSSW 


Number of entries in control system 
variable selector tables. Must be 

44 

2-6 

N 

Tcssva(i) 


Entries in CSSW table. Table en- 
tered with CSSW and usjd with all 
other tables in this COMMON segment. 

44 

7-11 

N 

TAKQS(I) 


Entries in pitch rate selector table, 
TAKQS(I)?<0. specifies that Q is fed- 
back in pitch control system. 

44 

12-16 

N 

TKPINT(I) 


Entries in integral of roll rate 
selector table. TKPINT(I)?^0. speci- 
fies that PINT is fedback in roll 
control system when CSSW=TCSSWl(l) . 

44 

22-26 

N 

TKPB(I) 


Entries in roll rate selector table 
TKPB(I) ?^0. specifies P is fedback 
in roll control system when CSSW= 
TCSSWl(I). TKPB(I)*0. specifies PS 
feedback when CSSW=TCSSWl(l) . 

44 

27-31 

N 

TKPS{I) 

— 

Not used in final version of program 

44 

32-36 

N 

TKQINT(I) 


Entries in integral of pitch rate 


selector table. TKQINT^^O. specifies 
that QINT is fedback in pitch control 
system when CSSW=Tf'SSWl(l) . 
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Coded 


Segment 

Element 

Type 

Name 

Units 

Definition 

UU 

37-41 

K 

TKTHETA(I) 

— 

Entries in pitch angle selector 
table. TKTHETA(I)?‘ 0 . specifies 
that THETA is fedback in pitch con- 
trol system '«dien CSSW=TCSSW 1 (I). 

UU 

42-46 

N 

TKRINT(I) 


Entries in integral of yaw rate 
selector table. TKRINT^^O. specifies 
that HINT is fedback in yaw control 
system when CSSW=TCSSW 1 (I) . 

UU 

47-51 

H 

TKPSI(I) 


Entries in yaw angle selector table. 
TKPSI(I)# 0 . specifies that PSI is 
fedback in yaw control system when 
CSSW=TCSSW 1 (I). 

4U 

52-56 

N 

TKRB(I) 


Entries in yaw rate selector table. 
TKRB(I)?^ 0 . specifies that R is fed- 
back in yaw control system when 
CSSW=TCSSW(i). TKRB(I)= 0 . specifies 
RS feedback when CSSW=TCSSWl(l). 

44 

57-61 

N 

TAKPCI(I) 

- 

Entries in roll input integrator 
selector table. TAKPCl(l)?* 0 . 
activates roll input integrator 
when CS 3 W=TCSSW 1 (I). 

44 

62-66 

N 

TAKPEI(I) 


Entries in roll error integrator 
selector table. TAKPEl(l)^^ 0 . 

Activltates roll error integration 
in rolL control system whan CMW 
TCSSWl(I) 

44 

67-71 

N 

TAKRCI(I) 


Entries in yaw input integrator se- 
lector table. TAKRCl(l)?* 0 . activates 
yaw input integrator when CSSW=TCSSWl(l) 

44 

72-76 

N 

TAKPEI(I) 


Entries in yaw error integrator 
selector table. TAKREl(l )?^0 activates 
yaw error integrator in yaw control 
system when CSSW=TCSSWl(l) . 

44 

77-81 

N 

T.\KQCI(I) 


Entries in pitch input integrator 
selector table. TAKQCl(l)# 0 . activates 
pitch input integrator when CSSW= 
TCSSWI(I). 

44 

82-86 

N 

TAKQEI(I) 


Entries in pitch error integrator 
selector table. TAKQEl(l)j^ 0 . acti- 
vates pitch error integrator in pitch 
control system when CSSW=TCSSVn ( I ) . 

44 

87-91 

N 

TOS(T) 

— 

Not used in fin^l version of program. 

44 

92 

N 

AKPHIRS 

varies 

AKPHIRS?« 0 . specifies (GC/VA)* PHI 


feedback in yaw control system (K^ ). 

•r 

s 
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CtHlfVl 

Segment Element Type Name Units Definition 


kk 

93 

C 

AKPINT 


Integral of roll rate feedback 
selector gain for roll control syston 
(Kj, ) 

kk 

9k 

C 

AKPHI 

— 

Roll angle feedback selector gain 
for roll control system, (Kjj) 

kk 

95 

c 

AKPB 


Roll rate feedback selector gain for 
roll control system. AKPB=0. selects 
PSj AKPB*1, selects P; AKPBii^O. or 1. 
mixes P and PS. (K^^) 

kk 

96 

c 

AKPS 

— 

Not used in final version of program. 

kk 

97 

c 

AKQINT 


Integral of pitch rate feedback 
selector gain for pitch control 
system (Kj^^) 

kk 

98 

c 

AKTHETA 


Pitch angle feedback selector gain 
for pitch control system 

Ul* 

99 

c 

AKEINT 


Integrsd of yaw rate feedback se- 
lector gain for yaw control system 
(Ksr). 

uu 

100 

c 

AKPSI 

— 

Yaw angle feedback selector gain for 
yaw control system (Ky/) 

kk 

101 

c 

AKPB 


Yaw rate feedback selector gain 
for yaw control system. AKRB=0. se- 
lects RS; AKRB*1. selects R; 

AKRBf^O. or 1. mixes R and RS. (K- ) 

kk 

102 

c 

AKRS 

— 

Not used in Cinal version of program 

uu 

103 

c 

AKQ3 

— 

Pitch rate feedback selector gain for 
pitch control system. 

A.UO 

Segment 1*5 

- Tables emd Constants for 
Functio*' i of Airspeed 

Control System Gains which Vary as 

Segment 

Element 


CrKlod 

Name 

Units 

Definition 

*♦5 

1 

N 

VELGATN 


VELGAINf*0. specifies that control 


system gains are functions of air- 
speed. VELGAIN=0. implies that 
control system gains are functions 
of pitch angle 
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Coded 

Segment Element Type Name- Units Definition 


U5 

2 

N 

AVEL 

” 

Number of entries in control system 
gains airspeed table (^). 


3-10 

N 

TVEL(I) 

ft /sec 

Airspeed table for determining con- 
trol system gains 

1*5 

11-18 

N 

VAKPC(I) 

rad 

Roll control -ystem input gain table 
if the gain is a function of air- 
speed. 

1*5 

19-26 

N 

VKPCII(I) 

rad/sec 

Roll control system input integrator 
gain table if the gain is a function 
of airspeed. 

1*5 

27-3' 

N 

’/KPEII(I) 

rad/sec/ 

rad 

Roll control system error integrator 
gain table if the gain is a function 
of airspesd. 

1*5 

35-1*2 

N 

VKP(I) 

red/rad/ 

sec 

Roll control system roll rate feed- 
back gain table if the gain is a 
function of airspeed 

1*5 

1*3-50 

N 

VAKPE(I) 

1/rad 

Roll control system error gain table 
if the gain is a function of air- 
speed 

1*5 

51-58 

N 

VAKRC(I) 

rad 

Yaw control system input gain table 
if the gain is a function of air- 
speed 

1*5 

59-66 

N 

VKRCII(I) 

rad/ sec 

Yaw control system input integrator 
gain table if the gain is a function 
of airspeed. 

1*5 

67 - 71 * 

N 

VKREII(I) 

rad/sec/ 

rad 

Yaw control system error integrator 
gain table if the gain is a function 
of airspeed 

1*5 

75-82 

N 

VKR(I) 

rad /red/ 
sec 

Yaw control system yaw rate feedback 
gain table if the gain is a function 
of airspeed. 

1*5 

83-90 

N 

VAKRE(I) 

1/rad 

Yaw control system error gain table 
if the gain is a function of airspeed 

1*5 

91-98 

N 

VAKQC(I) 

rad 

Pitch control system input gain table 
if the gain is a function of airspeed 

1*5 

99-106 

N 

VKQCII(I) 

rad/sec 

Pitch control system input integrator 
gain table if the gain is a function 
of airspeed. 
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Coded 


Segment 

Element 


Raaie . 

Unite 

Definition 


107-llU 

N 

VKQEII(I) 

rad/sec/ 

Pitch control system error inte- 
grator gain table if the gain is 
a function of airsi>eed. 

h3 

115-122 

N 

VKQ(I) 

rad/rad/ 

see 

Pitch control system pitch rate 
feedback gain table if the gain 
is a function of airspeed. 


123-130 

R 

VAKQfi(l) 

1/rad 

Pitch control system error gain 
table if the gain is a function 
of airspeed 


131-138 

R 

VAKZE(I) 

Ib/ft/sec Heave control system error gain table 
if the gain is a function of airspeed 

i*5 

139-11*6 

N 

VKZE'^ld) 

It/sec/ 
ft /sec 

Heave control system error integrator 
gain table if the gain is a function 
of airspeed 

i*5 

ll»7-15l* 

N 

VTPITCH(I) 

sec 

Pitch input shaping filter time con- 
stant table if the time constant is 
a function of airspeed 

k5 

155-162 

N 

VTROLL(I) 

sec 

Roll input shaping filter time con- 
stant table if the time coistant is 
a function of airspeed 

h5 

163-170 

N 

VTYAW(I) 

sec 

Yaw input shaping filter time con- 
stant title if the time constant 


is a :\'.ncL'i’n of airspeed. 
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